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Abstract. We provide a new method for treating free boundary problems in perfect 
fluids, and prove local-in-time well-posedness in Sobolev spaces for the free-surface in- 
compressible 3D Euler equations with or without surface tension for arbitrary initial 
data, and without any irrotationality assumption on the fluid. This is a free bound- 
ary problem for the motion of an incompressible perfect liquid in vacuum, wherein the 
motion of the fluid interacts with the motion of the free-surface at highest-order. 
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1. Introduction 

1.1. The problem statement and background. For ct > and for arbitrary initial data, 
we prove focal existence and uniqueness of solutions in Sobolev spaces to the free boundary 
incompressible Euler equations in vacuum: 



dtu + VuU + Vp = in 


Q, 


(l.la) 


div u = in 


Q, 


(1.1b) 


p = (7 H on 


dQ, 


(1.1c) 


(St + v„)|aoGr(aQ), 




(l.ld) 


u = uo on 


Qt=o J 


(Lie) 


Qt=o = ^ , 




(l.lf) 



where Q = Uo<t<T{t} x n{t), n{t) C R", n = 2 or 3, dQ = [Jo<t<T{t} x dn(t), V„;u = 
u^du^/dx^ , and where Einstein's summation convention is employed. The vector field u is 
the Eulerian or spatial velocity field defined on the time-dependent domain ^{t), p denotes 
the pressure function, H is twice the mean curvature of the boundary of the fluid dn{t), 
and a is the surface tension. Equation (l.la) is the conservation of momentum, (1.1b) is 
the conservation of mass, (1.1c) is the well-known Laplace- Young boundary condition for 
the pressure function, (Lid) states that the free boundary moves with the velocity of the 
fluid, (Lie) specifies the initial velocity, and (l.lf) fixes the initial domain O. 

Almost all prior well-posedness results were focused on irrotational fluids (potential flow), 
wherein the additional constraint curlu = is imposed; with the irrotationality constraint, 
the Euler equations (1.1) reduce to the well-known water-waves equations, wherein the 
motion of the interface is decoupled from the rest of the fluid and is governed by singular 
boundary integrals that arise from the use of complex variables and the equivalence of 
incompressibility and irrotationality with the Cauchy-Riemann equations. For 2D fluids 
(and hence ID interfaces), the earliest local existence results were obtained by Nalimov [14], 
Yosihara [22], and Craig [5] for initial data near equilibrium. Bcalc, Hou, & Lowcngrub [4] 
proved that the linearization of the 2D water wave problem is well-posed if a Taylor sign 
condition is added to the problem formulation, thus preventing Rayleigh- Taylor instabilities. 
Using the Taylor sign condition, Wu [20] proved local existence for the 2D water wave 
problem for arbitrary (sufficiently smooth) initial data. Later Ambrose [2] and Ambrose 
& Masmoudi [3], proved local well-posedness of the 2D water wave problem with surface 
tension on the boundary replacing the Taylor sign condition. 

In 3D, Wu [21] used Clifford analysis to prove local existence of the full water wave 
problem with infinite depth, showing that the Taylor sign condition is always satisfied in the 
irrotational case by virtue of the maximum principle holding for the potential flow. Lannes 
[11] provided a proof for the finite depth case with varying bottom by implementing a Nash- 
Moser iteration. The first well-posedness result for the full Euler equations with zero surface 
tension, cr = 0, is due to Lindblad [13] with the additional "physical condition" that 

S7p-n<0 ondQ, (1.2) 

where n denotes the exterior unit normal to 9f2(f). The condition (1.2) is equivalent to the 
Taylor sign condition, and provided Christodoulou & Lindblad [6] with enough boundary 
regularity to establish a priori estimates for smooth solutions to (1.1) together with (1.2) and 
(7 = 0. (Ebin [10] provided a counterexample to well-posedness when (1.2) is not satisfied.) 
Nevertheless, local existence did not follow in [6], as finding approximations of the Euler 
equations for which existence and uniqueness is known and which retain the transport-type 
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structure of the Euler equations is highly non-trivial, and this geometric transport-type 
structure is crucial for the a priori estimates. In [12], Lindblad proved well-posedness of the 
linearized Euler equations, but the estimates were not sufficient for well-posedness of the 
nonlinear problem. The estimates were improved in [13], wherein Lindblad implemented a 
Nash-Moscr iteration to deal with the manifest loss of regularity in his linearized model and 
thus established the well-posedness result in the case that (1.2) holds and a = 0. 

Local existence for the case of positive surface tension, a > 0, remained open, and 
although the Laplace- Young condition (1.1c) provides improved regularity for the boundary, 
the required nonlinear estimates are more difficult to close due to the complexity of the 
mean curvature operator, and the need to study time-differentiated problems, which do 
not arise in the a = case. It appears that the use of the time-differentiated problem in 
Lindblad's paper [13] is due to the use of certain tangential projection operators, but this is 
not necessary. We note that our energy function is different from that in [13], and provides 
better control of the Lagrangian coordinate. 

After completing this work, we were informed of the paper of Schweizer [16] who studies 
the Euler equations for u > in the case that the free-surface is a graph over the two-torus. 
In that paper, he obtains a priori estimates under a smallncss assumption for the initial 
surface; well-posedness follows under the additional assumption that there is no vorticity on 
the boundary. We also learned of the paper by Shatah and Zeng [15] who establish a priori 
estimates for both the a = and a > cases without any restrictions on the initial data. 



1.2. Main results. We prove two main theorems concerning the well-posedness of (1.1). 
The first theorem, for the case of positive surface tension cr > 0, is new; for our second 

theorem, corresponding to the zero surface tension case, wc present a new proof that does 
not require a Nash-Moser procedure, and has optimal regularity. 

Theorem 1.1 (Well-posedness with surface tension). Suppose that a > 0,r is of class H^-^, 
and uq e H^-^{Q). Then, there exists T > 0, and a solution (u{t),p{t),fl{t)) of (1.1) with 
u e L~(0,T;i?4 -^(fXi))); P e L^{0,T-H^{n{t))), andT{t) e H^-^ . The solution is unique 
ifuoeH^-^ andVeH^-^. 

Theorem 1.2 (Well-posedness with Taylor sign condition). Suppose that a — 0, d^l is 
of class , and uq € H^{fl) and condition (1.2) holds at t = 0. Then, there exists 
T > {), and a unique solution (u(t),p{t),9.{t)) of (1.1) with u € L'^{0,T- H^{n(t))), p e 
L°^{0,T;H^ '%n{t))), and dn(t) G H^. 

1.3. Lagrangian representation of the Euler equations. The Eulerian problem (1.1), 
set on the moving domain ^l{t), is converted to a PDE on the fixed domain fi, by the use 
of Lagrangian variables. Let r]{-,t) : Cl ^ Q{t) he the solution of 

dtr]{x,t) = u{r]{x,t),t), r]{x,0)=ld. 

and set 



v{x,t) := u{r]{x,t),t), q{x,t) := p{ri{x,t),t), and a{x,t) := [Vri{x,t)] 
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The variables v, q and a are functions of the fixed domain f2 and denote the material velocity, 
pressure, and pull-back, respectively. Thus, on the fixed domain, (1.1) transforms to 

?7 = Id+ / V mnx{0,T], (1.3a) 
Jo 

dtv + a\/q = in O x (0, T] , (1.3b) 

Tr(a Vt;) =0 in 17 x (0, T] , (1.3c) 

qa^N/\a^N\= -aAg{r]) onrx(0,T], (1.3d) 

(?7, v) = (Id, uo) on O X {t = 0} , (1.3e) 

where denotes the unit normal to T, and Ag is the surface Laplacian with respect to the 
induced metric g on F, written in local coordinates as 

= V9~^da[y/gg°''^dfj] , g"'^ = [fiia/s]"^ , 9ai3 = • V.is , and ^ = y/detg. (1.4) 

Theorem 1.3 (cr > 0). Suppose that a > 0, dil is of class H^-^ , and e H'^-^{Vt) 
with divwo = 0. Then, there exists T > 0, and a solution (v,q) of (1-3) with v G 
L°°(0,T;iJ4-5(fi)), q e L°°{0,T; H%n)), and T{t) e H^-^. The solution satisfies 

sup (|5n(t)|L + Ell^^Wlli5-i.5. + El|5^(i)llti.5J <A^o 

where Mq denotes a polynomial function of ||r||5.5 and ||wo||4.5- The solution is unique of 
Uo e H^-^{n) and T e iJ^-^ 

Remark 1 . Our theorem, is stated for a fl,uid in vacuum,, hut the analogous theorem holds 
for a vortex sheet, i.e., for the motion of the interface separating two inviscid immiscible 
incom,pressible fluids; the boundary condition (1.1c) is replaced by \p]± = aH, where \p]± 
denotes the jump in pressure across the interface. 

For the zero-surface-tension case, we have 

Theorem 1.4 (a = and condition (1.2)). Suppose that a = 0, T is of class H^, uq G 
H^{n), and condition (1.2) holds att = 0. Then, there exists T > 0, and a unique solution 
(v,q) of (1.3) with v e L°°{0,T;H^{Q)), q e L°°{0,T; H^{n)), and T{t) e H^. 

Because of the regularity of the solutions. Theorems 1.3 and 1.4 imply Theorems 1.1 and 
1.2, respectively. 

Remark 2. Note that in 3D, we require less regularity on the initial data than [13]. 

Remark 3. Since the vorticity satisfies the equation dtcurlu + i?„curlu = 0, where £u 
denotes the Lie derivative in the direction u, it follows that i/curlwo = 0, then curlw(i) = 0. 
Thus our result also covers the simplified case of irrotational flow. In particular, Theorem 
1.3 shows that the 3D irrotational water-wave problem with surface tension is well-posed. In 
the zero surface tension case, our result improves the regularity of the data required by Wu 
[21]. 

1.4. General methodology and outline of the paper. 

1.4.1. Artificial viscosity and the smoothed K-problem. Our methodology begins with the 

introduction of a smoothed or approximate problem (4.1), wherein two basic ideas are im- 
plemented: first, we smooth the transport velocity using a new tool which we call horizontal 
convolution by layers; second, we introduce an artificial viscosity term in the Laplace- Young 
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boundary condition (cr > 0) which simultaneously preserves the transport-type structure of 
the Euler equations, provides a PDE for which we can prove the existence of unique smooth 
solutions, and for which there exist a priori estimates which are independent of the artifi- 
cial viscosity parameter k. With the addition of the artificial viscosity term, the dispersive 
boundary condition is converted into a parabolic-type boimdary condition, and thus finding 
solutions of the smoothed problem becomes an easier matter. On the other hand, the a 
priori estimates for the k problem are more difficult than the formal estimates for the Euler 
equations. 

The horizontal convolution is defined in Section 2. The domain i7 is partitioned into co- 
ordinate charts, each the image of the unit cube in M^. A double convolution is performed 
in the horizontal direction only (this is equivalent to the tangential direction in coordinate 
patches near the boundary). While there is no smoothing in the vertical direction, our hor- 
izontal convolution commutes with the trace operator, and avoids the need to introduce an 
extension operator, the latter destroying the natural transport structure. The development 
of the horizontal convolution by layers is absolutely crucial in proving the regularity of the 
weak solutions that we discuss below. Furthermore, it is precisely this tool which enables us 
to prove Theorem 1.2 without the use of Nash-Moser iteration. To reiterate, this horizontal 
smoothing operator preserves the essential transport-type structure of the Euler equations. 

1.4.2. Weak solutions in a variational framework and a fixed-point, a > 0. The solution to 
the smoothed K-problem (4.1) is obtained via a topological fixed-point procedure, founded 
upon the analysis of the linear problem (4.2). To solve the linear problem, we introduce a few 
new ideas. First, we penalize the pressure function; in particular, with e > the penalization 
parameter, we introduce the penalized pressure function = iTr(a Vw). Second, we find a 
new class of [Hi (0)]'-weak solutions of the penalized and linearized smoothed K-problem in 
a variational formulation. The penalization allows us to perform difference quotient analysis 
in order to prove regularity of our weak solutions; without penalization, difference quotients 
of weak solutions do not satisfy the "divergence-free" constraint and as such cannot be used 
as test functions. Furthermore, the penalization of the pressure function avoids the need to 
analyze the highest-order time-derivative of the pressure, which would otherwise be highly 
problematic. In the setting of the penalized problem, we crucially rely on the horizontal 
convolution by layers to establish regularity of our weak penalized solution. Third, we 
introduce the Lagrange multiplier lemmas, which associate a pressure function to the weak 
solution of a variational problem for which the test-functions satisfy the incompressibility 
constraint. These lemmas allow us to pass to the limit as the penalization parameter tends 
to zero, and thus, together with the Tychonoff fixed-point theorem, establish solutions to 
the smoothed problem (4.1). At this stage, however, the time interval of existence and the 
bounds for the solution depend on the parameter k. 

1.4.3. Solutions of the K-prohlem, for (7 = via transport. For the cr = problem, we 
use horizontal convolution to smooth the transport velocity as well as the moving domain. 
Existence and uniqueness of this smoothed k problem (17.1) is found using simple transport- 
type arguments that rely on the pressure gaining regularity just as in the fixed-domain case. 
Once again, the time interval of existence and the bounds for the solution a priori depend 
on K. 

1.4.4. A priori estimates and n-asymptotics. We develop a priori estimates which show that 

the energy function £^^(0 in Definition 10.1 associated to our smoothed problem (4.1) is 
bounded by a constant depending only on the initial data and not on k. The estimates rely 
on the Hodge decomposition elliptic estimate (5.1). 
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In Section 10, wc obtain estimates for the divergence and curl of 77, v and there space and 
time derivatives. The main novelty lies in the curl estimate for 77. The remaining portion 
of the energy is obtained by studying boundary regularity via energy estimates. 

These nonlinear boundary estimates for the surface tension case a > are more com- 
plicated than the ones for the cr = case with the Taylor sign condition (1.2) since it is 
necessary to analyze the time-differentiated Euler equations, which is not essential in the 
a = case (unless optimal regularity is sought). 

We note that the use of the smoothing operator in Definition 2.1, where a double con- 
volution is employed, is necessary in order to find exact (or perfect) derivatives for the 
highest-order error terms. The idea is that one of the convolution operators is moved onto 
a function which is a priori not smoothed, and commutation-type lemmas are developed for 
this purpose. 

We obtain the a priori estimate 

sup E^{t) <Mo + TP{ sup E^{t)) , 
te[o,T] te[o,T] 

where Mq depends only on the data, and P is a polynomial. The addition of the artificial 
viscosity term allows us to prove that Ef^{t) is continuous; thus, following the develop- 
ment in [8], there exists a sufficiently small time T, which is independent of k, such that 
supig[o,T] E^{t) < Mo for Mq > Mq. 

We then find K-independent nonlinear estimates for the a = case for the energy function 
(20.1). 

Outline. Sections 2-15 are devoted to the case of positive surface tension ct > 0. Sections 
16-27 concern the problem with zero surface tension a = together with the Taylor sign 
condition (1.2) imposed. 

1.5. Notation. Throughout the paper, we shall use the Einstein convention with respect to 
repeated indices or exponents. We specify here our notation for certain vector and matrix 
operations. 

We write the Euclidean inner-product between two vectors x and y as x ■ y, so that 

X ■ y = x '' 

The transpose of a matrix A will be denoted by , i.e., {A'^Tj = A^. 
We write the product of a matrix A and a vector 6 as ^ 5, i.e, [A hf = A^^V . 
The product of two matrices A and S will be denoted by A-S, i.e., {A-Sfj = A\ . 
The trace of the product of two matrices A and S will be denoted hy A : S, i.e., 
{A : 5)5 = 4 S^. 

For O, a domain of class i?^ (s > 2), there exists a well-defined extension operator that 
we shall make use of later. 

Lemma 1.1. There exists E{Q,), a linear and continuous operator from 7J'"(f2) into i?'"(M^) 
(0<r< s), such that for any v € H''{n) (0<r<s), E{n){v) =v inO.. 

Wc will use the notation _ff*(ri) to denote either H^{Q.;M.) (for a pressure function, for 
instance) or H'^{il;M?) (for a velocity vector field) and we denote the standard norm of 
H'{0) (s > 0) by II • 11^. The H\Vi) inner-product will be denoted (•, •)«• 

We shall use the following notation for derivatives: di or denotes the partial time 
derivative, d denotes the tangential derivative on F (or in a small enough neighborhood of 
r), and V denotes the three-dimensional gradient. 
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Letting {x^,x^) denote a local coordinate system on F, for a = 1,2, we let either da or 
(•),„ denote We define: 

5" := gfdp , laVP = d"'d"- ■ ■ ■ d'^'da^da, ■■■3^, 

for integers A; > 0, where t/o = 5t=o is the (induced) metric on T. In particular, = \(f)\, 

\d^(j)\'^ = \d4>\'^ = d"(j)da(t>. d'^(j) will mean any fcth tangential derivative of (f). 

The area element on F in local coordinates is dSa = ^/godx^ A dx^ and the pull-back of 
the area element dS on F(i) = 77(F) is given by r]*{dS) — y^dSo- Let {J7i}|Li denote an 
open covering of F, and let {^i}fLi denote the partition of unity subordinate to this cover. 
The L'^{T) norm is 

\<t>\o := ||'/'IU^(r)= (^-/-'dSo 
and the if'^(F) norm for integers > 1 is 

m := ll'/'l|H'=(r)= (EEI^'^VIo 

\i=l 1=1 

Similarly, for the Hilbert space inner-products, we use 

» k K 

[(/),V']o := [</',V']L2(r) = / (jjipdSo, [4',^]k ■= [</», ^/']i/'=(r) = [4','<P]o+^^[^id'4',^id,tp]o ■ 

■^^ i=l 1=1 

Fractional-order spaces are defined via interpolation using the trace spaces of Lions (see, for 
example, [1]). 

The dual of a Banach space X is denoted by X', and the corresponding norm in X' will 
be denoted || • \\x'- For L G il*(0)' and v G H^{fl), the duality pairing between L and v is 
denoted by {L, v)s- 

Throughout the paper, we shall use C to denote a generic constant, which may possibly 
depend on the coefficient a, or on the initial geometry given by CI (such as a Sobolcv constant 
or an elliptic constant), and we use P(-) to denote a generic polynomial function of (•). For 
the sake of notational convenience, we will often write u{t) for u{t, •). 

2. Convolution by horizontal layers and the smoothed transport velocity 

Let O C K" denote an open subset of class H^, and let {f7,}^j^ denote an open covering 
of F := d^, such that for each i G {1, 2, K}, 

6i : (0, 1)^ X (-1, 1) — > Ui is an H'^ diffcomorphism , 

u^nn^e,{{o,if) and u^nr ^ e,{{o,if x {o}) , 

di{xi,X2,X3) = {xi,X2,ipi{xi,X2) + X3) and det V6'i = 1 in (0, 1)^ . 

Next, for L > K, let {CA}fiLK+i denote a family of open sets contained in fl such that 
{Ui}f=i is an open cover of CI. Let {ai}^i denote the partition of unity subordinate to this 
covering. 

Thus, each coordinate patch is locally represented by the unit cube (0, 1)^ and for the 
first K patches (near the boundary), the tangential (or horizontal) direction is represented 
by(0,l)2x{0}. 
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Definition 2.1 (Horizontal convolution) . LetO < p G 2?((0, 1)^) denote an even Friederich's 



mollifier, normalized so that / P = 1, with corresponding dilated function 

For w e H^{{0, 1)3) such that supp{w) C[S,1- S]"^ x (0, 1), set 

pi*hw{xH,X3)= PiixH -yH)w{yh,X3)dyH , yjf = (2/1,2/2) • 
We then have the following tangential integration by parts formula 



while 



Pi*hW,aixH,X3) = P^,a{xH - yH)w{yh,X3)dyH , a = 1,2 , 
Pi*hW,3{XH,X3) = / pi{XH - yH)w,3{yh,X3)dyH ■ 

It should be clear that -kh smooths w in the horizontal directions, but not in the vertical 
direction. Fubini's theorem ensures that 

*h w||s,(o,i)3 < C's||w||s,(o,i)3 for any s > , (2.1) 

and we shall often make implicit use of this inequality. 

Remark 4. The horizontal convolution -khW does not smooth w in the vertical direction, 
however, it does commute with the trace operator, so that 



(pi -khw) 



(0,1)2 x{0} 



— pi *h W'l(0,l)2x{0} 



which is essential for our methodology. Also, note that -kh smooths without the introduction 
of an extension operator, required by standard convolution operators on bounded domains; 
the extension to the full space would indeed be problematic for the transport structure of the 
divergence and curl of solutions to the Euler-type PDEs that we introduce. 

Definition 2.2 (Smoothing the velocity field). For v G L'^{fl) and any k e (0, ^) with 
Ko = min dist (supp{ai o Oi) , [(0, 1)^ x {0}]'= n d[0, 1]^) , 

i=l 

set 

K L 

= 51 Pi ->'h. [Pi *h {{Voiv) o 9i)] o 0-^ + ^ aiV. 

i=l i=K+l 

It follows from (2.1) that there exists a constant C > which is independent of k such 
that for any v E H'^{fl) for s > 0, 

\K\\s<C\\v\\s and K\s-i/2<C\vU_i/2. (2.2) 

The smoothed particle displacement field is given by 

77« = Id + / v^. (2.3) 
Jo 
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For each x £ Ui, let x ~ 6^ ^{x). The difference of tlie velocity field and its smoothed 
counterpart along the boundary F then takes the form 

v^{x)-v{x)=y^ [ [ Ci(a^)Pi(y)Pi(5)[(Ci^^)(^.(i-(y + 5)))-(Ci^;)(e.(i))]d5dy, 

(2.4) 

where C,i{x) = yjai{9i{x)). Combining (1.1a), (2.3), and (2.4), 

■n^{x) - r,{x) =y^ I I Q{x)pi{y)pi{z)[{CiVM{x-{y + z)))-{Cir]){ei{x)^^ 

(2.5) 

For any u e H^-^{r), and for y e B{x,k), where B{x,k) denotes the disk of radius k 
centered at x, the mean value theorem shows that 

\u{y) - u{x)\ < C\r~'^\L.i(^B(x,K))\du\Lv(B{x,K)), r = radial coordinate, 
so that in particular, with p = 4 and 5 = |, 

\u{y) - u{x)\ < C\fK\du\j,i < C«;|u|i.5 , 
the last inequality following from the Sobolev embedding theorem. Hence, for U e il^'^(F), 

\U^{x) - U{x)\loo < Cx/^|C/|i.5 . (2.6) 
Note that the constant C depends on maxjgj]^ ^^-j |^i|5.5- 

Letting = ^^'^ ^ ~ (^' ■'■)^' ^® have that for any </> G -^^^(r), 
. K . ^ r 

/ Vf,(f> = Y2 / Px*h pi- *hCiv{x)Ci(l){x) =Y2 / P^-^hQiv{x) pL-khC,i<i>{x) 

f ^ 

= / Yl^P^ *h iCiV o Oi)] o {Qcp o 9i)] o . (2.7) 

Finally, we need the following 

Lemma 2.1 (Commutation-type lemma). Suppose that g G i^(L) satisfies dist{supp{g),dR) < 
Ko) and that f € iJ*(F) for s > 1. Then independently of k G (0, kq), there exists a constant 
C > such that 



pi *h [fg] - fpi *hg <C K,\f\s+i^R \g\o,R . 

0,R 



We also have 



PL */i [fg] - fpi *hg ^ < Ck||/L+|,[o,i]3 llffllo,[o,i]3 
< mm{dist{supp fg, {1} x [0, 1]^), dist{supp fg, {0} x [0, 1]^)). 



whenever g € L^i^), f € 7f^(f2) and 

K 

2 

Proof Let A = p i ★fc [fg] - fpi -kh g- Then 



|A(x)| 
so that 



P±ix - y)[f{y) - f{x)]g{y)dy 
b{x,k) 



|A|o,ij< CK|/|,+i,fl 



Pi 151 



<C'K\f\s+i,R pi{x-y)\g{y)\dy, 
Jb{x,k) " 

< CK[f\s+i,R [g[o,R. 
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The inequality on [0, 1]^ follows the identical argument with an additional integration 
over the vertical coordinate. The hypothesis on the support of fg makes the integral well- 
defined. □ 



Remark 5. Higher-order commutation-type lemmas will he developed for the case of zero 
surface tension in Section 21. 



3. Closed convex set used for the fixed-point for > 

In order to construct solutions for our approximate model (4.1), we use a topological fixed- 
point argument which necessitates the use of high-regularity Sobolev spaces. In particular, 
we shall assume that the initial velocity uq is in H^^-^{^) and that ^ is of class C°°; after 
establishing our result for the smoothed initial domain and velocity, we will show that both 
Q. and uq can bo taken with the optimal regularity stated in Theorem 1.3. 

For T > 0, we define the following closed convex set of the Hilbert space 1/^(0, T; H^^-^{p,)): 

Ct = {v& L^(0, T; H^^-^{Q))\ sup ||i;||i3.5 < 2||wo||i3.5 + 1}, 

[0,T] 

It is clear that Ct is non-empty, since it contains the constant (in time) function 7io, and is 
a convex, bounded and closed subset of the separable Hilbert space L^(0, T; iJ^'^-^(r2)). 

Let V e Ct be given, and define t] by (1.3a), the Bochner integral being taken in the 
separable Hilbert space _ff ^•^ "''(fi). 

Henceforth, we assume that T > is given such that independently of the choice of 
V e Ct, we have the injectivity of r]{t) on O, the existence of a normal vector to ry(0, t) at 
any point of rjiV , t), and the invcrtibility of V77(i) for any point of fl and for any t € [0, T]. 
Such a condition can be achieved by selecting T small enough so that 

l|V?? - Id||i,»(o,T;Hi3.5(Q)) < eo , (3.1) 

for eo > taken sufficiently small. Condition (3.1) holds if T||Vuo||h2 < cq. Thus, 

a=[Vr/]-i (3.2) 

is well-defined. 

Then choosing T > even smaller, if necessary, there exists kq > such that for any 
K e (0, ^), we have the injectivity of riK,{t) on fl for any t G [0, T]; furthermore, Vrj^ satisfies 
the condition (3.1) with rj^ replacing r]. We let n^iilKix)) denote the exterior unit normal 
to r]K{fl) at rjtiix) with x gT. 

Our notational convention will be as follows: if we choose v G Ct, then fj is the flow map 
coming from (1.3a), and a is the associated pull-back, a = [Vfj]~^. Thus, a bar over the 
velocity field, will imply a bar over the Lagrangian variable and the associated pull-back. 

For a given v^, our notation is as follows: 

r]i^{t) =ld+ / and r/K(0)=Id, 
Jo 

a« = Cof Vr?K , J« = det Vr?^ , QKafS = ^cVk ■ d^rj^ . 
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We take T (which a priori depends on n) even smaUer if necessary to ensure that for 

ie [o,T], 

VgW^ < '^Va'o"^ , (3.3a) 



< 2v^" . (3.3b) 

I < Ut) < I ■ (3.3c) 

Lemma 3.1. For v e Ct, and for any s>0, we have independently of the choice ofvG Ct 
that 

sup \v^\s < C«,sf(||wo||i3.5). 

[0,T] 

Proof. By the standard properties of the convolution a.e in [0,T]: 

K\s < [-^ + l]b|i3 < + l][2||«o||i3.5 + 1], (3.4) 

where we have used the definition of Ct for the second inequahty. □ 

Recall that {Oi}f^^ is our open cover of T. Given v E Ct, wc define the matrix = 
[^{Vk ° f^i)]^^ , and assume that T > is sufficiently small so that independently oi v G Ct, 
we have the following determinant- type condition for 6^: 

l<(ti)ltM)ir, in (0,1)^. (3.5) 

i=l 

Such a condition is indeed possible since at time t = we have (^L)! Si=i[(^K;)f]^ = 1 + 
V'z,i+V';,2- 



4. The smoothed k-problem and its linear fixed-point formulation 

Unlike the case of zero surface tension, for a > there docs not appear to be a simple 
sequence of approximate problems for the Euler equations (1.1) which can be solved only 
with simple transport-type arguments. For the surface tension case, the problem is crucially 
variational in nature, and the addition of an artificial viscosity term on the boundary T 
seems unavoidable in order to be able to construct a sequence of approximate or smoothed 
solutions. 

As we shall make precise below, our construction of the approximating sequence of prob- 
lems is based on smoothing the transport velocity by use of the horizontal convolution by 
layers (see Definition 2.2), and hence smoothing the Lagrangian flow map and associated 
pull-back. Simultaneously, we introduce a new type of parabolic-type artificial viscosity 
boundary operator on F (of the same order in space as the surface tension operator). Note 
that unlike the case of interface motion in the fluid-structure interaction problem that we 
studied in [8], there is not a unique choice of the artificial viscosity term; in particular, other 
choices of artificial viscosity are possible for the asymptotic limit as the artificial viscosity 
is taken to zero. 
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We can now define our sequence of smoothed K-problems. For our artificial viscosity 
parameter /t e (0, let {v,q) be the solution of 

77 = Id + /" V in r2 X (0, T] , (4.1a) 
Jo 

dtv + J^^a^Vq = Q inOx(0,T], (4.1b) 

Tr(a« Vv)=Q in Q x (0, T] , (4.1c) 

-o^ -^^^giv) ■ nkivii) rikivK) - kAo[v ■ nK(77K)]nK(??K) = gnK(r?„) on T x (0,T] , (4.1d) 
v 9 k 

{t], v) = (Id, uo) on f2 X {t = 0} , 

(4.1e) 

where n^{r]^) = and Aq = ^/gZ~^^a[^/go9o^^/3]. Note that on T, ^ = \alN\, 

and that (g„)a/3 = rj^,^ ■ri^,0. 

In order to obtain solutions to the sequence of approximate K-problems (4.1), we study 
a linear problem whose fixed-point will provide the desired solutions. If we denote by v 
an arbitrary element of Ct, and 77^, Ok, and are the associated smoothed Lagrangian 
variables given by Definition 2.2, then we define w to be the solution of 

dtw + J-^a^ = innx (0, T] , (4.2a) 

Tr(a«, Vw) =0 in x (0, T] , (4.2b) 

-CT-^[Ag(r7) • n«(^«)]n«(r7«) - kAo[m; • n«(77K)]n«;(77«;) = gn^ on T x (0,T], (4.2c) 
V 9k, 

(77, w) = (Id, Uo) on f2 X = 0} , 

(4.2d) 

where 5^/3 = j?,a -J?,/?, and Aq = Vg^^^da[^9o^df3]. 

For a solution w to (4.2), a fixed point of the map v 1-^ w provides a solution of our 
smoothed problem (4.1). 

In the following sections, we assume that v e Ct is given, and k is in (0, ^). Until 
Section 10, wherein we study the asymptotic behavior of the problem (4.1) as k — > 0, the 
parameter k is fixed. 

5. Hodge decomposition elliptic estimates 

Our estimates are based on the following standard elliptic estimate: 

Proposition 5.1. For an H"" domain Q, r > 3, if v e L'^{fl) with curlt; e H''~^{Q), 
divv € if''~^(ri), and ii • iV|r € H^~^{T) for 1 < s < r, then there exists a constant C > 
depending only on ft such that 

\\v\\s < C [\\v\\q + II curlu||s_i + II divu||s_i + \v ■ N\^_i J , 

> (5.1) 

ll^'lls < C M|f ||o + II curlulls-i + II divt;||s_i + \v ■ Ta\s-^ \ , 

where Ta, a = 1,2 are the tangent vectors to F. 

The first estimate with V ■ N is standard (see, for example, [19]), while the second with 
V ■ Try follows from the fact that T„-N = 0. 
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6. Weak solutions for the penalized problem and their regularity 

The aim of this section is to estabhsh the existence of the solution Wc to the penaUzed 
version (of the divergence- free condition) of hncarizcd and smoothed K-problem (4.2). In 
particular, we study the weak form of this problem with the pressure function q, approxi- 
mated by the penalized pressure 

q" = -^Trid^Vw) for < e « 1 . 

In this section, as well as in Sections 8 and 9, we let 

N{uo,x,y) = P{\\uo\\i3.5,x,y) (6.1) 

denote a generic polynomial function of ||uo||i3.5, x, and y, where x and y will typically 
denote norms of various quantities. 

6.1. Step 1. Galerkin sequence. By introducing a basis (e;)^^ of H^{fl) and L^{fl), 

I 

and taking the approximation at rank I > 2 under the form wi{t,x) = ^^2/fe(f) ek{x) , 

fe=i 

satisfying on [0, T], the system of ordinary differential equations 

(i) ( Wit, 'P)o + k[wi ■ n^{f]^),cj) ■ n^{f]^)\i - a[Lgfj ■ n^{fjf,), (j) ■ nK(r^„)]o 

- {{^n)Uh 4>\i )o = 0, V(/) e span(ei, ei) , 

(ii) wi{0) = {uo)i, in O , 
1 

,\- r,. - - 

span{ei, ...,e;), we see that the Cauchy-Lipschitz theorem gives us the local well-posedness 
for Wl on some [0,Tmaa:]- The use of the test function wi in this system of ODEs (which is 
allowed as it belongs to span{ei, e;)) gives us in turn the energy law for any t G (0, Tmax), 

i||J|u;i(i)||o + K / [wi ■ nK{'nK),wi ■ fiKiVK)]! + e I \\qi\\l 
^ Jo Jo 



where Lg = -^Ag, qi = — {aK,)iwl,j, and {uo)i denotes the i^(ri) projection of uq on 



~2j (("^«)*"'''"'')o " 2"*'"°'''"° ^ [LgVni^{'nK),wi-n^{r]^)]o, 
which, with the control of ry" provided by the definition of Ct, gives the bound 

\\\wi{t)\\l + CKj^ \wi-n^{fj^)\l + ej^ \\qi\\l<CN{uo). (6.2) 

6.2. Step 2. Weak solution We of the penalized problem. We then infer from (6.2) 
that Wl is defined on [0, T], and that there is a subsequence (still denoted with the subscript 
/) satisfying 

wi^w^ in L^{0,T:L\Q)), (6.3a) 
qi^q, mL\0,T;L\n)) , (6.3b) 



where 



We can also rewrite (6.3) as 



ge = -^(a«)>e4- • (6-4) 



Wl We in L^(0, T; L^(f2)), (6.5a) 
div(w;, o fj-'){f},) - div(w;, o fj-'){f}^) in L\0, T; L2(0)) , (6.5b) 
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which with the bound (6.2) and the definition of the normal provides 

wi ■ fi^f]^) We ■ n^fj^) mL^{0,T;H\r)). (6.6) 

It follows from standard arguments and the ODE defining wi, that w^t € ^^(0, T; H^{iY)'), 
We G C°{[0,T]■,H^ny) with We{0) = uo, and that for G L^{0,T; {Q}), 

{Ji^Wet,(j>)i+ K I [9(we •«■«(??«)), 9(<;i-n«(r?«))]o 
Jo 

T i-T 

{qe,{a^yi(l)\j)i =a [Lgfj ■ n^fj^) , cj) ■ n47]^)]Q . (6.7) 
Jo 

Since by definition = CofV?7K, this implies that in fl, 

Wet+VpeiV^) ^0, (6.8) 

where Pe o f]K, = Qe in ^- Since VpeifjK) G L^(0, T; J?~^(f2)), this equality is true in 
L2(0,T;i?-i(f7)) as well. 

6.3. Step 3. is bounded in L^(0, T; iJ'^(O)) independently of e. Denoting = 
We o fj~^, by integrating (6.8) in time from to t, we obtain the important formula 

cmlueifj^) = cmluo +[ B{u^,Ue){fj^) mL'^{0,T;H-\n)), (6.9) 
Jo 

with 

B{UK,Ue) = -«,2 We3,i — U^,3Ue2« , <,3 Wei « -<,1 "eSn , <,1 We2n -<,2 W^i.j ). 

Remark 6. Note well that our approximated and penalized n-prohlem preserves the structure 
of the original Euler equations as can he seen hy (6.8). As a result, (6.9) contains only first- 
order derivatives of the velocity. 

Our next task is to prove that We in i^(0, T; H^(p.)). For suppose that this was the case; 
then, (6.9) together with bounds on the divergence of We and We ■ N on T, provide bounds 
for We in L'^{0,T; H^{Q)) (by the Hodge elliptic estimate (5.1)) which are independent of 
e > 0. 

We proceed by showing that appropriately convolved velocity fields are bounded inde- 
pendently of the parameter of convolution in L^(0, T; H^{ft)). This is the first instance that 
our horizontal convolution by layers is crucially required. 

6.3.1. For any subdomain w CC fi, We G L'^{0, T; H^{lo)). Wc analyze the third component 
of (6.9), the other components being treated similarly. This leads us to the following equality 

in L^iO,T; H-^{n)): 

(a«)^u;e,] -{ai^){we,'j = - curlug + / [ -w,} (a«)^u;e,/ (««)' + v,) (a„)ju;e,f (a«)^]. 

Jo 

Our goal is to prove that We G H^{Q). To proceed, we let ap denote a standard sequence 
of Priederich's mollifier in with support B{0, 1/p), and establish that cjp-kWe is bounded 
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in H^{uj) for any uj CC il. For this purpose, we choose ijj E and find that 



2 

ft 



Jo 

(6.10) 

In order to proceed, we shall need to identify curl-type structures (in Lagrangian vari- 
ables) for (jp-kw^; this requires the following: for - < dist(supp'^, O''), and / e C°°{Q.), we 
have the equality in H~^{Q.) 

Op -k [f{ij;we),j ]- fcTp-k [{i^We),j ] = / (crp),j {x - y){f{y) ~ f{x))ipwe{y)dy 

JR3 



/ 



<^p{x - y)f,j {y)ipwe{y)dy , 

showing that ap -k [f{tpWe),j ] — / * [{ipWe),j ] S L'^{^), with 

\\ap-k[f{il)We),j] - f crp*[(V'We),j]||o < C[ \\a,j ||o,M3 + IkHo.iRs] || V/||i,c»(a) ||u;e ||o. (6.11) 
We thus infer from (6.10) and (6.11) that the vorticity structure satisfies 

(a^)^o-p ★ (V'wJ,] -{d^){(Jp * {ipw^),'l = / [ -u,} {ai^)lap ★ {il^w^),] (a«,)' 

Jo 

+ V,) {a^){op ★ (?A«;,),2 {a^)\] + (6.12) 

with \\Ri\\l'^(o^t;L'^{q.)) < N{uo), where N{uo) is defined in (6.1). Next, we infer from (6.5) 
and (6.11) that the divergence structure satisfies 

{d^)ic7pk{iPw,),i = R2, (6.13) 

with ||ii2||L2(o, T:L2(o)) < N{uo). Since we also have ipWe = on F, so that with (5.1), we 
have a.e. in (0, T) 

\\apk{lljWe){t)\\i<\\Ri{t)\\o + \\R2{t)\\o+N{uo) ( ||(Tp*(^W;e)||l, 

Jo 

and thus 

|2 



/ 

Jo 



\\(Tp^{il,w,)\\i<N{uo). (6.14) 



Since this inequality does not depend on p, this implies that ipw^ € L'^{0,T; H'^{Q.)), and 
therefore G L'^{0,T; H^{ll))), with an estimate depending a priori on u) CC fl. 

6.3.2. r/ie horizontal convolved-by-layers velocity fields are in L^{0,T; H^{^l)). 
Fix I e {1, K}, and set 

VF(Z) = w,oei and = [V(^« o 61)]-^ 

Hence, in (0, 1)^ x (i, 1) for p > 1, the Lagrangian "divergence-free" constraint is given by 

{bi)i{aiW{l)),i = -(V^{ai,j Wiiy - aieqeiOt) , (6.15) 
where the crucial observation is that the right-hand side of (6.15) is in L^{0, T; i^([0, 1]^))- 
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Now for — < dist(suppQ;;, 9(0, 1) x (0, 1)^), and / smooth in [0, 1]^, we have by Lemma 
m 

2.1 that for /3 - 1,2, p„ *^ [f{aiW{l)),^\ - fp^*h [{aiW{l)),0] G ^^((0,1)2), with the 

estimate a.e in (-,1): 
P 

< C'/9|V/|l,oo((0,l)2x{j;}) |W'(/)|o,(0,l)2x{j/}- 

This leads to 

\\pm*h [f{aiW{l)),i3 ]- f Pm*h [{aiW{l)),,3 ] |lo,(0,l)3 

<Cp||V/|U»((o,i)3)||W-(0llo,(o,i)- (6-16) 
Now, for the case of the vertical derivative, we will need to express W{1),3 in terms of W^(/),i, 
W{1),2, cur^^oe, W{1) and divfj^oe, W{1), where 

div^.oe, W{1) = {bi)iW{l),i , 
cmll^o, W{1) = {bl)iW{l),l -{bl)iW{l),', , 
cur4„,, W{1) = {bi)iW{l),] -{bi)\W{l),^ , 
cur4„,, W{1) = {bl)\W{l),', -{bDiWil),] . 
Notice that the first three lines above can be written as the following vector field: 

3 

(div^^o^, W{1), cm\l,,^ W{1), cmll^e, W{1)) = ^ M,^W{IU , 

where the are smooth matrix fields depending on bj^. Prom condition (3.5), since 

detMS = {Fjl ^[(&'J|P> J, 
i=i ^ 

we see that is invertible on [0, T] (regardless of the choice of v G Ct)- Therefore, 

2 

W{1),3 = div^^oe, W{1) V- + M« curl^^oe, W{1) + ^ A'yW{l),^ , (6.17) 

i=l 

where M'^ and the Af arc smooth matrix fields depending on b\^, and is a vector field 
depending on 6j,. From (6.10), we have that 

curl^-^o0, W{1) = cmluoiei) + V / TVf W(Z),i . 

where the are smooth matrix fields depending on 6^. By using (6.17) and the fact that 
divj^^ogi e Z/^(0,T; L'^{fl)) from (6.5), we obtain after time differentiating that 

2 

[curl^^o,, W{l)]t - N^M^ curl^„o^, W{1) = ^ P^WiJU , +NSV^ div^.o^, W{1) 

0=1 

where P^, = 1,2, are smooth matrix fields depending on b^.. Therefore, 

curl^^oe, W{1) = A'' cm\uo{ei) + A'' [ {B^W , + diYfj^oe^ W{1)) (6.18) 

^0 



18 



D. COUTAND AND S. SHKOLLER 



where A'^ and Bj^, (3 = 1,2, are smooth matrix fields depending on With (6.16) and 
(6.18), we infer in a similar way as for (6.13) that on (0, 1)^ x (-, 1) we have 



2 /.* 



curl^^oe, Pm *h [aiW{l)] = A'^S^ j B^p^ *h [ai o 9iW{l)y +R3, (6.19) 

/3=i-^o 

with ||-R3||l2(o,t;L2((o,i)3)) < ^(wo). Therefore, with (6.17) and (6.18), we have that 

2 „t 2 

[ai{9i)W{l)U = M-A''Yl / B-0[ai{ei)W{l)],i} + J2^0[MOi)W{l)],i}+R4, (6.20) 

/3=l-^0 (3=1 

with ||-R4||l2(q .p.^2((o,i)3)) < N{uo). Thus, for any test function G -f^^((0, 1)"^), 

/ ai{ei)W{l) ■ip= f [ai{ei)W{l)U I ai{ei)W{l) ■ <p,3 . 

AoAVxj; J(0,l)2x(i,0) J(0,l)2x(i,0) 

Now, since for /3 = 1,2, we have 



[ai{ei)W{l)y = - / [ai{9i)W{l)] ■ ip,p , 

(0,l)^x(i,0) ./(0,l)2x(i,0) 

using (6.20), we infer that 

ai{9i)W{l) • <p| < C (||W(/)||o,(o,i)3 + ||J?4||o,(o,i)3)|k||i,(o,i)3, 

(0,l)=xi 

implying (independently of p > 1) the following trace estimate for W{1) (not just its normal 
component): 

\ai{ei) T4-(OI'_i,(o,i)2xi < ^M- (6.21) 
Similarly as (6.19), we also have the divergence relation 

2 .t 

dwfi^oe, Pm *h [aiW{l)] = C"^ / ^0P"^ *h h ° 0iW{l)],i3 +R5, (6.22) 

with ||i?5|lL2(o,T:L2((o,i)3)) < N (uq) , and and _D^, 13 — 1,2, are smooth matrix fields in 
terms of From (6.19) and (6.22), we then infer, just as in (6.14), that 

/ \\pm ★fc [aiW{l)] o {fj^ o 9i)-^\\l < N{uo) + [ \pm *h [aiW{l)] o (f?« o Oi)-^ ■ n«|i 

Jo ' " Jo 2' p 

where = ^K(0, 1)' x (i, 1)). Thus, 

£ WPm^h [";(^;)^(0]IImo,i)2x(|,i) < ^("o) + 1'^-*'' iMOi)wm\,^o,ir^^- 

Now, from the properties of the convolution, 

^ \pm*h [ai{9l)W{l)] I 1,(0,1)2x1 < C\pm*h h(^0^(0]l-i,(O,l)=xi ' 

which, with (6.21), leads us (independently of p > 1) to 

^ £ \\Pm *h [«/(^OW^(0]ll?,(0,l)2x(i,l) < ^(«0), 
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for any < — < dist(suppa;(0;), 9(0, l)'^ x (0, 1)). Since this estimate holds for any p > 1, 
m 

we then infer that 

Jo 

for any < — < dist(suppa;(^;): ^(0, 1)^ x (0,1)). Therefore, pm [ai{9i)W{l)] e 

if^((0, 1)^) (which was not a priori known since our convolution smooths only in the hori- 
zontal directions), with a bound depending a priori on m. 

6.3.3. Control of the horizontal convolved-by-layers velocity fields independently of m. 
Prom (6.19) and (6.22), we infer that 

Jo 

< N{uo)+ [ \Pm*h[ai{ei)W{l)]o{f],oei)-' -n^ll^Q-^^^^, 



and thus. 



[ \\Pm*h[ai{ei)W{l)]\\l^o,i)- (6.24) 

< N{uo)+j^ |Pm*/i[a;(6li)M^(0]-"«feo^i)l|,(o,i)2x{o}- 
Next, we have for any x € (0, 1)^ x {0}: 

Pm *h [ai{ei)W{l)] ■ n^{fji^ o ei){x) = Pm *h [ai{ei)W{l) ■ n^{fj^ o ei)]{x) + f{x), 

with 

f{x) = / pm{xH - yH)ai{ei)W{l){yH,X3) ■ [««(??« o ei){xH,X3) - n^ifj^ o 9i){yH,X3)]dyH- 
Jm? 

Therefore, with (6.24), we obtain 

rWPm^h [ai{ei)wm\l^o,ir 

Jo 

< N{uo) + I/I i,(o,i)2x{o} + / \Pm *h [ai{ei)W{l) ■ n^fj^ o 6'i)]||,(o,i)2x{o} 

J 

< N{uo) + I/I 

i,(0,l)2x{0} ,(0,l)2x{0} 

< N{uo) + I/I i,(o,l)2x{0} + ^ lo^i^e ■ n«(f7«)||_p 

<^(«o) + |/|i,(o,i)2x{o}, (6-25) 

where we have used the trace control (6.6) for the last inequality in (6.25). We now turn 
our attention to |/| i.(o,i)2x{o}- We first have that 

||/||o,(o,i)3 < ^\\n.in.)\\HHn)\\prn^hai{ei)\Wm\o,io,ir ^ ^^(^o)' (^-^S) 
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where we have used the definition of Ct to bound ||nK(j?K)||//3(Q). Next, we have for /3 = 1, 2 
that 

(x) = / P-m,f3 {xh - yH)oii{0i)W{l){yH,X3) ■ [n„(f?„ o ei){xH,xs) - n^ifj^ o ei){yH,X3)]dyH 
+ / Pm{xH-yH)ai{9i)W{l){yH,X3)dy-n^{fji^oei),0{x), 
showing that 

C 

\\f,0 llo,(o,i)3 < — ll?i«fe)llff3(o) ^ II \pm,f3 1 *h ai{ei)\w {ly \\\o,(o,ir 



m . 

1=1 



+ ll"«fe)||i/3(n) J2 \\pm *h ai{0i)\W{iy\\\o,(o,i)3 

i=l 

3 

< C\\n^iv.)\\HHn) ^ \\\{p,p U *h ai{0i)\W {If \\\oxo,i)^ 

i=l 

3 

+ ll^«^fe)l|ff3(o) ^||pm */i a/(6'i)|W^(01llo,(o,i)3 

3 3 

<cY,\\\iP,0U*h a;(^0im0illo,(o,i)3+c5^||p™^fta/(^0im0illo,(o,i)3 



< C||a,(0OW^(Ollo,(o,i)^ < ^Wo). (6.27) 
Next, for the vertical derivative, 

/,3 (a;) = p„(ar^^ - yH)[aK^0W^(0].3 iyH,xs) ■ o ^O]^;^'^?'^^ 

+ j^^ Pm{xH - yH)[aii9i)W{l)]{yH,X3) ■ [n«(?7« o 6'i),3 ll^^^^^jdy^ , 

where [•ll^"'!^!,^) = [jixHjXa) — [-liyHjXs). Notice that for a smooth matrix field A in (0, 1)^ 
and for $ = 1,2, 

G{x) = [ pm{xH - yH)A{yh,X3)[ai{ei)W{l)],0 {yH,xs) ■ [n,{fj. o 0i)t"'Z)dyH, 
satisfies 



G{x) = J {p,!3)m{xH - yH)A{yh,x:i)[ai{ei)W{l)]{yH,X'i) ■ [^^(t]^ o eify^^y'^'JJ^dyH 

- f Pm{xH - yH)A,D {yh,xz)[ai{ei)W{l)\{yH,xs) ■ [n«fe o ei)][l^f ]dyH 

- / pmixH - yH)A{yh,X3)[ai{ei)W{l)]{yH,X3) ■ [n^ifj^ o ei),f3&"'2)'iyH, 

showing, just as for (6.27), that ||G||o,(o,i)3 < N{uo)- Therefore, with (6.20), we see that the 
first integral term appearing in the expression of /,3 is bounded in a similar way, implying 
that 

ll/,3||o,(o,i)3 < A^(wo). (6.28) 



WELL-POSEDNESS OF THE FREE-SURFACE EULER EQUATIONS 



21 



Consequently, with (6.26), (6.27), (6.28), we obtain that 

||/||i,(o,i)3 <7V(uo). (6.29) 

Therefore, (6.25) imphes that 

/ \\Pm*h[ai{ei)W{l)]\\l^o^,)s<N{uo). (6.30) 
6.3.4. Control of We in L'^{0,T;H^{Q)). 

Since (6.30) holds independently of m sufficiently large, this implies that 

£ \\ai{ei)W{l)\\l^,^,^s < N{uo). 

Since we proved in subsection 6.3.1 that is bounded in L'^{0,T;H^{w)) independently of 
e for each domain u> CC fl, this provides us with the estimate 

/ \\we\\i<N{uo), (6.31) 
Jo 

independently of e > 0. 

Remark 7. In the two-dimensional case, a simpler proof of Step 3 is possible, founded upon 
a scalar potential function for the velocity field. For conciseness, we consider a simply- 
connected domain, the non-simply connected case being treated similarly by local charts. 
Once again, we let = ofj^^- From (6.5b) and (6.6), let w% G L^(0, T; i/^(0)) such 
that 

div«(r7-i ))(,?«) = dW{w,{f)-'m,) m L\0,T; L^n)), 
w% ■ n,^{f]^) ^ ■ ni.,{f]i.,) in L^(0, T; H'^{T)) . 

We infer the existence oftp"^ G ^^(0, T; iJg (?^„(n))) such that = w^{fi~^) + {—il}^,2 , "4^^,! )■ 
Now, from (6.10), we see that in L'^{0,T; H~^{fl)), we have for 'ip^ = o fj^, 

-{a^)Hia.)lr,j ),k =fr- f AWni , (6.32) 

Jo 

where f^ is bounded in L'^{O.T: L'^( fl)). It is readily seen that ijf is the unique solution 
of this equation in L^{Q,T] Hq{Q.)) . We now establish that this uniqueness provides extra 
regularity for 'ip^. By defining the mapping © from L^{0,T; H^{Q) n Hq{CI)) into itself by 
associating to any ^ in this space, the solution Q(, (for almost all t G [0, T]) of 

-{a,)H{a.){e^,j ),k =fr- f A^iUj , 

Jo 

we see that for ti small enough (depending on Sobolev constants and on \\uK\\L'^(a.T-H^{n))) 
9 is contractive from L^(0,ti; II^{Q) n IlQ{n)) into itself, which provides a fixed-point for 
© in this space. It is thus a solution of (6.32) on [0, ii]. By uniqueness of such a solution, 
we have that tp^ G L^{0,ti; H^{^1)) and thus that E L^{Q,ti; H^{n)) . By defining a 

mapping similar to ©, but this time starting from t2 G ["^'^i] such that Wc{t2) G H^{Q) 

instead ofuo (which ensures that the new fr is still in L^{0,t2; L'^{0,))), we obtain the same 

conclusion on [t2,t2 + ti], leading us to G L^{0, -ti; (fl)) . By induction, we then find 

WeGL^{0,T;H^Q)). 
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Remark 8. Whereas Hodge decompositions with vector potentials ip are possible in higher 
dimension, it turns out that a, Dirichlet condition 1^ = for the associaied elliptic problem 
is not possible. This in turn is problematic for any uniqueness argument in L^(0,T; H^{Q)) 
for tp, since it does not seem possible to find a boundary condition that would be naturally 
associated to the second order operators appearing on both sides of the three-dimensional 
analogous of (6.32). 

7. Pressure as a Lagrange multiplier 

We will need two Lagrange multiplier lemmas for our pressure function in our analysis 
as the penalization parameter e ^ 0. We begin with a lemma that is necessary for a new 
Hodgc-typc decomposition of the velocity field. 

Lemma 7.1. For all I G H^{^), t G [0,r], there exists a constant C > and (p{l) G Hi{n) 
such that {aK,)l{t)(j}^ ,j = I in Cl and 

\\m\i<c\\i\\i. (7.1) 

Proof. Let ^p{l) be the solution of 

{a,)l[{a,)':m,k],j=linn (7.2a) 
i>{l) = on r. (7.2b) 

We then see that </>'(?) = {aK)'l'fp{l),j satisfies the statement of the lemma. The inequality 
(7.1) is a simple consequence of the properties of / and of the condition v G Ct- □ 

We can now follow [18]. For p G H^^Q)', define the linear functional on H^{Q) by 
(P) (^K;)i (0'^*'j )|' where f G il2(0). By the Riesz representation theorem, there is a 
bounded linear operator Q{t) : {H^{Cl)y Hi{fl) such that 

yip€Hi{n), {p, {a,){{t)ip\j)i={Q{t)p, ^)s. 

Letting ip = Q{t)p shows that 

IIQ(*)p|||<C|b||^i(^^, (7.3) 
for some constant C > 0. Using Lemma 7.1, we see that 

yiGH^n), {p, i).={Q{t)p, (/.(0)3, 

and thus 

lbl|^i(,,), <C||Q(t)p|||, (7.4) 

which shows that R{Q{t)) is closed in iJ 1(0). Let Ve(<) = {v G L'^{n) \ {a^)l{t)v\j {t) = 0}. 
Since Vy{t) n Hi (f2) = R{Q{t))-^, it follows that 

Hi{n) = R{Q{t)) e^i^^^ v,{t) n h'^q). (7.5) 

We can now introduce our first Lagrange multiplier 

Lemma 7.2. Let £{t) G ili(O)' be such that £{t)(p = for any (p G Vs(t) n Hi{n). Then 
there exists a unique q{t) G H^{Q)' , which is termed the pressure function, satisfying 

V<^ G ifi(n), Z{t){,p) = {q{t), {a,)W\j)i^. 
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Moreover, there is aC > (which does not depend ont ^ [0, T] and on the choice ofv^ C't ) 
such that 

Proof. By the decomposition (7.5), for tp e M^^^ ^et (p = vi + V2, where vi € 

Vv{t) n iJ2 (O) and V2 € R{Q{t)). From our assumption, it follows that 

for a unique tplt) G R{Q{t)). 

From the definition of Q{t) we then get the existence of a unique q{t) e H2{ny such 
that 

e Z{t){^) = {qit), id,)iv\j).. 

The estimate stated in the lemma is then a simple consequence of (7.4). □ 

We also need the case where the pressure function is in Hi{fl). We start, as above, with 
a simple elliptic result: 

Lemma 7.3. For all I G Hi{QY, t G [0,r], there exists a constant C > and (f>{l) G i?5(0) 
such that {an)l{t)(l>'^ ,j = I in fl and 

110(0111 <qi^ll^^,,. (7.6) 

Proof. Let ■ip{l) be the solution of (7.2). Since tjj is linear and continuous from H^{ny into 
H^{n) and from L^{ft) into H^{ft), by interpolation, we have that ip is linear and continuous 
from H^{Q.)' into H^{Q.). We then see that ^'(Z) = {aK,)l'4'{l),j satisfies the statement of 
the lemma. □ 

For p e H2{n), we define the linear functional on X{t) by ((are)j(t)<^%j ,p) i , where 
(p G X{t) = {ijj G H^{n)\ {a^yii^'^^j G i/^(0)'}. By the Riesz representation theorem, there 
is a bounded linear operator Q{t) : H^{^) X{i) such that 

G X{t), {{a^)l{t)^\j ,p)i = {Q{t)p, ip)xit). 

Letting = Q{t)p shows that 

\\Qm\x(,)<C\\p\\^^^^^^ (7.7) 
for some constant C > 0. Using Lemma 7.3, wc see that 

V/ G Hi{ny, {I, p)i = {Q{t)p, m)xit), 

and thus 

\\p\\^i^^^<C\\Qm\x,, (7.8) 

which shows that R{Q{t)) is closed in X{t). Since Vj,{t) n X{t) = R{Q{t))^, it follows that 

X{t) = R{Q{i)) ®x(t) Vv{t) n X{t). (7.9) 

Our second Lagrange multiplier lemma can now be stated. 

Lemma 7.4. Let £(<) G X{ty he such that £,{t)ip = for any ip G Vv{t) n H^{n). Then 
there exists a unique q{t) G H^{Q), which is termed the pressure function, satisfying 

V^GX(t), Z{t){v) = {{a,)y,^,q{t)).. 
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Moreover, there is aC > (which does not depend ont ^ [0, T] and on the choice ofv^ C't ) 
such that 

Proof. By the decomposition (7.9), for (fi G X{t), we let (fi = vi + V2, where vi G Vv{t) fl 
H^{Cl) and V2 € R{Q{t)). Prom our assumption, it follows that 

£(i)(<^) = = {m,v2)x{t) = {m,v)x(t), 

for a unique ijj{t) G R{Q{t)). 
From the definition of Q{t) wc then get the existence of a unique q(t) G (Q) such that 

The estimate stated in the lemma is then a simple consequence of (7.8). □ 

8. Existence of a solution to the linearized smoothed k-problem (4.2) 

In this section, we prove the existence of a solution w to the linear problem (4.2), con- 
structed as the limit e ^ 0. 

The analysis requires establishing the regularity of the weak solution. Note that the extra 
regularity on uq is needed in order to ensure the regularity property for w, q, and their time 
derivatives as stated in the next theorem, without having to consider the variational limits 
of the time differentiated penalized problems. 

Theorem 8.1. Suppose that uq G H^^-^{Cl) and CI is of class C°° . Then, there exists a 
unique weak solution w to the linear problem (4-2), which is moreover in L'^{Q,T;H^^-^{Cl)). 
Furthermore, 

diw G L^{0,T; H^^-^-'^'in)) n L°°{0,T; H^'^-^-^\n)), i = 1,2,3,4, 
diq G L^(0, T; H^^-^-^'{Q)) n L°°(0, T; H^°-^-^'{fl)), i = 0, 1, 2, 3 . 

Proof. Step 1. The limit as e ^ 0. 

Let e = ^; we first pass to the weak limit as m ^ oo. The inequality (6.2) provides the 
following bound, independent of e: 

/ -\\{a,)iw,\j \\l + \w, ■ n,{f].)\i + \\w,\\l dt < N{uo) 
Jo ^ 

which provides a subsequence {w^ } such that 

w^^w in L'^{0,T;L^{n)) , (8.1a) 

(a«)>i^„ - (a,)>\, in L^^g, T; ^^(j^)) ^ (8.1b) 

mi 

■ n«(??«) w ■ njrj^) in L^O, T; H\r)) . (8.1c) 

The justification for w ■ ^^(t]^) being the third weak limit in (8.1) comes from the identity 
{anYiWe^jj = diy{we o f]i^^){f]K.) and the fact that is the normal to f]K,{Cl). 

Moreover, since (6.2) also shows that \\{aK)i'w'^j_ ,j \\L^(o,T;L'^(n)) — > as m — > oo, we then 

have \\{aK)iW\j \\L^{o,T;L^{n)) = 0, i.e. 

{a^)iw\j = in L2(o, T; L\ft)). (8.2) 
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Now, let US denote u = w 077^^1 so that thanks to (8.2) and (6.9) we have 

divu = in ?7K(il) , (8.3a) 

curlu(ry«) = curluo + / B(Vw„, Vu) in i?"^(0). (8.3b) 
Jo 

By proceeding as in Step 3 of Section 6, the trace regularity {u ■ nK){^K) € -^^(0,T; iJ^(r)) 
and the system (8.3) then yield 

ll"'llL=(0,T;Hi(a))^^("0), 

where N{uq) is defined in (6.1). 

Step 2. The equation for w and the pressure. 

Now, for any y G L^(0, T; iJ 2 (fj)) and I = {aK)\y^,j, we see that for a solution almost 
everywhere on (0, T) of the elliptic problem 

{^K)l[J;^'^{a^fi'P,k],j = / in O 

= on r, 

if we let e* = J~^{aK)i^-,k-, and set v = y—e, we have that e and v are both in i^(0, T; Hi{Q)), 
with 

pT pT 



2 



/ [\\e\\\ + \\v\\W<C I \\y\ 
Jo ^ ^ Jo 

Since (a«,)>%j = in 1,2(0, T; L^{n)), we infer that (a^)>j,j = -[(a,.) e i^(0, T; if ^ (f])), 
and that 

(J« u!t,e)3 = ([(a«)j]t«;\j,<i?)o- 
But V also satisfies the variational equation 

/ {JKWet,v)3+K [tWe • n„(^K)),i; • n„(^K)]i = a / [ig^ •««(??«), u • n„(?7K)]o , 
Jo Jo Jo 

leading to 

lim / {J^Wet,y)3= {[{aK.)i]tw\j ,ip)o+(T [ig?? • n«(^«), u • n«(^«)]o 
«^Ojo Jo Jo 

Jo 

We then see that as e — > 0, 



By standard arguments, we infer that Wet — ^ Wt in ^^(0, T; if 2 (O)'). This ensures that 

w £ C^{[Q,T]]L'^{yi)), and the condition w^{Q) = uq provides w(0) = uq. Furthermore, we 



also have for any </> G ^^(0, T;H^ (Q)) such that (aK)i^%j = in (0, T) x fJ, the variational 



equation 

r-T 



/ (JnWt,4>)3+K [w ■ nf,{r]^),(j)- nf,{r]K)]i = cr [ig?? • ««(/?«), ^ • nK(7?K)]o • 
Jo Jo Jo 
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Next, since wt G L^(0, T; -ff ' (f2)'), the Lagrange multiplier lemma 7.2 shows that there 
exists q e ^^(O, T; {Q)') such that for any (j) e L^(0, T; ifi (O)), 



Jo Jo 

~ lo ^'^'^""^^^i'^^'^h ^ [Lgr]-n^{f]^),(l)-n^{fj^)]o . (8.5) 

Now, if we have another solution w e L^{Q,T; (U)) such that w{Q) — uq and Wt e 
L^(0,r;il2(O)'), we then see, by using w — ■iy as a test function in the difference between 
(8.5) and its counterpart with w, that we get w — w = 0, ensuring uniqueness to the weak 
solution of (4.2). 

Step 3. Regularity of w. Wc can now study the regularity of w via difference quotient 
techniques. We will denote Rl = {x e R^] 2:3 > 0}, 5*0 = 5(0,1) n {x e R^\ X3 = 0} 
and B+{0,r) = B{0,r) n . We denote by 6* a diffeomorphism from B{0, 1) into a 

neighborhood F of a point xo€T such that 6»(B(0, 1) n M^) = Vnn, with det = 1. We 

1 

consider the smooth cut-off function tp{x) = e i"''^" 2 if a; € B{0, i), and ipix) = elsewhere, 
and with the use of the test function [D^hbpDhiw o 0)]] o 6''^ G L'^{0,T; Hi {Q)) in (8.5), 
with h = \h\ea{a = 1,2), we obtain: 

Ii + Kl2+h=cr / [LgT]-nK{fjK),[D-h[ipDh{wo9)]]o9~^ -nKifjKjjo , 
Jo 



with 



ii = [ {J.wt,[D_h[iljDh{woe)]]oe-^)3, 

Jo 



-T 



/2= / [d{wn,{n^)),d{[D_h['<pDh{wo9)]]oe-^ -n^inMo. 







- / {q,{a^)i[D_h[i>Dh{wo9):^]]o9-'). 
Jo 

For 7i, we simply have 

h =\\V^wo 6{t)\\l2i^B+(o,i)) -\\\/^uoo 9\\l2^B+{o,i)) 

+ [ {Dh[M9)]wto9,^Dh{wo9)). 
Jo 

>\\Vi^woe{t)\\l.^s^^o,i))-NM- [ \\wt\\^. \\Df,[M9)]i;Dh{woe)\\s 

J 

>||V^«;oe(f)||2 ( C;,iV(t.o)-^ f \\^R'Dt,{wo9)\\, (8.6) 

Jo 

where we have used (8.4) for the last inequality, and where the choice of 5 > will be made 
precise later. 
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For I2, we have, if wo define in -B+(0, l),W = wo9 and N'^ = nK,{riK){6), 
h = r I ^[W- Ar«],„ [D-h[^DhW] ■ N%is 

Jo JSo V'^O 

= r [ ^[DhW-N%^[i;DhW-N%0 
Jo JSo V°'0 

+ r f m^[W-N%„]-^[DHW-N%„][i,DHW-N%fs 

J a J Sn 



+ 



/' / [^[W ■ N%^]{. + h)[i^DHW ■ DhN%0 , (8.7) 

Jo J So 

where Ga^ = 6, a -0,0 and ao = det G. 

In this section, we will denote by || • \\s^e and | • \s,d@ the standard norms of iJ*(9) and 

For the first term appearing in the right-hand side of the second inequality, we have 



ao ^/ao 



and thus, since ip, and are chosen smooth, we infer that 

Jo JSo Va Jo 

- N{uo). 

The other terms in (8.7) are easily estimated leading to the estimate: 

l2>C r \\^/lpDhW-N^\\lso-N{uo). (8.8) 
Jo 

Concerning J3, we have 

h = - j\q,M[D_h[i'Dn{W)Y,joe-'])i, 
with = [V(77k o 6*)]^^ Now since (5„)^W*,j = 0, we obtain 

- D-hiihYii- + h)]DhW\j , 

and thus 

I/3I <c£ M^i^^yi m,W\j II .,5^(0,1) + 11^ Vi^D,W\j II .,^^(0,1) 

<CsN{uo)+sj^ iyV^^?.W^'.il|,B+(o,i)' (8-9) 



28 



D. COUTAND AND S. SHKOLLER 



where (5 > is arbitrary. Now, let be a smooth domain included in -B+(0, 1) and containing 
B+(0, i). The inequalities (8.6), (8.8) and (8.9) yield 



(8.10) 



/ \y/^DhW-N%aQ<C,N{uo)+6[ \\y^D,,W\\l ^. 
Jo Jo ^' 

We now define in -B+(0, 1) 

div^^oe W = div(W o d~'^ o f]~^){fjK o 0) = div(u)(^K o 
curlj^^oe W = curl(M)(7?K o 6*). 
Thus, (8.3) translates in B+(0, 1) to 
divfj^oe W = 0, 

[curl^^oe W]{t) = [curluo] 06+ [ B{Vu^, Vu){f]^ o 9), 

Jo 

= [curl«o]oe + / B{Wu^,\7W{fj^oe)-'W{fi^oe)-^){fj^o0), 
Jo 

and thus 

dWfiMV^DhW) = -^Dh{b^)i W\j {■ + h) + MDhW\ (8.11a) 

[cmlfj^oeiVi^DhWm = R{W) + [ B{Vu^{fj^ o 0), \7[^DhW][V{fjK o 0)]-'), (8.11b) 

Jo 



with 



\R{W)\\l Q < CN{uo). 



With the trace estimate (8.10) and the control of W in L'^{0,T; H2 (6)), we can then infer 
as we did in Step 3 of Section 6 that 

/ WViiDhWWl e < N{uo)+CJ [ \\ViiDhW\\l q, 
Jo ^' Jo ^' 

and thus with a choice of 6 small enough, 

[ \\ViiDhW\\lQ<C^N{uo), 
Jo ^' 

yielding 



/ \^/iiDhW\lQQ < N{uo). 
Jo 



10 

Since this estimate is independent of h, we get the trace estimate 

\Vi'W\lQQ<C,N{uo), 

and thus with this trace estimate and the div and curl system (8.11), still with arguments 
similar as in Step 2 of Section 6, 

r \^/i^W\\Q<C^N{uo). 
Jo ^' 
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By patching together aU the estimates obtained on each chart defining ^l, we thus deduce 
that 



T 

2 



||w||i < C« N{uo). (8.12) 

10 

Now, for the pressure, we see that for any y € (t) = {(j) e {an)'l{t)4i^ ,k e (fi)'}, 

for if a solution of the eUiptic problem 

{a.ma.tv\k],j = {a.)Ht)y\k in (if')'(O) 
= on r, 

we have by interpolation that cp G Hi{n). If we once again let e = {di^)'^ip,k, and set 
V := y - e, we have that e G Hi{n), v € V{t) = {0 S Hi{n)\ (a^)f(t)0%fc= 0}, with 
||e|| 1 + llt^ll 1 < C||y|| 1 , , • Now, by proceeding in the same fashion as in Step 2 above, we 

see that thanks to our decomposition and the regularity (8.12), Wt G L^{0,T;Xi{ty) with 



Jo 



*ll^i,,<A^(uo). 
By the Lagrange multiplier Lemma 7.4, we then infer 

\\q\\l < N{uo). (8.13) 

/o 

Next, by using D-hDhitpD-hDhw] as a test function in (8.5), we infer, similarly to how we 
obtained (8.12), that the estimates (8.12) and (8.13) imply that 

i-T 

hll 7 < N{uo). (8.14) 

2 

We now explain the additional estimates employed for this higher-order differencing. We 
need the fact that independently of any horizontal vector h, there exists a constant C > 
such that for SuppV* + ft. C 9, we have that 

V/Gffi(e), ||V?^„/||i,e<C ||/||3,e, 

V/ e H^Q), II V^I>./||^i(e), <C 11/11 i,e- (8.15) 

The first inequality easily follows by interpolation. For the second one, if / G L^{Q) we 
notice that for any (j> G -ff^(G), since the difference quotients are in an horizontal direction, 

/ ^f4,Dhf ct>= I ^fD^h<P+ I D_hVi^f^{- - h) 
Je Je Je 

<C||/||o,e||0||i,e, 
which shows that there exists C > such that 

V/ G L^e), \\V^Dh.f\\„i(oy < C ||/||o,0. 
By interpolating with the obvious inequality (for some C > 0) 

V/ G H'{e), \\Vi^Dhno,e < C ||/||i,e, 

we then get (8.15). 

Now, the pressure solves the elliptic equation 

Ap= -(u„)',jU.',i in ^^(n), (8.16a) 

p = -[cTAgT] ■ n^{f]^) ««(??«) + kAo{w ■ ri«(?7«)) n„(?7„)](?7^^) (8.16b) 
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Using the same change of variables that provides the pressure estimate (18.4), and using 
the eUiptic estimates for coefficients with Sobolev class regularity as in [9], we find that 



I h\\\ < iV(Mo,sup|w;«|4, / 
Jo ^ [o,t] Jo 



\w\\\ 



where the right-hand side is defined in (6.1). Therefore with (8.14), 



I \m\i 

Jo ^ [0,1 



2 < N{uo,snp\\w^\\4,). 



Higher-order regularity results follow successively by appropriate higher-order difference 
quotients, leading to, for n > 1, 

/* Ikll^+i + t IkllL 1 < CN{uo, sup \w,\n+2). (8.17) 
Jo ^ Jo ^ [0,t] 

Now, since Wt = — {a'^)lq,j in f2, we then infer that for n > 2, 

/ \\wt\\l_^<N{uo,SU^\w^\n+2). (8.18) 
Jo ^ [0,t] 

and thus in [0,t], 

Ik(i)||l3.5 < ||mo||i3.5 + Vt iV(Mo,SUp \Wn\l7)- 

[0,t] 

By Lemma 3.1 (for the smoothing operation given in Definition 2.2 on Ct), we have that 

||«'(i)||i3.5 < ||«o||i3.5 + \/i A^o(wo, C°), (8.19) 
where we use to denote a fixed (nongeneric) constant which depends on k. 

9. Existence of a fixed-point solution of the smoothed k-problem with 

SURFACE tension 

Lot A : {w G -Bo) w, with w a solution of (4.2). By the relation (8.19), we see that if 

we take e (0, T) such that 

^/t'^No{uo,C°)<1, 

then 

A(CtJcCt„. (9.1) 
We now prove that A is weakly lower semi-continuous in Ct^ . To this end, let ('tw")^Q be 
a weakly convergent sequence (in -L^(0, Tk; H^^-^{V,))) toward a weak limit w. Necessarily, 

By the usual compactness theorems, wc have the successive strong convergent sequences 

fj"-^fj inL2(0,T,;i/i2.5(^))^ 
(??")« inL2(o,T,;Fi2-5(f^)). 

Now, if we let = A{w"), we obtain from the stability of Ct^ by A and (8.18) the following 
bounds: 

''^ M\\lo.5<CN{uo), 



sup ||w„||l3.5 < 2||uo||l3.5 + 1- 
[0,T] 
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We thus have the existence of a weakly convergent subsequence {w"^"^) in the space L^{0, T^] H'^''^-^{Q)), 
to a limit / € Ct„. By compactness, from our bound on w", 



Prom the strong convergence of {rf'Y, we then infer from the relations (a"), J = in 
il, that 

(a,),;./,} = inn. (9.2) 



Moreover, we see that 



with p the solution of 

Ap = -{u^)\jP,i mf]^{n), 
p = -[crAgT] ■ nn{f]^) + kAo(u) • n^{fi^)]{f]~^). 

From the relations (8.5) for each n, we see from the previous weak and strong convergence 
that 

' ^ ^ ^ ^ ^Jl 



Jo Jo 

- {<l,{aKyi(f>\j)i =(T [Lgri-n,^{f]^),(f>-n^{f]^)]o , (9.3) 
Jo Jo 

which together with (9.2), and the fact that I G Ct^ implies that I = A{'w). By uniqueness 
of the limit, we then infer that 

^ w in L2(o, T«; H^'^-^iO)), 

By the Tychonoff fixed-point theorem, we then conclude the existence of a fixed-point 
w = w'm the closed convex set Ct^ of the separable Banach space i^(0, T^; H^^-^{Q.)). This 

fixed-point satisfies the smoothed system (4.1), if we denote tj = Id-|- / w and u = wori'^~^. 

Jo 

It is also readily seen that w, q and their time derivatives have the regularity stated in 
Theorem 8.1. □ 



10. Estimates for the divergence and curl 
Definition 10.1 (Energy function for the smoothed K-problem). We set 

Erit) = ^ \\d!^r,mh-k + Wvtuml and El^it) = ^ \\d^m\L-k + htmiml ■ 

k=0 k=0 

We use Ef^{t) to denote the energy function when the dimension is clear. 

We use these energy functions to construct solutions for the Euler equations. The increase 
in the derivative count from the 2D case to the 3D case is necessitated by the Sobolev 
embedding theorem. We will show that solutions of the K-problem (4.1) have bounded 
energy E^lt) for t G [0,T] when T is taken sufficiently small, and that the bound is, in fact, 
independent of k; as such, wc will prove that the limit as k ^ of the sequence of solutions 
to the K-problem converges to a solution of the Euler equations. 

Our estimates begin with the following 
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Lemma 10.1 (Divergence and curl estimates). Let n := dim{il) —2 or 3. Letting Li = curl 
and L2 = div, and let r]o ■= r]{0) and 

Mo := P(||Mo||2.5+n, |r|4+n, \/K||'Wo||l.5+3n, \/K|r|i+3n) 

denote a polynomial function of its arguments. Then for j = 1, 2, 

sup ||^/^LJr/(^)||i5+„ + ^ ( sup IIMtSWIlf.S+n-fc + / llx/^^i^NllL+n-fc ) 
te[o,T] \*e[o,T] Jo J 

<Mo + CTP{ sup E^{t)). 
te[o,T] 

Proof. In Eulerian variables, equation (4.1b) is written as + u\i (mk)' +p,i = 0, where the 
transport velocity is the horizontally smoothed vector u^. Taking the curl of this equation 
and using the formula (curlw)* = eiji~u^ ,j with e^j^ denoting the permutation symbol, we see 
that Eijkldtu'^ ,j +u'',ji + u'',i ] = 0. Thus, defining the bilinear form B*(Vu, Vu^) = 
iijku'^ ,1 {ukY , we can write the vorticity equation as -^curlu = _B(Vu,V?iK). (When 
the transport velocity is divergence-free, then B is the familiar vortex-stretching term.) 
Composing this equation with 77^, switching to Lagrangian variables via the chain rule, and 
integrating this from to t, we have 

Sijkv'', r a^^'j =cmluo+ BaST)dT, B^^ := eijkJ^^v'' ,r aJiiv^Y ,maKf- (10.1) 
Jo 

This is the time-integrated Lagrangian form of the vorticity equation. We will need to 
space-differentiate this equation once more for the estimate on curl 77. Hence, 

eijkSv^,r ajj = Vcurluj) + Sikjv'',r Va^ + / VBa,{T)dT . (10.2) 

Jo 

We begin with the estimates for the case that n = 2; we set E^{t) = E'^^{t) and proceed 
with the estimate for curlr/. Using that Vu'^,ra„J = 9t(Vr/'^,r OkJ) ~ Vr?'^,^ ^ta^J, we see 
that 

Sijkdt{Vr]'',raJj) = V cm\ + SijkVrj'^ ,r dta J, + Sikjv'' ,rV a J, + / VBa^{T)dT . 

Jo 

Integrating once again in time from to f yields 

Jo Jo Jo 

where 

+J-^v\r vj,m (Va4a«™ + a^lVa^f) + (V 7"^) v\r vj,m (10-3) 

and 

ata«7 = J-\a4a^^ - a«™a«j)^;«^, , (10.4) 
VJk = a4^\/r]J,s . 
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Since < C'||^^||s (and similarly for 77^), we will write (10.3) and (10.4) in the following 

way: 



Vi?a„ ^ J-^al\7vV\ + J-'al{VvfV^7j, 
V ~ a^V^r? , 

where note that we are not distinguishing between rj^ and rj or between and v in the 
highest-order terms. The point is that the precise structure of these equations is not impor- 
tant for our estimates; wc need only be careful with the derivative count appearing in these 
expressions. The power on each expression is merely to indicate the number of times such 
a term appears. 

Next, with the fundamental theorem of calculus, 



eufcVr?'',^a«J = Vcurl77* +£y7;V7?'',r- / 9ta«^ 

Jo 

so that 

V(curlr/'— f curlug) = Sijk 
Let 



Vr/' 



'',r / dtajj + 
Jo Jo 



i pt' 

Jo 



(10.5) 



F:= P{J-\a^,Wv) andFi := P{F,W^t],V^v) 
denote polynomial functions of their arguments. We then express (10.5) as 

Vcurlr?' ~iVcurl4-FV^77 / F + [ FV'^r]+ f f F (V^r? + V^t;) , 

Jo Jo Jo Jo 

and taking two more spatial derivatives yields 

V^curlr?' ~tV^curl4-K V^r? / F + V^r}[ Fi + V'^rjf Fi 

Jo Jo Jo 

+ ( [ I + [ ] [Fi{V^r] + V\) + FV^+ f I [FiV^i; + FVS] . 
\Jo Jo Jo J Jo Jo 



Since jl f; FV'v = - /J F, V^r? + F, V^r? 



t rf 



curl r?* ~tV^ curl 4-1- V^?7 / f -h (V^t? -h V^r?) / Fi 



Jo 

I I + / I [^1 ('^^^ + "^^^^ + {F + Ft) V\] + / / FiV\. (10.6) 
^Jo Jo Jo I Jo Jo 
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We use interpolation to compute HV"^ curl?7||o.5 — || curlr/Ha.s. We begin with the highest- 
order term: 



[\F + Ft)V^r, 
Jo 

f\F + Ft)V^V 
Jo 



< sup + 
te[o,T] 







<C sup WF + FtW^ 

4 te[0,T] 



< sup {\\F + Ft\\Loo 
1 te[o,T] 



7? + sup \\F + Ft\\L. 

5 *e[o,T] 

t 

<C sup ||F + ft||2 
te[o,T] 



V 



Jo 



Since \\F + Ft\\2 < C||F||i:,oo ||wt||2.5, by the interpolation theorem 7.17 in [1], 



Jo 



Ft)W\ 



< C sup lli^ll \\vth.^ 
0.5 te[o,T] 



Jo 



<CT sup ||i^||||7;t||2.5b||4.5 

4.5 te[o,T] 



The other terms have similar estimates in the if'^'^(f2)-norm, so that 

sup ||curlr?|||.5 <T||wo||i5 + i sup || / v\\l^+T sup Ht^tlk.slhllL • 
te[o,T] t'e[o,t] Jo *6[o,T] 



(10.7) 



By differentiating (10.6) once more in space, the same interpolation estimates show that 

sup ||V^CUrlr/||i5<T||^/;^Uo||5.5 + ^ / 1^111.5 + ^ sup |bt||2.5||\/^?7||5.5 • 

tG[0,T] Jo te[o,T] 

Next, we rewrite (10.1) as 

curl t; = curl Mo + £yfei^'',r / dtajj+ . (10.8) 

Jo Jo 

Using the fact that H'^{Vl) is a multiplicative algebra for s > 1, It follows from (10.3) and 
(10.4) that supjg[o,T] II curlt;(t)||2.5 < ||uo||4.5 + C^'-P(supj£[o,T] E^{t)). Differentiating the 
above expression for curl v yields 

curlt;t = eijkV^,r dtaiJj + B^^ + Sijkdtv'^ ,r / dtaj^ + Ba^ , 

Jo 

so with the fundamental theorem of calculus and our generic polynomial function F, 



cvivXvt ~ P(Vwo) + Vt;t / F+ ( Ft, Ft ^FVvt. 
Jo Jo 



(10.9) 



Again using the properties of the multiplicative algebra, we see that sup^gjo t] II *^^^l''^«(0lli.5 — 
P{\\uo\\4.5) + CTP(sup^g[o.T] E^{t)). From the time differentiation of (10.9), 

curlt;„ ~ VvtF + Wvtt I F 
Jo 

~ Vt;t(0)P(Vuo) + Vf„ f F + [ [FVvtt + FVvtVvt] (10.10) 
Jo Jo 

We must estimate the if°-^(r2)-norm of the three terms on the right-hand side of (10.10) by 
using interpolation. Let L denote the linear form given by L{w) = /o Fw. Then 

||i(«')||o<Co|| /*«;||o, Co=sup||P|| 
Jo [0,t] 
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Letting Fi :— P{J^^,af,, Vv, V^i], V^w), it is easy to check that 

IliHIll <Cl|| / W\\i, Cl =SUp||i^i||Loo. 
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By the interpolation theorem 7.17 in [1], 



Jo 



FVVu\\o.5<VCoVCl II / Vu\\l.5, 



[0,t] 



/' 

Jo 



so that by Jensen's inequality and Sobolev embedding, || Jg -FVutt||o.5 < C'^-P(supj£[o,T] 
All of the other timc-dcpcndcnt terms in (10.10) have the same boimd by the same inter- 
polation procedure. For the time t — term, interpolation provides the estimate 

||P(VMo)Vi;t(0)l|o.5 < CP{\\ua\\i.r.)\\vim\i.r, < CP(||wo||4.5)|k(0)||2.5 < Mq. 
The initial pressure q(0) solves the Dirichlct problem 



Ag(0) =io := K)kF, 



n. 



q{0) = bo := ^no}5o'^7?o%a/3 K + «AoK • AT,) on T . 
V 90k 

Since ||g(0)||2.5 < C'(||io||o.5 + jbob) < Mo, we see that sup^gjoT] II curlwtt(t)||o 5 < Mq + 

CTP{supt^^o,T]E.{t)). 

Differentiating (10.10) with respect to time, we see that cmlvut Vi'ttt Jq F + F'Vvtt + 
FVvt Vwt so that by the fundamental theorem of calculus 

FVvtt + FVvt Vvt = F(o)[Vwtt(0) + Vvt{{)) Vwt(O)] 

+ / [FWvttt + FWvtWvtt+FWvtWvtWvt], 
Jo 

so that 



\\V^cm\vttt\\lz< / \\V^F{{))[Vvtt{0) + Vvt{Q)Vvtm\L. 
Jo 

/ VKVvttt F + / [Fy/]iWvttt+FS7vty/KVvtt + FVvtVvt\/KVvt] 

Jo Jo 0.5 -^0 Jo 0.5 

(10.11) 

We repeat the interpolation estimates between L'^{fl) and iJ^(O) just as for the estimates 
for curlwtt; for example. 



/ F^/UVvttt 
Jo 



< VCoVCi 



0.5 



'KVttt 



1.5 



SO that by Jensen's inequality and Sobolev embedding. 



I / i^x/;^Vt;«t||o.5 < Ct SUpS, \\^/^Vtu\\l,r, . 

Jo [0,t] 



Thus, integrating from to T gives the estimate 



[ II / F^Vvttt\\i5<CTP{ sup E,{t)) [ \\V^vttt\\i5<CTP{ sup 
Jo Jo te[o,T] Jo «e[o,T] 



(*))• 
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The other time-dependent terms in (10.11) have the same bound by the same argument. 
The time t = terms require analysis of the elliptic problem for qt{0): 

Agt(O) ~ ti := V[P(Vmo) Vg(0)] +FVvt{0) in Q , 
qt{0) ~ bi := Q{drjo)d\o + Q{d7]a)duad''rjo 

+ kAoK(O) • N^) + kAo(wo • QidrjoJduoJ on T. 
By interpolation estimates (as above), 

^||^/^F(0)[V^;«(0) + Vt;t(0)V^;t(0)]||;;,</vT||P(V^^o)||ioc||MO)||?.5, 
Jo 

and since time-differentiation of the Euler equations shows that 

t^tt(O) =-VwoV(z(0)- Vgt(O), 
interpolation provides the estimate 

V^||ftt(0)||l.5 < V^IIV'uo V(7(0)||o.5 + v^IIVmo V\{0)\\o.5 + V^|kt(0)||2.5 

< Mo + \/K||il||o.5 + \/K|bl|2 

< Mo + \/K|bl|2 

where we have used the elliptic estimate ||(j'(0)||2.5 < Mo (from above) for both the second 
and third inequalities. (The remaining estimate for |bi|2 places the regularity constraints 
on the polynomial function Mq in the hypothesis of the lemma.) Because H'^{T) is a multi- 
plicative algebra, the bound for -v/K|bi|2 is controlled by the highest-order terms \/K|^t(0)|4 
and v^IuokIs < V^C'||mo||5.5- Now, 

\/K\vt{0)U < \/K||g(0)||5.5 < \/K||io||3.5 + \/k|&o|5 , 

and ||io||3.5 is bounded by P(||uo||4.5) while the highest-order terms in y^lbols require bounds 
on \/k||^o||7.5 and \/k||wo||7.5 With our definition of Mq, we see that 

KT||P(Vuo)||ioo||?;„(0)||f.5<Mo 

and hence \\^/KcwLlvttt\\o.5 < ^o CrP(supjg[o E^{t)). 

The proof that \\^/Kcmlvtt\\i,5 < Mo-|-CTP(supjg[o,T] is essentially identical. 

The divergence estimates begin with the fundamental equation a^^^v"^ ,j = 0. By taking 
one derivative of this equation and integrating-by-parts in time, we find that 

Vdiv7? = Vr?'„- fdta4+ [\dta,iyr)\j-Wa4v\j). 
Jo Jo 

Computing the if^-^(f2)-norm of this equation yields the estimate sup^gjg || div?7(t)||3 5 < 

Mo -I- CrP(supjgjo,T] Ei^{t)). The divergence estimates for v, Vt, Vu, y/KVu, and y/KVut 
follow the same argument as the corresponding curl estimates. 

In the case that n = 3, the estimates are found in the same way, with one minor change. 
Set E^{t) = E^^{t). The estimates for curl?7, which rely on Sobolev embedding, require 
greater regularity on vt- The estimate (10.7 becomes 

sup ||curlry||2 5<T||uo||L+* sup || / v\\l^ + T sup Ht^tlUlhllL > 
te[o,T] t'e[o,t] Jo te[o,T] 
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and similarly, 

sup ||v^curlr?||i5<T||v^uo||L+I^ / IMls + T sup ||t;t||3||v^v?||L • 
te[o,T] Jo «e[o,T] 

□ 

11. Some geometric identities 

We will usually omit writing dSo in our surface integrals, and for convenience we set 
a = 1. Let n* denote the projection operator onto the direction normal to //(F), defined as 
= 5j — g^^V^^jillffj where g is the induced metric on r]{T) defined in (1.4). The mean 
curvature vector motivates us to introduce the projection operator 11. In particular, we have 
the important formula 

-^/gHno'q = ^/gAgirj) 

= ^g'^'^uy . (11.1) 

where the last equality follows since {g^^g'^" — g°"'g'^^)'n^,/3 V'' w ^ncx = 0. For a vector field 
F on F, nF = [n- F]n, i.e., H = n O n. 
We let 

Q{dr,) = h{dji)/h{Vg), (11.2) 

denote a generic rational function where /i and /2 are smooth functions. We record for 
later use that that n = = |^;^^^;^| and that |a'^Ar| = ,/detg on F, as 

= {SjrSks - SjsSkr)v'',i'rf,sV'^,iV\2= \va \'^\V,2 f - [va -^,2]^ = detfif, 

where Sijk denotes the permutation symbol of (1, 2, 3). We will use the symbol Q to denote 
any smooth (tensor) function that can be represented as (11.2). 

Remark 9. The L°°-norm of the numerator of Q is hounded by a polynomial of the energy 
function, while the L°°-norm of the denominator of Q is uniformly controlled by (3.3a). 
Thus, the generic constant C which appears in the following inequalities may depend on a 
polynomial 0/ dot go- In particular, \\Q(dri)\\jj^ < C(det (7o)||^'(977)||i,oo. 

For a vector field F on F, • A?^ = • n + • (A^ - n) and 

\N-n\Loo<[ \nt\Loo= [ \Q{dr])dv\Loo <CtP{E^{t)), 
Jo Jo 

the last inequality following from (3.3a). If \UF\s < Mq + CP{E^{t)), then \F ■ N\s satisfies 
the same bound. 

12. k-independent estimates for the smoothed problem and existence of 

solutions in 2d 

All of the variables in the smoothed K-problem (4.1) implicitly depend on the parameter 
K. In this section, where we study the asymptotic behavior of the solutions to (4.1) as 
K — > 0, we will make this dependence explicit by placing a over each of the variables. We 

set Er^it) = El^it). 

Remark 10. The only difference between the 2D and 3D cases arise from the embedding of 
vt e Z/°°(0). In 2D, vt e H^-^{Q.) is sufficient, while in 3D, we needvt e H^{Q.). 
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The pressure function q can be formulated to solve either a Dirichlet problem with bound- 
ary condition (12.5) or a Neumann problem found by taking the inner-product of the Euler 
equations with a^N. We use the latter. 

Lemma 12.1 (Pressure estimates). With {v,q) a solution of the K-problem (4-1) 

wmwL + umL + Uuimi < cpiE^) ■ (12.1) 

Proof. Denoting by A, we define the divergence-form elliptic operator La and corre- 
sponding Neumann boundary operator Ba as 

LA = d,{J-'AiA{di), BA = J-'AiAlNjdi. 

For fc = 0, 1, 2, we analyze the Neumann problems 

La{d^q) = fk in 0. Ba{d^q)=gk on T (12.2) 

with 

/o = dtAl v\j go = -vt ■ V^n« 

/i = -LaS) - d^Al v\j -dtAl vl,j gi = BaS) - vt ■ (v^n«)t - vtt ■ V^^« 

/2 = -2LAt{qt) - LAttiq) - 9tAl v\j 52 = 2BAt{qt) + BAtt{q) - vm ■ -\/^"'« 

- 2dtAl vl,j -dtAl vlt,j - 2vtt ■ (-\/^n«)t - Vt ■ {\^h^)tt ■ 

For s > 1, elliptic estimates provide the inequality 

Wdhms < Cs[P{\\v\U.5)\\fk\\s-2 + \gk\s-3/2 + \\q\\o] , (12.3) 

where \\ ■ denotes the norm on [H^{il)]' . We remark that the usual ff* elliptic estimates 
require that coefficient have the regularity d''-'^{AlAi) € L°°{n), however (9^-i(A^^^') e 
L^{fl) is sufficient. See see [9] or the quasilinear estimates in [19]. 

As we cannot guarantee that solutions q to the K-problcm (4.1) have zero average, we use 
ll^llo < C'llo'lli and the elliptic estimate for the Dirichlet problem LA{q) = /o in with 
-q = Agfj ■ Uk + kAo{v ■ Hk) on F. Thus, < C(||/o||o + \A.gfj ■ Uk + kAo{v ■ nK)|o.5) < 
CP{E4t)). 

From (10.4), it is clear that II/0II1.5 + I.90I2 < C!P{E,^{t)); thus, from the elliptic estimate. 

Ml, < CP{E^{t)). (12.4) 

Next, we must show that ||/i||g.5 + \gi\l < C P{E^{t)). But /i ~ P( J-\ A, Vw„)([V5i]2 + 
V^q) so that with (12.4), ||/i||o.5 < C P{Ei^{t)), with the same bound for \gi\l, so that 
Il9t||i.5 < CP{E^(t)). Using this, we find, in the same fashion, that II/2II0 < C P{E^{t)). 
The normal trace theorem, read in Lagrangian variables, states that if Vttt G L^(ri) with 
W^l^tttij llo e L^(0), then Vtu ■ Vd^n^ e if~°-^(r) with the estimate \vttt ■ v^'t-kI-os < 
CP{Ejt)). Since |1 Tr(A V«„t)llo - \\Tr{3AtVvtt + SAtVvt + Attt'^iM < C P{E^(t)), 
and using the above estimates for q and qt, we find that I52I-0.5 < C P{EK{t)), thus com- 
pleting the proof. □ 

Our smoothed K-problem (4.1) uses the boundary condition (4.1e) which we write as 

gn^ = -^ffn • n^n^ - kAo({; • n^) , (12-5) 

where (we remind the reader) k > is the artificial viscosity, Aq = \/g^ ^ dai^/gag^^ dp) , n 
is the unit normal along the boundary T){t)(T) and is the unit normal along the smoothed 
K-boundary 77K(f)(F). 
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We begin with an energy estimate for the third time-differentiated problem. Although 
we are doing the estimates for the 2D domain O, we keep the notation of the 3D problem as 
well as terms that only arise in 3D when differentiating the mean curvature vector. Thus, 
when we turn to the 3D problem in Section 14, the modifications will be trivial. 

Lemma 12.2 (Energy estimates for the third time-differentiated K-problem). For Mq taken 
as in Lemma 10.1 and ^ > 0, solutions of the K-problem (4-1) satisfy: 

sup [\\vttt\\l + \vtfn\l]+ I Iv^^t^w-n^l? < Mo + rP( sup E^{t)) + 5 sup E^{t) 
te[o,T] Jo te[o,T] te[o,T] 

+ C sup [P(||t)t||i.5) + J"(l|t5|li5)+nil'7lll5)] + C'P(||^€tt||i.(o,T;i/-(a)))- (12-6) 
te[o,T] 

Proof. Letting A = a^, and testing df{Ji^vl) + df{A^q,k ) = with with dfv'^ shows that 

ri [ d^J^vDd^v^- r [ d^{Atq)dtv\u=- f [ d^,{^.qn,{f,.)) ■ d^^v dSo . 
Jo ^ Jn Jo Jn Jo Jr 

(12.7) 

Step 1. Boundary integral term. We rewrite the modified boundary condition (12.5) 
as 

V9 



qriK 



Hh + Hh ■ {hn — h)n + Hh ■ (n^ — h) — kAo{v ■ ««) hn . 



(12.8) 



We first consider the boundary integral on the right-hand side of (12.7) with only the first 
term on the right-hand side of (12.8): 



Jr 

T 



d'}{./^Hn' ori)dlv' dSo 



Qf{dft,dv)d^v'^dfv]^ 

la Jr 
I + II + III + IV . 



Qf{dfi,dv)dii] 



(12.9) 



The first term / on the right-hand side of (12.9) is given by 



-| T T 

^ Jr J Jo Jr 



where we use the notation /]q = /(T) — /(O). Since IljV^^,/}= {IljV^^),^ —lIj,0V^^ and 



n},/3 = Q'^i{df,)f,\^p with Q{dfi) defined by (11.2), for 5 > 0, 

-\l ^/l{n)d^,vlp)r^{W^d'iv^^) < -\\mu\l + 5\mu\l + (1 + Cs) \Q'f{df,) f,\^0 5; 



where the constant Cs depends inversely on 5. Since for any t G [0, T\ 



it follows that 



I<-7; sup \Uvu\l + Mo{6) + S sup E^{t) + CTP{ sup E^{t)) . 



te[o,T] 



te[o,T] 



te[o,T] 
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The second term // requires some care (in the way in which the terms are grouped 
together). Letting 



V,2 ■ dfv^i u 2 • dfv^2 



(12.10) 



we find that 

i-T 



(12.11) 



(12.12) 



// = / / det g-^ {dt det - det A"^ - det A^) 
Jo Jr 

n-(dets'"5)t det^l^ -det5"5(det^2 + det>l^)+ / det 5" ^ det 
Jr 

For a = 1, 2, let Va = ^,c« • d^v; thus Va,f3 = + ??*,a/3 ^tt so that 
det^^^ = det(F„,^ -??\c,/3 vlt) 

= detVo„0 - det{fj\^0 vit) + Pij{d^fi)vivi, + P^j {8^)^14, a • 
With A = Co{{dV), detdV = Al^°'Va,p. It follows that 

J detg-^detdV = - J (det 5-^),/3^^"^,a , 
as A^",/} = since A is the cofactor matrix. Hence, 

^ det ~g-^ det ^1 = ^ Pijid^mMt + P^j{d'mMt,a , 
so that 

II< r [ Qtf{dv,dv)d^v:^d^v^^0+ [[PiAd^)iltiit + m9'mA,a]\ . 
Jo Jr Jr Jo 

By the fundamental theorem of calculus and Young's inequality, for 5 > 0, 

J^[prjid'fj)vUvi,,^]{T) = j[p^^{d^)vim 4,, (t) + £ mo'vyu], atvi,,^ (t) 

< Mo(5) + 5||5«(T)||f.5 + T ( / sup \{P^^{d^)vl,)fdx\ \\vi{T)\\i,. 

Since [Ptj{d^fi)vlt\ G L°^{Q,T- L^{V)), wc conclude that 

II<Mo{5)+5 sup E^(t) + CTP{ sup E^{t)) . 

t£[0,T] te[o,T] 

A temporal integration by parts in the third and fourth terms on the right-hand side of 
(12.9) yields 



III + IV = 



= -l I iQtfi^V, dv)vi,0 +Qndfj, dv)]^ vi,0 

+ Imfidv, dv)vi,0 +Qt{df,, dv)]vl 
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which has the same bound as term //; it follows that 



Jo Jt 



= sup \mu\\ + MQ{5) + 5 sup E^{t)+CTP{ sup E^{t)) . 
^ te[o,Ti te[o,T] te[o,T] 

(12.13) 

Remark 11. The determinant structure which appears in (12.11) is crucial in order to 

obtain the desired estimate. In particular, the term det Ai is linear in the highest-order 
derivative ddfv rather than quadratic (as it a priori appears). 

There are three remaining boundary integral terms appearing on the right hand side of 
(12.7) arising from (12.8); the terms involving k are 

- K {[dt{v ■ n«), dfv • n„] + 3 [d^{v ■ h^),dlv ■ dtfi^,] + 3 [dt{v ■ n«), d^v ■ n«] 
Jo 

+ [v-n,,d^tv-d^,h^]^) . (12.14) 

The first term in (12.14) provides both the energy contribution jj" \^/Kd^v ■ n^l^ as well as 
error terms. We start the analysis with the most difficult error term, 

k( {[v-dd^n^,d{h^-d^v)]o, (12.15) 
Jo 

whose highest-order contribution has an integrand (modulo L°° terms) of the form d'^s/Hvnn 

\fKdvttt- 

With fi^ = {dxf}^ X dif]^)! ^ =: Qidf]^)^ Q given by (11.2), the highest-order term in 
ddfuK is Q{dfjK)d'^d^VK, so that with TZ^ denoting a lower-order remainder term, and using 
(3.3b), we have that 

- K I [i- ddfn^, d{n^ ■ dfv)]Q < C sup \P{v, df)^)\L-^ / \^/K^^v ■ n«|i \^/lld1vi^\2 + Ui 
Jo te[o,T] Jo 



< C sup |P(u,9^«)|l~|v^9(^w • n«|i2(o,T;Hi(r))|\/K5t^5«|i2(o,T;H2(r)) -|-7^l 

te[o,T] 

r n ^ 

<Cs sup \P{v,^f|K)\Loo\\^/KV|^^^\\L2^Q^T.H2.5^Q,)) ' - 

te[o.T] 

<Ma + CTP{ sup E^it)) + \\y^v^tt\\hroT-m-Hn))+S sup E^{t) , 
te[a,T] te[o,T] 

where TZi also satisfies TZi < CrP(supjg[o t] E^it)) + (Jsupj^jQ Ei^{t). The second term 
in (12.14) has a highest-order contribution with the same type of integrand, and its analysis 
(and bound) is identical. The third and fourth terms in (12.14) are effectively lower-order 
by one derivative with respect to the worst case analyzed above. 

Next, we estimate df{Hn ■ {fiK — fij} dfv. Since = Q{f]K) and since — < 

sup^ \dQ{dr]K)\ ■ |5?7k — dr]\, then our assumed bounds (3.3) together with (2.6) imply that 

|fi«-fi|ioo < C s/^\P{dfi,d'^f})\Loo\fi\2.5 < CKP{E^{t)). (12.16) 

Similarly, 

\dn^ -dh\Loo<Cy/K \P{dfl, d''f,)\Loo |?7|3.5 . (12.17) 
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Also by (2.6), for k = 1,2,3, 

\d!^h, - d^h\L^ < cV^\p{dfj,d^)\L^ IdfvU.^ , (12.18) 

and 

\d^v^-d^v\o < CV^\vtt\i.5- (12.19) 

Taking three time-derivatives of formula (11.1), we see that the highest order term in 
df{\/gHn) is Q{dfi)d'^Vtt- Thus, the highest-order term in the integral 

r [ dU^Hnrn% « - h^)d!v^ 

JQ Jr 

is estimated using an integration by parts in space. The highest derivative count occurs 
when the tangential derivative is moved onto the Vttt term giving us 

r [ Qt{dv)dvlK {< - n^diltt <C r \P{dfj)\L^ \vu\i |n« - nk- \d^vU 
Jo Jr Jo 

<C I \P{dfi,dff)\L<» \ f1\2.5 \vtt\i \^vut\i 
Jo 

<CsTP{ sup E,{t)) + S f \V^Vttt\l, 
te[o,T] Jo 

where (12.17) is used for the second inequality. If instead, integration by parts places the 
tangential derivative on fi^ — n, then (12.17) provides the same estimate for this term. The 
other terms are clearly lower-order. 
Thanks to (12.18), 



Jo Jr Jo 



<CsTP{ sup E^{t)) + 5 sup E^{t). 
te[o,T] te[o,T] 



We next consider the integral 

i-T 



= / ^/§Hh ■ d^n^ (riK - n) ■ Vttt + / / y/^Hn ■ dfih^ - n) ■ vttt + 7^2 
Jo Jr Jo Jr 

=:I + II + 1l2, (12.20) 

where 7?.2 is a lower-order term. For term /, we use the estimate jn^ — n|L°° < C* k|?7|3.5 One 
^/k goes with dfh^ and the other ^/k goes with vttt- Thus, \I\ < CTP(sup(g[Q E^{t)) + 

^SUPtg[o,T]-E^«(*)- 

To study //, we set / = ^/gHh ■ n^, and consider the term Jp / httt ■ vttt- We expand 
V into its normal and tangential components: set Tq = 77,^, so that 

V = v'^ T + n, where t = {v ■ To) Ta and v" = v ■ h. 

Then 

Vttt = Vfti-T + Svl^Tt + SvjTtt + v'^Tttt + Vfttf^ + Si^r*^* + 3wf ntt + v^-fittt ■ 
The most difficult term to estimate comes from the term v^^^t, which gives the integral 
lo frf^ttf'^'^ttf 
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First, notice that hut ■ t is equal to —h ■ Tut, plus lower order terms that have at most 
two time derivatives on cither h or r, and h ■ ntt = h ■ dpvtt for /? = 1 or 2. Next, the k 
problem states that wj^j = {A^q,^ )tt, where recall that 

vl = a„ . 

We have the formula 

r?\/3 A'^dkqtt = J^dfjqtt (no sum on /?) for /3 = 1, 2 , 

where 

(In the case that k = 0, = J = 1.) Using this, we see that the highest-order term in our 
integral is given by 



// 



Jf {n ■ dpvtt)d0qtt ■ 



Second, write qtt as 
qtt = - 



^ [A§(r?) • n + Ag(r?) • (n« - n)] + kAo{v ■ n«) 



(12.21) 



(12.22) 



We begin by substituting the first term on the right-hand side of (12.22) into (12.21); the 

highest-order contribution comes from dpdf/\g{ri) = Q{dfj,dfjfi)g'^'^h ■ div^^i^p. Integrating 
by parts with respect to dv, the highest-order term in our integral is given by 

/ / Q{d'n,d'ni^)f{n-vt,i^0)g'^''{h-vtt,u0). 
Jo Jr 

Letting := Q{dfj, dfj^) f hi hj = Q{dfi, dfiK)d'^r], integration by parts in time yields 

<CTP{ sup K(t))+Mo + C sup IGtUoo llvtlll^ 

te[o,T] te[o,T] 

<Mo + CTPi sup E4t))+C sup [P(||i)t||i5)+^(ll^ll3.5) + ^(ll'7ll4.5)]- 
te[o,T] te[o,T] 

(12.23) 

For the second term on the RHS of (12.22), the highest-order term gives the integral 

I fQ{dfi,dfi^){h-vtt,0)r''^vt,^u0-{h^-h)<CTP{ sup E^{t)) + 5 [ W^/^vtWh, 
Jr te[o,T] Jo 

where we used {hn — h\L°° < C k|7?|3.5 again. 

For the third term on the RHS of (12.22), the highest-order term gives the integral 



« r / fQ{dfi,d7u){ 
Jo Jr 



h ■ d vu) (h^ ■ d vu) 
<Mo + CTP{ sup E^t)) + ||v^i)tt||i2(o,T;H2.5(a) • 

*6[0,T] 

We have thus estimated the integral /p df{Hh ■ {h^ — h)} dfv. The remaining integral 
So Sr 9f{Hh ■ {hfz — h) dfv has the same bound. 
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Step 2. The pressure term. Wc next consider the pressure term in (12.7): 

- r [ d!{A^q)d^v\k = - f j dfv\k [dfA^q + Sd^AUt + SdtAUtt + A^dU] 

Jo JQ Jo JQ 

=: I + 11 + III + IV . (12.24) 
We record the following identities: 



dtA'y 








(12.25a) 






-A>;Af)dtv:,s- 


^P^{J-\A,Vi,), 


(12.25b) 




= Jk. ^{A-rA'l 


- AiAt)d^ff,,s 


+P^{J-\A,\/v^)dtvU . 


(12.25c) 



With (12.25c), and := J-^{Af.Al - A'^Af)q term / is written as 

1= r [ [f:^d^v\kdX,s+d!v\,dtdi,iP^HJ;^\AVv,)] 

Jo Jn 
=: la + h- 

We fix a chart 6i in a neighborhood of the boundary F and let ^ = y/al^ where once again, 
we remind the reader that {ai}^^^ denotes the partition of unity associated to the charts 
With la denoting the restriction Ia\uii where f/; fl f2 = ^^((0,1)^), and letting 
p := pi and 9 := 9i, we have that 

r [ f:H&)d?v\k{9)^p^HP*h^d^y\sid)+fmitttGi^,^0vtt 

Jo J(0,l)3 
=: 2^0 ^ + Ia 2 1 

where G{^, V^) is a bilinear function which arises when the gradient acts on ^ rather than 
vtt- The term Xai is the difficult term which requires forming an exact derivative, and this, 
in turn, requires commuting the convolution with f^l^. We let 

v{e) = p*h 

so that using the symmetry property (2.7), we see that 



^ai = f f P*h [f:Hd!v\k {eMv^,s (e) 

Jo ^(0,1)3 
Jo J(0.1)3 



'0 ^(0,1)3 



Since is symmetric with respect to i and r and k and s, we see that 

lau = \f I frHm [d^V\k {9)d^V\, {9)] 
^ Jo i(0,l)3 

= -\f \ dtn\{9)d1V\^ {9)d1V\, {9) + \ [ f?H&)d?V\k {9)d^V\, (9) 

^ Jo J(0,l)3 ^ ^(0,1)3 



T 

2Trr 
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We sum over / = 1, L. The spacetime integral is bounded by CTP(suptg[o,T] EK,{t)). For 
the the space integral at time t = T,we employ the fundamental theorem of calculus: 



n 

T 



tWL" 



<\\Vu{T)\\\mh + \\Vu{T)\\\\\ j I/, 

JQ 

<lk«(?^)ll?:'5 lk«(T^)liy'll9o||2 + CTP( sup E,{t)) 

te[o,T] 

<\\\MT)\\l, + CmMT)\\y^\\q4T + CTP{ sup E,{t)) 

<5\\vu{T)\\l^ + C(5)\\q4l + CTP{ sup E,,(t)) 

/e[o,T] 

< 5 sup E^{t) + Mo{S) + CTP{ sup E^{t)) , 
te[o,T] te[o,T] 

where Young's inequality has been used. For Iain, the commutation lemma 2.1 shows that 

IaUi<CK^ r \\fh\\^/^im\\i\\iu\\i 
Jo 

'[ WV^VtttWl + Cs [ ||/||3lNt||I 

Jo Jo 
Summing over I = 1, L, we integrate-by-parts in time and write the term Ia2 as 

Ia2 = - [ [ m„G{^,VO{fitt)t+ f fVvuG{^,WOfht . 
Jo Jq Jn Jo 

This is estimated in the same way as term Tan- The term is handled in the identical 
fashion with the same bound. Thus, wc have shown that 

I<5 sup E^{t) + Moid) + OTP { sup E^{t)) . 
te[o,T] te[o,T] 

Using (12.25b) for term //, integration by parts in time gives the identical bound as for 
term /. For term ///, a different approach is employed; we use (12.25a) and integration by 
parts in space rather than time, and let F^^ :— 3J^^(^*^^' — A'^^.Af) to find that 

in = - f f vUv:,sF:t<luU+ f I vl,F^tNkquv:,s 
Jo Jo. Jo Jr 

The Cauchy-Schwarz inequality together with the pressure estimate (12.1) give the bound 
CTP{supf.fz[o T] EK,{t)) for the first term on the right-hand side. The boundary integral term 
requires integration by parts in time: 



'0 JT JT 

= : Ilia + Ilh ■ 



Jo Jt Jt Jn Jo Jt 



First, note that 

F^fNk = 3 J-^A^ [{qAilNkh - 2^*5* TV, - g^^Af AT,] 

- 3J-^At [{qA^,Nk)tt - 2qtdtA^,Nk - qd^A^.Nk] . (12.26) 



46 



D. COUTAND AND S. SHKOLLER 



Next, substitute the boundary condition (12.8), written as 

-qA'UNk = ^/^Ag{fi^) [Sji + {{h^)j - fij)hi + nj{{hf,)i - fn)] + k(V9o5o^['5 • n],[j ),q fii , 

(12.27) 

into (12.26). The two bracketed terms in (12.26) are essentially the same, so it suffices to 
analyze just the first term. We begin by considering the term v^Ag(?7') in (12.27). 
Then can be written as 

(12.28) 

= 3 ^(^/^r^^^/3 MJ-'v:,s A^^U -2J-'qtdtA^Niv:,s vi - J-\d^A^Niv:,s vi 

<6 sup E^{t) + Mo{6) + CTP{ sup E^{t)) , 

te[o,T] te[o,T] 

the last inequality following from the fundamental theorem of calculus and the same argu- 
ment we have used above. 

In order to estimate 1 1 lb, because we do not have a trace estimate for dfA, we let 
Qr{df]n) := v^(^K)r and compute 

Since ^/g^{nn)rqu = {■\fg^{nK)rq)tt - 2dtQrqt - d^Qrq, it follows that 

r-T 



(12.29) 



Illb = -3 



VttW9'^in^)rQ)ttK,s At + vlt{Qrq)tt{vl,s )t 



il,{2Q:^ilo. il,s Alq, + Q'^i'^il^ A^q + Q^AiU y^^s Atq)t 



12.30) 



Using the pressure estimates and by definition of our energy function, for t G (0,T), 

te(i) e F0-5(r), Qr{t),dtQr{t)eL^{T), d^Qrit) & l\v), 
dvt{t) e H'-'{r), dvn{t) e l\t) ■ 

thus, all of the terms, except the first, on right-hand side of (12.30) can be easily bounded 
by CTP(supj£[Q E^{t)). Integrating by parts in space, the first term in (12.30) has the 
following estimate: 

3 r I ^^.r.^'ali^tt,0{ilA,sJ^'A:lU+V¥M.^ 

Jo Jr 

+ [P,f{dfj, dv)vU +Pndn, dv)]{v>ls J-'A^U < CTP{ sup E^t)) . 

te[o,T] 

(12.31) 

The remaining three terms in the boundary condition (12.27) are now considered. The 
additional integrals which arise in the (12.30) are given by 

- 3 / / d^{^/^Hh ■ (n„ -h)hr + \/^Hn ■ ((n«)r - fir) + k^/YkAo{v ■ nK)n„} v'^,s a-wj^ 
Jo Jr 

=: Ji+J2 + J3- 
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Term J3 with the artificial viscosity provides the integral k jQ[dfiv ■ Uk), P(yf]^ Vv) ■ d1v\i. 
The highest-order terms in this integral are estimated as 

Jo 

<C{5) f \\P{\/f,,\/i)\\L^\d^v\l + K6 [ {\V^d^v\l + \V^d^v,\l) 

Ja Jo 
<C{5)TP{ sup E^{t)) + S sup E^{t). 
te[o,T] te[o,T] 

The lower-order terms in J3 also have the same bound. As to terms Ji and J2, the es- 
timates for the terms with {^Hn)ttt are obtained exactly as in (12.31). For the terms 
that contain hut^ we use the formula (12.29) for the third time derivative of the unit nor- 
mal; it immediately follows that terms J\ and J2 are also bounded by ^suptg[o,T] ^^{t) + 

CTP{sup,^[o,T]E4t))- 

We need only consider now the additional terms in (12.28) from the remaining three 
terms in the boundary condition (12.27). The only novelty is in the highest-order integral 
coming from integration by parts in space in the k term: 



K [ Pf{T)vl,^ {T)vU{T) 
J ^ 



'1 

where i^f(0 8.ud dtP^^" {t) are both in L°° (F) for each t & [0, T] . Using the fundamental the- 
orem of calculus and the fact that ^/nvut S ^^(0, T; H^-^{Q.)) together with Jensen's inequal- 
ity shows that this term is bounded by Mq{5) + ^ supjgjg.T] -^kC^) + C"rP(supjg[Q 7^] E,^{t)). 

To study term IV, we use (12.25) together with the incompressibility condition (5*,^ Af)ut 
to find that 

IV^~ r [ [(3«j„fe dtA'l + v\k d!Al)qtu + d^A^^qtu] =: + m . 
Ja Jn 

For IVb, we integrate by parts in time: 



IV, = f f SivU d^A^)tqu- I M^kdU^t 
Jo Jn Ja 



Since dfA is bounded in H°-^{Q), the spacetime integral is easily bounded by CTP(sup(g[o,T] Enit)); 
meanwhile, the remaining space integral satisfies 



/H2 < 3 / [vlk iO)d^ A'^mqttiT) + I [ [vlk d^^A% dt qu{T) + Mo 
Jn JnJo 

<S\\quiT)\\l+Mo{5) + T sup \\{vU d^AlUo\\qu{T)\\o 

te[o,T] 

<5 sup E4t)+Mo{S) + CTP{ sup E^t)) . 

te[o,T] te[o,T] 

With := J,7i(A^ylf - A'^Af), IVa is written as 

IVa = - f I {Zd^^v^s F^^vU +dy^,s F:tv\k )qm + v\k P^{J-\ A, Vv.)dtvi,i 
Jo Jn 

=:IVa,+IVa2 + IV„s- (12.32) 
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Term IVa^ is estimated in the same way as term IVb- For term IVa2, we integrate by parts 
in space to find that 



Jr 



Jn 



The first integral IVa2i handled identically to term IIIj, to give the bound CTP(sup(£[Q 
We write the second integral as 

Jo Ju 

integrate by parts in time, and obtain 

IVa2u = - r [ mvlF^riV\k))tqtusHdK{KiV^ 
Jo Jn 

+ [ [{d^v:F^tv\k)qtt,sHd^v:{F^tv\k),sq^ 
Jn 

Since O^Vk is bounded in and (FVf)t is bounded in i°°(0), the spacetime inte- 

gral is boimded by CTP(supfg[Q E,-{f)). Next, we analyze the highest-order term in the 
remaining temporal boundary integral in IVa2'- 



j\ 

Jn 



[{dtvlF^^v\u)qu,s]{T) 

I id?v:{0)F^HO)v\kmqu,s{T)+ [ (d^vlF^tv^Mus (T) 
Jn Jn Jo 

<Mo{5) + 5\\quiT)\\i + T sup \\{d^i:F^ti\k)th\\qtt{T)h 

te[o,T] 

<Mo{5)+5 sup E4t) + CTP{ sup E^t)) . 
te[o,T] te[o,T] 

The remaining term is analyzed in the same fashion. 

Step 3. The inertia term. Finally, the inertia term in (12.7) satisfies 



Jn 



+ 



fl 

Jo Jn 



Since dtJ,, = Trace ( o„ Vi;„), d"^!,, = Trace (a^ V5t{j„) + P{Ji^'^ ,an,\/v^,), and 9|J«; = 
Trace(oK Vd^Vn) + P( J~^, a^, Vw„) S/dtv^, then using condition (3.3c), we see that 

sup \\vttt\\l<\\vtttm\l + CTP{ sup E,{t)). 

t£[0,T] te[0,T] 

From (12.2) evaluated at t = 0, we see that ||vttt(0)||o < Mq, so the lemma is proved. □ 

Lemma 12.3 (Energy estimates for the second time-differentiated K-problem). For Mq taken 
as in Lemma 10.1 and 5 > Q, solutions of the K-problem (4-1) satisfy: 

i-T 



sup \d2Vt-n\l + / • n^lo 

te[o,T] Jo 



<Mo + TP{ sup E^{t))+6 sup E^{t)+CP{\\V^vt\\h,oT-HSHn)))- (12.33) 
te[o,T] te[o,T] 
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Proof. We let ddf act on (4.1b) and test with QfOvn where = Ui, and ai is an element 
of our partition of unity. This localizes the analysis to a neighborhood of the boundary 
r where the tangential derivative is well-defined. In this neighborhood, we use a normal 
coordinate system spanned by {dirjo, d^rjo^ N). 

We follow the proof of Lemma 12.2 and replace with dd^. There are only two dif- 
ferences between the analysis of the second and third time-differentiated problems. The 
first difference can be found in the analogue of term III in (12.9), which now reads 

Jj, P{dri,dv) dvtd^vtt- After integration by parts in space, this term is bounded by 
CrP(supjg[Q Eii{t)): however, this term requires a bound on which requires us to 

study the third time-differentiated problem. (Compare this with the third time-differentiated 
problem wherein integration by parts in time forms an exact derivative which closes the es- 
timate.) 

The second difference is significant. Because the energy function places 
iui e L^{0,T;L^in)) and m G L°°{0,T; H'-%n)), 

there is a one-half derivative improvement that accounts for (12.33) being better than (12.6). 

In particular, the analogue of term // in (12.20) is /g /p \/gHn ■ 99j (n^ — h) ■ dvu, 
and since \ddf{hK. — n)\o < C |P(977)|loo |ut|2, then this integral is easily seen to be boimded 
by CTP(supj£[o,T] -^k(O)- (This is in sharp contrast to the difficult analysis required at 
the level of the third time-differentiated problem which follows equation (12.20).) 

All of the other estimates follow identically to the proof of Lemma 12.2 with being 
replaced by dd^. □ 

Lemma 12.4 (Energy estimates for the time-differentiated K-problem). For Mq taken as in 
Lemma 10.1 and ^ > 0, solutions of the K-problem (4.1) satisfy: 

sup \d^v-n\l-\- I \\/Kd3Vt ■ ni^\l 
te[o,T] Jo 

<Mo + rP(sup E^{t)) + 5 sup E^{t)+CP{\\^/llv\\l^,^T.H...<n)))- (12-34) 

te[o,T] ie[o,T] 

Proof. After replacing dd^ with d^dt, the proof is the same as the proof of Lemma 12.3. □ 

Lemma 12.5 (Energy estimates for the K-problem). For Mq taken as in Lemma 10.1 and 
6 > 0, solutions of the K-problem (4.1) satisfy: 

sup \d'^ri-ri\l+ [ \^/K^4V■n^\l< Mo + TP{ sup E^{t)) + 5 sup E^{t) . (12.35) 
te[o,T] Jo te[o,T] te[o,T] 

Proof. Let act on (4.1b) and test with d^v. All of the terms are estimated as in Lemma 
12.3, except the analogue of (12.15) which reads, after replacing df with d^, as 

K [v- d^n^,h^ ■ d^v)]o . 
Jo 

Since the energy function places ^/Kfj G L°°(0, T; ff^'^(f2)), we see that this integral is 
bounded by S sup^gjo^y] E^it) -|- CTP(sup(g[o,T] E^it)). □ 

To the above energy estimates, we add one elliptic estimate arising from the modified 
boundary condition (12.8). We will make use of the following identity: 

[V9\ivi]n = [V99''''^rf,0j +V9{9''''9"'' ' Q'^" g^^ri^p rf rf ,^.-y U ■ (12-36) 



50 



D. COUTAND AND S. SHKOLLER 



Lemma 12.6 (Elliptic estimate for ^fj). Let Mq be given as in Lemma 10.1. Then for 
5>Q, 

sup \^fi\l{t) < Mo{5) + 5 sup E^{t) + CTP{ sup E^)), (12.37) 

te[o,T] te[o,T] te[o,T] 

Proof. Letting Q"* := Q'^^^{df]) denote a smooth function of Sry, from (12.36), wc see that 

+ [dQf'd^^^p, +d'Qf'dfi\p, +d'Qf'v',0,U 

The estimate (12.37) is obtained by letting d'^d^ act on the modified boundary condition 
(12.5) and then testing this with (^"^g'^^Uld^fi^s where = ^i- 

For convenience we drop the subscript i from fli and Fj. (Recall that ai denotes the 
partition of unity introduced in Section 2.) For the surface tension term, integration by 
parts with respect to da yields 

r[d\^gHn o fi),^ , ]o < r [-|CV5'n95j?|2 + C\F\l^ \rjU \mv\o\ , 

Jo Jo 

where F := P(^, dr], dv, d^rf) is a polynomial of its arguments. To get this estimate we have 
used the fact that 

since t?*,/?!!; = 0. (This ensures that the error term is linear in |n5^77|o rather than qua- 
dratic.) 

We next analyze the artificial viscosity term. The testing procedure gives us the integral 

- / K[d^d^Aoiv ■ n^)h^,C''9^^'nd\5]o- 
Jo 

The positive term comes from d'^d-y acting on v. This gives, after integration by parts in 
space, the highest-order integrand K{d^v,a'y ■nK.)g'^^gQ^{nK-^d^fj,i3s ), where 11 = n(g)n. We 
can write this term as 

k(S3s,„^ •n«)5^^5o^(ri« • d%0s ) + ^(a^i;,^^ ■n,)g^^gf{h^ ■ (E - n,)d%0s ), 

where 11^ = nn^fiK- The first term is an exact derivative in time, and yields 

The space integral of the second term is estimated by C |F|i,oo |-0tw|5 |y^^|5|nre — ni^oo 
and |n« -niioo < C«;|j?|3.5. Prom (12.5) 

Kd^Ao{v ■ n«) = -d^{^/^~^^/§Ag{f]) ■ n«) + d^q. 

Thus, 

Cy/K\KV\5 < K;|F|i,oo(|{;|3 -I- \y/Kv\4:) + (k^ + )|i^|L~ I^U + | V^??|5 + V^l^ls , 

SO that 

/ / Kid^v,a,^ ■h^)g^^gfih, ■ (n - n^)d^fi,0s ) 
Jo Jr 

<C [ [V^\F\Loo{\v\3 + \y/Kv\4 + \fiU + \q\3)\V^fi\5+ \F\L^\V^fl\l] ■ 

Jo 
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Having finished the estimates for the terms leading to the positive energy contribution, we 
next consider the most difBcult of the error terms. This occurs when d^d-^ acts on n^, 
producing the integral jj" w • d'^n,^, ■ S^iyjo- 

To analyze this term, we let v = v'^d.yT)^ + VnfiK', we remind the reader that for 7 = 1,2, 
d^fjn spans the tangent space of r?K(i)(r) at fiK{x,t), so that d^fjn is orthogonal to n^- It 
follows that v'^d-yfji^ ■ (9^n„ is equal to —v'^fii^ ■ d^d^fj,^ plus lower order terms v"'TZj{fj), which 
have at most only five tangential derivative of 77^- Note also that since hK,-hn = 1, • d^n^ 
is also a sum of lower order terms which have at most only five tangential derivative of fjn- 

Thus, 

/ «;[C^ V ■ d^h^, ■ d^fi]o = / k[C^ v^Hk ■ d^d^fj^ + v''TZj{fi), n« • d^fi]o , 
Jo Jo 

where the remainder term satisfies 



[ [Kvmjifi) , ■ d'r]]o <C [ \F\Lo<.[\V^v\4\y/^v\5 + \V^v\ 
Jo Jo 



We must form an exact derivative from the remaining highest order term 



/ 

Jo 



(12.38) 



and this will require commuting the horizontal convolution operator, so that the fj on the 
right side of the (F) inner-product also has a convolution operator, and is hence converted 
to a • d^PL *h V term. With this accomplished, we will be able to pull-out the operator 
and form an exact derivative, which can be bounded by our energy function. 

Noting that on T the horizontal convolution -k^ restricts to the usual convolution * on 
M?, we have that 

K 

??K = ^ Vol pi * [pi * {y/cHV o 0i)] 

For notational convenience, we set p = Pi/k, C = \/^' ^^'^ ^ ~ ■'■l^ ~ ^^^O-' ^ ^i)- 
follows that (12.38) can be expressed as 

I I {v^n^) o Oi ■ d^d^ \C{ei)p * p * {Qf}) o Oi] (n« • 8^) o Oi (12.39) 

Jo Jr 

With :— d^p * Q{9i)fi{6i), we see that 

Kn^ • d^d^ [C,{9i)p * p*{C,f)o Oi)] = ku^ ■ {d^d^C{ei))p* p* (C?7(6'i) + nCfi^ ■p*d^g^ + KJl^if)) , 

(12.40) 

where the remainder TZ^if]) has at most five tangential derivatives on rj. Substitution of 
(12.40) into (12.39) yields three terms, corresponding to the three terms on the right-hand 
side of (12.40). For the first term, we see that 

TkV / id^d'aOi))i''n^- P*P*mei) {h,-d^fi{9i)) <c r \\V^eih.5mi)\L^\y/^vU. 

Jo Jr Jo 

The second term on the right-hand side of (12.40) gives the integral 



^0 Jr 



P * d^g^) (n« • g^) + TZeifi), 
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where the remainder TZeifj) is lower-order containing terms which have at most four tan- 
gential derivatives on fj and five on ({9i). 

We fix i e {1, K}, drop the explicit composition with 9i, we set 

^ii^Xd'^v =nK- P* Cd^V - P*{nK- (d^v) , 
and analyze the following integral: 

/ {C«« • P * 5«} {«K • d^n} = / {P * • 9k)} {n^ ■ Cd^fj} + / An^^g^ {n« • Cd^fj} , 

JR JR JR 

JR JR JR 

where the remainder TZrifj) comes from commuting with the cut-off function C, and has 
the same bound as TZ4{fi). The first term on the right-hand side is a perfect derivative and 
for the remaining terms we use Lemma 2.1 together with the estimate kI^kIo.a < C||??||5.5 
to find that 



K f {Cn. ■P*9.} {n« • d% < C \F\l^ \V^fj\l 

JR 

Thus, summing over i G {1, ...,K}, 

k[ [C^v-d^h^,n^-d^fi]o<C [ |i=^|Loo(|^/;ir|6 + |\/;^i;|4 + |\/^??|5)|\/)^^|£ 
Jo Jo 



where |\/Kr|6 := maxjg^i^...^^} |\/K^i|6- 
It is easy to see that 



Jo Jo 



c;|2 
5 



with the same bound for [d'^d^{^Ag{fj) ■ nK{hK — n), C'*5'''*^n9^77,5 ]. With (3.3), we infer 
that 

W^d'f, ■ h,\l{t) <Mo + c [ \q\l + c [ (I V^rie + iV^iU + \V^fi\5) IV^vU ■ 

Jo Jo 

Adding to this inequality the curl estimate (10.7) for the divergence estimate 

(which has the same bound as the curl estimate), and using Young's inequality to get 
\F\l^ \^/Kv\i \y/Kfi\5 < 5 \y/Kv\l + Cs Jq Wi^vll, we See that 

sup \y^v\l<Mo + CT sup (|kt||L + MIL + hills) 
te[o,T] te[o,T] 



+ CT sup {\F\l^{\^T\l + \^fi\l)+5 [ \^v\l, 

tG[0,Tl Jo 



-_[o:. 

from which the lemma follows. □ 

12.1. Removing the additional regularity assumptions on the initial data. At this 
stage, we explain how we can remove the extra regularity assumptions on the initial data, 

Uo and f^, so that the constant Mq depends on HuqIU.s and Irjs.s rather than \/k||'^o||io.5 
and v^lrlT as stated in Lemma 10.1. The modification requires the following regularization 
of the initial data: set uq = p^-^ * Eq^{u{0)), where fl^ is obtained by smoothing Cl via 
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convolution with Pe-'^j i-e., we use p^-^ * 9i as our family of charts. We make use of the fact 
that 

P{\\V^m\\io.5) < cp||uo||4.5) , P(lv^riio) < CP(|r|5.5) , 

which follows by integration by parts of six tangential derivatives onto the moUifier Pe-'^i 
this and results in the constant C > being independent of k. 

12.2. The limit as k ^ 0. 

Proposition 12.1. With Mo = P(||wo||4.5, jrjs.s) a polynomial of its arguments and for 
Mo > Mo, 

sup E,{t) < Mo , (12.41) 

t£[0,T] 

where T depends on the data, but not on k. 

Proof Summing the inequalities (12.6), (12.33), (12.34), (12.35), and (12.37), and using 
Lemma 10.1 and Proposition 5.1, wc find that 

sup E^{t)<Mo + CTP{ sup E,{t))+6 sup E^{t) , 
te[o,T] te[o,T] te[o,T] 

where the polynomial P and the constant Mo do not depend on k. Choose 5 < 1. Then, 
from the continuity of the left-hand side in T, we may choose T sufficiently small and 
independent of k, to ensure that (12.41) holds. (See [8] for a detailed account of such 
polynomial inequalities.) □ 

Proposition 12.1 provides the weak convergence as k of subsequences of (w, q) toward 
a limit which we denote by {v, q) in the same space. We then set rj = ld + J^v, and 
u = V oT]~^. It is obvious that v^, arising from the double horizontal convolution by layers 
of w, satisfies i^^vin L^{0, T; H^-^{n)), and therefore ^ r; in ^^(o, T; ^^(f]). It follows 
that divu = in r]{Cl) in the limit as k ^ in (1.3c). Thus, the limit {v,q) is a solution 
to the problem (1.3), and satisfies Eo{t) < Mo. We then take T even small, if necessary, to 
ensure that (3.3) holds, which follows from the fundamental theorem of calculus. 

13. A POSTERIORI ELLIPTIC ESTIMATES 

Solutions of the Euler equations gain regularity with respect to the EK,{t) from elliptic 
estimates of the boundary condition (1.3d), which we write as s/gHn{ri) = s/gqn. 
Replacing d-^ with dt in (12.36), we have the identities 

dti^Hnorjy = -[^/55"''n>^^+V5(rt"^ - 5"V'^y,/3»?^.^^^],a • (13.1) 

and 

(13.2) 

where Qf^ = Q{dr]) is a rational function of drj. 

Lemma 13.1. Taking Mo a,s in Proposition 12.1, and letting M.o denote a polynomial func- 
tion of Mo, for T taken sufficiently small, 

sup [|r(t)|5.5 + \Ht)\\i.5 + \\Vt{t)h] < Mo . 

te[o,T] 
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Proof. We begin with the estimate for vt- FoUowing the proof Lemma 12.6, wo lot d-ydf 
act on the boundary condition (1.3d) and test with —(^g^^Ill.v^,s, where = a,, an 
element of our partition of unity. Using (13.2), we sec that — Jp d^/df[^Hn{ri)] -(^^g^^vtjS = 
J-p{^qn)tt ■ [C^g^^vt,s],-y Using (12.41), letting C denote a constant that depends on Mq, 
and summing over the partition of unity, we find that 

\d\t • n|o < C [\v\2 + \v\2 {\vt\i + \v\i) + \{q^n)tt\o] \d\ ■ n|o + C\vt\l\v\3 ■ (13.3) 
This follow since 

"V5(5''^5"''-5"V^)'?%/3V,.^^^, 

while 

[{g^'lll),„v^,s] < C \vt\i {\d\t ■ n|o + kt|i Ivls) • 

Applying Young's inequality to (13.3) yields, after adjusting the constant, 

\d^Vfn\l<CMl + \v\l + \vt\l + \qu\l]. 
A similar computation shows that 

\d\ ■n\l<C [\r]\l + \v\l + \vt\l + \qtt\l] ■ 

Thus, by interpolation 

sup \d^vt ■ n\lj^ < C sup [\r]\lj^ + \v\lj^ + \vt\l,5 + \qtt\l.5\< -Mo- 

te[0,T] t(£[0,T] 

Computing the iI^(0)-norm of (10.9), we find that 

sup \\ cmlvt\\2 < Mo + C T sup \\vt\\3, 
te[o,T] te[o,T] 

with the same estimate for supjgjg || divui||2. Hence, for T taken sufficiently small, we 
infer from Proposition 5.1 that 

sup ||^;t||3<^o. (13.4) 
te[o,T] 

Next, we let d^d^dt act on the boundary condition (1.3d) and test with —C,'^g'^^Yl\d'^v'^,s. 
Using (12.36), we find that 

sup \d\-n\l<C sup Ml + \v\l + \(lt\l]<Mo. 

*6[0,T] te[O.T] 

Computing the Jf"^-^(ri)-norm of (10.8), and again taking T sufficiently small, we see that 

suPte[o,T] II^^IU.s < Mq. 

In order to prove our remaining estimate, we need a convenient reparameterization of 
r(t) via a height function h in the normal bundle over F. 

Consider the isometric immersion 770 : (F, .go) (K'',ld). Let 6 = F x (— e, e) where e 
is chosen sufficiently small so that the map _B : 6 ^ K"* : (y, z) y + zN{y) is itself an 
immersion, defining a tubular neighborhood of 7?o(F) in M?. We can choose a coordinate 
system a = 1,2 and ^. 

Let G = i3*(Id), denote the induced metric on and note that G{y, z) = Gz{y) + dz®dz, 
where Gz is the metric on the surface F x {z}, and that Gq = go- Let h : T ^ (— e, e) be 
a smooth function and consider the graph of /i in B, parameterized by ^ : F — > B : ?/ 1— > 



[g-^'li^^vl,^] =0, 



1/3 V )i/ 
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{y, h{y)). The tangent space to graph(h), considered as a submanifold of B, is spanned at a 
point (l){x) by the vectors 

9 d(f) d dh d 
^*^dy^' ^ dy^ ^ dy^ ^ dy^dz' 
and the normal to graph(h) is given by 

nfe) = ^'(»)(-G;f.,|i^ + |) (13.5) 

where J/j = (1 + h^aG'^^y-^h^f))^/'^ . Therefore, twice the mean curvature H is defined to be 
the trace of Vn while 

(Vn).,=G(V!^n,^) 
where — for a = 1, 2 and = Substituting the formula (13.5) for n, we see 



that 



—) 

dy^ ' dzV dyl^J 
{GMJ-'Gth^^),o. + Flp{y,h,dh); 



(Vn)33 = Gfv_ 
V a 



'dz) 



= Flp{y,h,dh) 

for some functions Fl,^^ and -Ff^-i, a.0 = 1- 2. Letting 70 denote the Christoffel symbols 
associated to the metric go on T, wc find that the curvature of graph(/i) is given by 

L,ih) ■.^H = -{J-'GthrXs + JuWl,]], - Gth,,[^o\]s) ■ (13-6) 

Note that the metric Gh = P{h), and that the highest-order term is in divergence form, 
while the lower-order term is a polynomial in dh. The function h determines the height, 
and hence shape, of the surface r(t) above F. 

Given a signed height function /i : Fq x [0, T) M, for each t G [0, T), define the normal 
map 

r/^ : To X [0, T) ^ T{t), {y, t) ^ y + h{y, t)N{y) . 

Then, there exists a unique tangential map t?'^ : Fq x [0, T) — » Fq (a diffeomorphism as long 
as h remains a graph) such that ?7|r(i) has the decomposition 

r?|r(-, t) = f^^i; t) o r^^{-,t), r7|r(y, t) = v^{y, t) + hW{y, t), t)N{j^^{y, t)) . 

The boundary condition (1.1c) can be written as aLh{h) = qo {r]'^)~^ . 

The operator Lh is a quasilinear elliptic operator; from the standard regularity theory 
for quasilinear elliptic operators with coefiScients on a compact manifold, we have the 
elliptic estimate 

\h\5.5<C\qo{rj-)-^\s,^<C\\q\\i. 
By (12.2), we see that for all t e [0, T], 

h\\4<C\\a\\l\\v\\l + C\v\3\vt\2<Mo, 

the last inequality following from (13.4). 

Since T{t) =graph/i(t), this estimate shows that T{t) maintains its if^-^-class regularity 
for t<T. □ 
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14. K-INDEPENDENT ESTIMATES FOR THE SMOOTHED PROBLEM AND EXISTENCE OF 

SOLUTIONS IN 3D 

The 3D analysis of the K-problem requires assuming that the initial data uq e H^'^{0) 
and r of class H^-^. This is necessitated by Sobolev embedding H'OfUL^ < C'll'^tHa- 

By replacing the third-time differentiated problem with the fourth time-differentiated 
problem the identical analysis as in Section 12 yields 

£;^^(i)<Mo, 

where Mq is a polynomial of ||uo|l5.5 and Irjg.s. (In fact, our analysis in Section 12 used all 
of the 3D terms and notation, so no changes are required other than raising the regularity 
by one derivative.) 

We let [v^q] again denote the limit of {v^q) as k ^ 0. The identical limit process as in 
2D, shows that {v, q) is a solution of the Euler equations. 

Having a solution {v, q) to the Euler equation, we can use the a posteriori estimates (13.1) 
as a priori estimates for solutions of the Euler equations. We see that 

sup [El'^it) + |r(i)|5.5 + |k(t)||4.5 + \\vt{t)h] < Mo , 

te[o,T] 

where Mo is polynomial function of ||wo||4.5 and |r|5.5. This key point here, is that the 
elliptic estimate for vt E H^{n) improves the regularity given by E'^^{t) and allows for the 
required Sobolev embedding theorem to hold. 

Since our initial data is a priori assumed regularized as in Subsection 12.1, we see that 
solutions of the Euler equations in 3D only depend on Mo- 

15. Uniqueness of solutions to (1.3) 

Suppose that {rj^,v^,q^) and {rj^^v^^q^) are both solutions of (1.3) with initial data 
uo e H^-^{9. and T e H^-^. Setting 

e,{t) = Y.\\dHt)\\l,-k, 

k=0 

by the method of Section 12 with k = 0, we infer that both (t) and £"^2 (t) are bounded 
by a constant Mo depending on the data uq and F on a time interval < t < T for T small 
enough. 
Let 

w := — , r := q^ — q^, and ^ ■= il^ — ij'^ ■ 
Then (^, w, r) satisfies 





J ^ 


in r2 X (0, T] , 


(15.1a) 






in 17 X (0, T] , 


(15.1b) 


(o^)>'u = 




in Q X (0, T] , 


(15.1c) 


rni = 




on r X (0, T] , 


(15. Id) 




(0,0) 


on f2 X {t = 0} . 


(15. le) 



Set ^ 

fe=0 
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We will show that E{t) = 0, which shows that w — 0. To do so, we analyze the forcing terms 
on the right-hand side of (15.1b) and (15.1c), as well as the term aAgi_g2{r]'^) in (15. Id). 
Wc begin with the third timc-difFcrcntiated problem, and study the integral df[{a^ — 

a^) Vg^] Wttt- The highest-order term is 

^ ^(fli - a^)VqltWttt < \\a^ - a'Hioo + ^ ^ \\wut\\l < CtP{E{t)) . 

The third space differentiated and mixed-derivative problems have forcing terms that can 
be similarly bounded. 

The difference in pressure r satisfies, using the notation of (12.2), satisfies the following 
Neumann problem: 

LaAr) = -dta^lw\j +[a^l{a'^ - a^)iq^,k],j infl, 

Ba, (r) = -wt ■ y7ni + a^lia^ - a^t<f.k N, on T . 

Since ^'dl^y^lU.s) is bounded by some constant C = P(Alo)7 (12-3) provides the estimate 

lk||3.5 < C[\\al\W^ \\wh.^ + lla^lb.s lk'||3.5 Ik' - a'l|2.5 + Ix/^n'b Iktlb.s] • 

Since \\a} — a^||2.5 < C'llClls-s, and ||aj||i.5, ||a^||2.5, lk^||3.5, and \\/g^'n}\2 are aU bounded 
A^O) we see that ||r(t)||3.5 < CP{E{t)). Similar estimates for the time derivatives of r show 
that ||r(t)||3.5 + ||rt(t)||2.5 + ||r«(t)||i < CP{E{t)). 

This shows that the energy estimates of Section 12 go through unchanged for equation 
(15.1); therefore, using (15. le), we see that sup(g[o,T] -^(0 ^ CTP(sup(g[o,T] ^(i))- 

16. The zero surface tension case a = 

In this, the second part of the paper, we use our methodology to prove well-posedness 
of the free-surface Euler equations with a = and the Taylor sign condition (1.2) imposed, 
previously established by Lindblad in [13]. The main advantages of our method over the 
Nash-Moser approach of [13] are the significantly shorter proof and the fact that we provide 
directly the optimal space in which the problem is set, instead of having to separately 
perform an optimal energy study once a solution is known as in [6]. If one uses a Nash- 
Moser approach without performing the analysis of [6] , then one obtains results with much 
higher regularity requirements than necessary, as for instance in [11] for the irrotational 
water-wave problem without surface tension. We also obtain lower regularity results than 
those given by the functional framework of [6] for the 3D case. 

We will extensively make use of the horizontal convolution by layers defined in Section 
3, and just as in the first part of the paper, for v G L'^{fl) and k G (0, kq), we define the 
smoothed velocity v'^ by 

K L 

v'^ = ^V^i[p^*h[pi.*h{{Va'iv)oei)]]oe:[^ + ^ aiv. 

i=l i=K+\ 

The horizontal convolution by layers is of crucial importance for defining an approximate 
problem whose asymptotic behavior will be compatible with the formal energy laws for 
smooth solutions of the original (unsmoothed) problem (1.1), since the regularity of the 
moving domain will appear as a surface integral term. In this second part of the paper, the 
properties of these horizontal convolutions will be featured in a more extensive way than in 
the surface tension case of the first part of the paper. 
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We remind the reader that this type of smoothing satisfies the usual properties of the 
standard convolution; in particular, independently of k, we have the existence of C > such 
that for any v G H''{Q): 

Mis < C and i+p < Ck~p\v\,_i for p > 0. 

We will denote for any I G {1, K}, the following transformed functions from v and 77 that 
will naturally arise at the variational level: 

Definition 16.1. 

= Pi *h iVrnv o 61) in (0, 1)^, 

T^"^ = Id+ [ m il, 







Jo 



= 0i+ vu in (0, 1)^ 



Remark 12. The regularity of the moving free surface will he provided by control of each 
riiK in a suitable norm independently of the parameter k. 

17. The smoothed k-problem and its linear fixed-point formulation 

As it turns out, the smoothed problem associated to the zero surface tension Euler equa- 
tions can be found quite simply and naturally, and involves only transport-type arguments 
in an Eulerian framework. Also, the construction of a solution is easier if we assume more 
regularity on the domain and initial velocity than in Theorem 1.4. We shall therefore as- 
sume until Section 26 that Q is of class and uq £ H^{fl). In Section 26, we will show 
how this restriction can be removed. 

Lotting u = V o 7]'^^^ . wc consider the following sequence of approximate problems in 
which the transport velocity u'^ is smoothed: 

Ut + Vu-u + Vp = in r]''{t,n), (17.1a) 

divM = in ?7'^(t,0), (17.1b) 

p = onri'^{t,r), (17.1c) 

u(0) = uo in fl. (17.1d) 

In order to solve this smoothed problem, we will use a linear problem whose fixed point 
will provide the desired solution. If wc denote by v an arbitrary element of Ct defined in 
Section 18, and f]'^ the corresponding Lagrangian flow defined above, then we search for w 
such that if u = w o {f]'^)~^ and u'^ = v'^ o (r?'*)~^, we have that 

ut + V u-u + Vp = in f]''{t,fl), (17.2a) 

diYU = Oinf]'^{t,n), (17.2b) 

p = on f;'=(t,r), (17.2c) 

m(0) = Uo in ^. (17.2d) 

A fixed point w = v to this problem then provides a solution to (17.1). In the following 
sections, v G Ct is assumed given, and k is in (0, kq). k is fixed until Section 20 where we 
study the asymptotic behavior of the problem (17.1) as k 0. 

Remark 13. Note that, for this problem, we do not add any parabolic artificial viscosity, in 
order to keep the transport-type structure of the Euler equations and to preserve the condition 
p = on the free boundary. 
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18. Existence of a solution to (17.1) 

18.1. A Closed convex set. 

Definition 18.1. For T > 0, we define the following closed convex set of the Hilbert space 
L\0,T;Hi{n)): 

CT = {vGL\0,T;Hi{n))\ sup 7 < 2||wolU + 1}, 

[0,T] " " 

It is clear that Ct is non-empty (since it contains the constant in time function uq), and 
is a convex, bounded and closed set of the separable Hilbert space i^(0, T; iJ 2 {Q)). 

By choosing TdlVuoHi +1) < Cneo, condition (3.1) holds for rj = Id + J^v and any 
V G Ct and thus (3.2) is well-defined. 

We then see that, by taking T smaller if necessary, we have the existence of ki > such 
that for any k G (0,Ki), we have the injectivity of r]'^{t) on Cl for any t G [0, T], and Vr/'^ 
satisfies condition (3.1). We then denote a'^ = [Vri'^]~^, and we let n'^{rj'^{x)) denote the 
exterior unit normal to ■r]'^{^l) at r]'^{x), with x € T. We now set K2 = min(Ko,«;i), and 
assume in the following that KG (0, K2)- 

18.2. Existence cind uniqueness for the smoothed problems (17.2) and (17.1). 

Suppose that v G Ct is given. Now, for v G Ct given, we define p on f]'^{t, SI) by 

Ap = -Ui,j Uj,i in fj'^{t,Q.), (18.1a) 
p = on ^'"(t,r), (18.1b) 

where u = v o (r?")"^. We next define in i7 by 

v{t) =uo+ f [Vp]{tf,r{tf,-))dtf, (18.2) 
Jo 

and we now explain why the mapping v ^ v has a fixed point in Ct for T > small enough. 
For each t G [0,T], let "^{t) denote the solution of A*(i) = in with 5'(t) = fj'^it) on 
r. For K and T taken sufficiently small \fj'^ — Id|4 << 1 so that ^'(t) is an embedding, and 
satisfies 

||*(i)ll4.5<q?7"(i)|4- (18.3) 
Letting Q{x,t) =p{^{x,t),t) and A{x,t) = [V*(a;,t)]~^, (18.1) can be written as 

[A'^4l,k],j = -[u-,ju,„]{^{x,t),t) in O, 
Q = on r . 

By elliptic regularity (with Sobolev class regularity on the coefficients [9]) almost everywhere 
in (0,r) and using (18.3), 

Mln-M < CPim,)\\v\\i\\v\\zPimi), (18.4) 

where P denotes a generic polynomial. Now, with the definition of T and Ct, along with 
the properties of the convolution that allow us to state that 

ru < -iris, 

K 

(since the derivatives involved are along the boundary, allowing our convolution by layers 
to smooth in these directions), this provides the following estimate: 
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where we have used the definition of Ct, and where denotes a generic constant depending 
on K. Consequently, we get in [0,T], 




< llwolli +C« / (18.5) 

With the definition of Ct, this yields: 

sup 11^(^)11 7 < \\uoh +C^T. 

Now, for G (0, T) such that T^C^ < 1, we see that v G Ct^, which ensures that the closed 
convex set Ct^ is stable under the mapping v ^ v. We could also show that this mapping is 
also sequentially weakly continuous in L^(0, T„; _ff 2 (fj)). Therefore, by the Tyclionoff fixed 
point theorem, there exists a fixed point v = v in. Ct^ ■ Now, to see the uniqueness of this 
fixed point, we see that if another fixed point v existed, we would have by the linearity of 
the mapping u — > p and the estimates (18.4) and (18.5) an inequality of the type: 

\\{v-v){t)\\i<C J\\v-v\\^, 

which establishes the uniqueness of the fixed point. By construction, if we denote u = 
V o (r?'*)~^, this fixed point satisfies the equation on (0, T^): 

ut + UiU,i +Vp = in f]'^{t, CI). 

Besides from the definition of p in (18.1), we have 

div Ut + Ui div u,i = in f]'^{t, fl), 

i.e. 

div u{t,f]'^{t,x)) — divMo(.T) = in f7. 

This precisely shows that u = v o (r^")"^ is the unique solution of the linear system (17.2) 
on (0,T,). 

Now, wc sec that wc again have a mapping v v from Ct^ into itself, which is also 
sequentially weakly lower semi-continuous. It therefore also has a fixed point in Ct^, 
which is a solution of (17.1). 

In the following we study the limit as k — !■ of the time of existence and of v^. We 
will also denote for the sake of conciseness Vk, Uk = VnOr}'^~^ and {u^Y respectively by v, 
u and u'^. 

19. Conventions about constants, the time of existence T^, and the 

dimension of the space 

Prom now on, until Section 26, we shall stay in M? for the sake of notational convenience. 

In Section 26, we shall explain the differences for the three dimensional case. In the re- 
mainder of the paper, we will denote any constant depending on ||mo|| a ^-s N(uq). So, for 
instance, with Qq solution of 

Ago = -uo,] uo,l in O, (19.1a) 
go = on r, (19.1b) 

we have by elliptic regularity ||golU — -^(^0) (since Q, is assumed in until Section 26). 
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We will also denote generic constants by the letter C. Moreover, we will denote 

K 

ii^^iu = Eii^^ii«.(o,i)- 

i=l 

Furthermore, the time > will be chosen small enough so that on [0,Tk], we have for 
our solution v given by Section 18: 

1 3 

- < detV??" < - in n, (19.2a) 

ll^l|3<NI + l, (19.2b) 
Ikl3< IkolU + l (19.2c) 
ll^lll < lko||| + l, (19.2d) 

The right-hand sides appearing in the three last inequalities shall be denoted by a generic 
constant C in the estimates that wc will perform. In what follows, we will prove that this 
can be achieved on a time independent of k. 



20. A CONTINUOUS IN TIME SPACE ENERGY APPROPRIATE FOR THE ASYMPTOTIC 

PROCESS 

Definition 20.1. We choose < ^ eC°°{Tl.) such that SuppS, C nf^j^^-i^lSuppai]" and£, = 1 
in a neighborhood ofT. We then pick < (3 £ I?(0) such that (3 = \ on [Supp^Y- We then 
define on [0, T^].- 

K 

E{t) =sup[E \\V^i{Oi)f,uh (o^^y + ll^'^lli + + ||t;||3 + 

[0,t] 

K 

+ supE[K;||v^i5oei||| (o ;^)2 +«;5||v^5o^;||4 (o_i)2] +1. (20.1) 

Remark 14. Note the presence of k- dependent coefficients in E{t), that indeed arise as a 
necessity for our asymptotic study. The corresponding terms, without the k, would of course 
not he asymptotically controlled. 

Remark 15. The 1 is added to the norm to ensure that E > 1, which mil sometimes he 
convenient, whereas not necessary. 

Now, since from Section 18, v € C'^(0, T„; Hi (fi)) (in a way not controlled asymptotically, 
which does not matter for our purpose), we have 77 e C°([0,T«];ff2 (rj)). Next with the 
definition (18.1), and the definition of our fixed point v, we have for p = g o (ry")"^: 

^P=-ut,3Uj„ in Tj''{t,n), 
p = Oon 77'=(t,r), 

which shows that q S C°([0, Tf^]; ifs (fj)). Consequently, is a continuous function on 

[o,r,]. 

Wc will then prove that this continuous in time energy is controlled by the same type 
of polynomial law as (41) of [8], which will provide a control independent of k on a time 
independent of k. 
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21. A COMMUTATION-TYPE LEMMA. 

We will need the following Lemma in order to later on identify exact in time energy laws 
from terms arising from our convolution by horizontal layers: 

Lemma 21.1. Let So > be given. Independently of k G {0,5q), there exists C > such that 
for any g e ((0, 1)^) and f £ hI ((0, l)^) such that 

So < mm{dist{supp fg, {1} x [0, 1]), dist{supp fg, {0} x [0, 1])), 

we have, 

*/i [/ff] - /pi ^/iffll 1 (0 1)2 <C ||K5l|i,(0,l)2|l/ll|,(0,l)^ 115110,(0,1)211/111, (0,1)=- 

Proof. Let A = p i [fg] — fpi-kh 9- Then, we have: 

A(x) = / pi.{xi-y\)[f{y\,X2)- f{xi,X2)] g{yi,X2) dyx, 

J X\—K 

this integral being well-defined because of our condition on the support of fg. We then 
have, since ]s embedded in L°° in 2d, 

PXI+K 

\A{x)\ < Ck||/|||_(o_i)2 / pi{xi-yi) \g{yi,X2)\ dyi, 

showing that 

||A||o,(o,i)2 < Ck||/||5 (o_i)2||pi l5lllo,(o,i)2 

<C«||/|||,(o,i).||5l|o,(o,i)- (21.1) 

Now, let p e {1,2}. Then, 

A,p = PL -kh [fg,p ]- f PL -*h 9,p +px *h [f,p g] - f,p pi- */i g- 

The difference between the two first terms of the right-hand side of this identity can be 
treated in a similar fashion as (21.1), leading us to: 

||0,(0,1)2 < C'^ll/llf ,(0,1)2 115111,(0,1)2 + *h [/,p5]l|0,(0,l)2 + ||/,p Pi -kh 5110,(0,1)2 

< C'k||/|||,(o,i)21I.91Ii. (0.1)2 + l|/,p,9llo,(o,i)2 + ||/,p IU~((o,i)2)||pi */t 5110,(0,1)2 

< C'kII/II 5 (Q i)2||.9||i^(o,i)2 + 2\\f,p ||l==((o,i)2)||5||o,(o,i)2 

< ^^«II/Ilf,(0,l)2|l5l|l,(0,l)2 + C||/||5,(o,i)2||g||o,(o,i)2. (21.2) 

Consequently, we obtain by interpolation from (21.1) and (21.2): 

l|A||i,(o,i)2 < C«ll/ll|,(o,i)2||5l|i,(o,i)2 +Ck^II/II|,(o,i)2||5||o,(o,i)2. 

□ 

We then infer the following result, whose proof follows the same patterns as the previous 
one: 

Lemma 21.2. Let So > be given. Independently of k £ (0, So), there exists C > such that 
for any g e H^{{0, if) (s = |, |; and for any f e Hi{{0, 1)^) such that 

So < mm{dist{supp fg, {1} x [0, 1]), dist{supp fg, {0} x [0, 1])), 

we have 

*h [fg]- fpx 9118,(0,1)2 < ||k5|U,(o,i)2||/|||,(o,i)2 +C'«^||5lL-i,(o,i)2||/ll|, (0,1)2- 
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22. Asymptotic regularity of the divergence and curl of fii,^. 

In this Section, we state the necessary a priori controls that we have on the divergence 
and curl of various transformations of v and fj. This process has to be justified again, since 
the functional framework substantially differs from the case with surface tension, with in 
this case one time derivative on the velocity corresponding to half a space derivative. 

We will base our argument on the fact that the divergence and curl of u satisfy the 
following transport type equations: 

Adivu = 0, (22.1a) 
Acurlu + ut,i -u^,2 u\i = 0. (22.1b) 

We now study the consequences of these relations on the divergence and curl of fj in the 
interior of ft, and of each 77;^, {1 < I < N). 

22.1. Estimate for div(^^),s. From (22.1a), we then infer in CI that: (a'')^{j,* = 0. Thus, 
for s = 1,2 

and by integration in time, 

{anim^ij it) = (o) + [\-P{ani,s +{{ay,m,ij -fi(r)iKj )] 

Jo 

+ [\ayitm,ij. (22.2) 

Jo 

Consequently, 

dwm,s it) =[-(a«)i + 6i]m;ij it) + {a'^)m),ij (o) + /* {a%m,ij 

Jo 

Jo 

Jo Jo 

+ f\-p{ani,s v,i -p(r)i<j )], (22.3) 

Jo 

showing that 

||div(/37?),4f)|l3 < Ct sup[ |la«||3 Wpf^Wr] +C + Ct snp[ < CtEit)+C, 

[0,t] [0,t] 

(22.4) 

where we have used our convention stated in Section 19. 

22.2. Estimate for div[riiK,s °^r^ ° ('?'') "^l* Since deWOi — 1, we then infer in (0, 1)^, with 
bf = [Vfj'^oei]-^ that 

{bf)iivoei),i = o. 
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Therefore, as for (22.2), 

{bmiV^iv) o ei),i^ it) = {imv^iv o (0) 

Jo Jo 

+ fmmv^ii o 00,^, -v^imimv « oo^,]. (22.5) 

Jo 

Consequently, 

+ f Pi M [(hmv^ii ° (^iVsi -v^imimv o ^o,:,- ] 

Jo 

'Vi M [{bf)UV^iflo9i),l^]+R, 
with 3 < Ct£'(t) + C. Next, thanks to Lemma 21.2, 

t/ 

t/ 

< CK^E{t) + CtE{tf + C. (22.6) 
By successively integrating by parts in time and using Lemma 21.2, 

\\f^ PL [mUV^in o ei),% it)] - ^' mlp^ [(^/^^ o ei),% ] II 3,(0,,). 

< 11^* Pi [{b^)i{^ivoei),i^{t)\- j\ht)\pL *h [(V^* 0^04] II 1,(0,1)= 

+ \\[pi [{bmV^lfl°(^lVso{t)\-(ht){pL [(%/^^° ^04- ]]oll 1,(0,1) = 

< Ck5 E{t) + CtE{tf + C. (22.7) 
Consequently, with (22.5), (22.7) and (22.6), we infer 

l|div[r);.„ o [rr'w o m) - mupi *h [iv^iv ° ^z),:,- ill 3,(0,1)= 

< CK^E{t) + CtE{tf + C, 

showing that 

||divfc,«o0-i o (^'^)-'] II 3, -.(,,((0,1).)) < CK^E{t) + CtE{tf + C. (22.8) 
We now study the curl of the same vector fields as in the two previous subsections. 
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22.3. Estimate for curl(/3?7),s. From (22.1b), we obtain: 

(a'^)it;,2 -{ay^v,] = curlti(O) + /*[ {dy,v,l {a'^)^ + ^;^} {ay^v,l [a^^]. 

Jo 

Tlierefore, for s = 1,2, 

{a'^)i{/3v),% -{a'^)i{Pv),lj =/?curln(0)„ -/3[(a«)i„ {v),'^ -{a%,, {v),] ] 

Jo 

wliicii implies by integration in time, 

Jo 

-P [\i~ay,s{i),^^-i~ay„,{v),]]+ [\f + g] 
Jo Jo 



Jo Jo 
rt' 



with/(iO= / p[-v'^,){dyv,l{d-)^U andg{t') = [ (a^^^^,^ (a«)f]„. Now, 

Jo Jo 
since is a, Banach algebra in 2d, 



+ .' 

'o 

<N{u^) + CtE{t)+ j \\f + g\\.. (22.9) 



We now notice that 



m = - [ piv'^ZjidyeA^n'^+v^ii^yi^'ski^n^]- I Pvy^vi[{dy{dy,\, 

Jo Jo 

= f I3[ [i~ayi,l {d^^ + K,;- Cayia^^] - f iivy, v,l [(a'^)i(a-): 
Jo Jo 

+ [Pi iayiv,i {d-)^+f,,i,v^,) {dn{{d^)^]]i, 

which allows us to infer that 

ll/(*)lli < [\mi + \m\i + ll^lli)[l|a?||3||a'^||3||t;||| + ||a«||3||a«||3||i;,||5] 

J 

+ r ii^iiiiiaiii + iwrwi + imii + Mim\\a^miMi+N{uo) 

2 2 2 2 2 2 2 

< CtE{tf + N{uo) + CE{t). 

Since g{t} can be estimated in a similar fashion, (22.9) provides us with 

W^m'sj -{dym,ij II 3 < ctEitf + N{uo). 
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showing that, 

\\cml{pf,),, II 3 < II [ - [ [\at)i] m,% II 3 + CtE{tf + iV(uo) 

<CtE{tf +N{uq). (22.10) 

22.4. Estimate for curl[?7;K,s oO'^^ o (t?*^)"^]. In a similar fashion as we obtained (22.8), we 
also have here 

||curl[77;„, o^-i o (77'^)-!] II 3 ^.^(j,) < CtE{tf + iV(wo). (22.11) 

22.5. Estimate for KA\Y[{^/aiVt o6i),s o9f^ o (jj'*)"^]. By time differentiating (22.5) twice 

in time, wc find: 

+ mav^iv o -v^i{9iM)i{v o ^04- (22.12) 

Therefore, 

'^m)ii{^Mvt)oOi),ij (t)|||,(o,i)^ < CnEitf 

1,(0,1)^ 

+ ii'j[(^r)i',<*]tiii,(o,i)=ii(v^^^o^^o,j III, (0,1)^ 

< CE{tr + \\4ibf)lMiAoA)4iV^ii'°()i)^]\\i(o,ir- 

(22.13) 

Next, we for instance have: 

/Ik \1 V2o9i,2s Ij 4.r-K /dm ?72 ° ^^2 



detV(77« o 6*0 ' ' det^V(r)'^ o 610 

K;[det(77 o 6';)]* — 5 — 7: 7 ~ K;[det(r; o 6';^' 



' detV(^« oOi) ^ ^ detV(^« o9i)' 

which shows that, 

Mb?,s)t\\l(o,iy<CE{t), 

and with (22.13) this implies: 

'^\mi{{v^ivt)o9i),ij (t)ii|,(o,i)2 < cE{tf. 

and thus, still by writing f/'^it) = r?'*(0) + / t;'*, we finally have 

Jo 

K||div[(V^5to^^,),s°^?^'° (^r)"']|||,,-,.(o) < CE{t)^. (22.14) 
With the same type of arguments, wc also have the following asymptotic estimates: 

22.6. Estimate for k cuT\[{^aivt o 61), s o9^^ o (yy**)"^]. 

4cur\[{^iit o 9i),s o^r' ° (r )"'] II ^ ^^^W'' (22-15) 

22.7. Estimate for div[(^{St o 9i),s 06^^ o {fj")-'^]. 

K^\\div[{^iVto9i),so9i' o (r)"']||5,^.(o) < C-^)'- (22.16) 
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22.8. Estimate for curl[(ya7'5t o 9i),s o0-'^ o (77'')"^]. 

K^\\cviTl[{^iVtoei),,o0-^o{,^^)-^]\\,^^,^^^^^<CE{t)^. (22.17) 

Remark 16. Since we will time integrate the previous quantities, the absence of a small 
parameter in front of E{t)'^ is not problematic. 

22.9. Asymptotic control of divu'*(^'^). We have 

Now, thanks to Lemma 21.2, this leads us to 
with 

llrilli < C\\Vn.{\\KVv\\.+Ki\\Vvh) < CE{tf. 
Next, we notice that 

{btUE]ri'i = mipj^ iiv^vm],'; 

with in virtue of Lemma 21.2, 

Ik2||| < C||Vr]'"||5(||KVw||5 +K.^||Vf;||2) < CE{tf. 
Now, since ibi)[v{0i),i = in (0, 1)^, this finally provides us with 

\\divu^{fj'')\\^<CE{tf. (22.18) 

23. Asymptotic regularity of KVt and of K^vt 

This Section is devoted to the asymptotic control of nvt and K'^Vt in spaces smoother 
than the natural regularity ffa (fi) for Vt, the idea still being that one degree in the power 
of K allows one more degree of space regularity. 

23.1. Asymptotic control of KVt in Hi (CI). Our starting point will be the fact that since 
g = on r, we have for any I e {1, K} on (0, 1) x {0}: 

"*°^' + d^IWW' 

where a;"*" = {—X2,xi). Therefore, we have on (0, 1) x {0}: 

{y/alvto0i),ni-v 01,1 = -[ V ° &i,n " V ° 



ai{0i)q 01,2 ~K a 1- ~« /3 



dctVr]'"(6'; 

o 0i,i-^^^ ■ fj'' o 01,1 



ai{0i)qo0i,2 - 



detVr?«(eO 
showing that 

K\\{y/mvto0i),ni-v'' o0i,i\\o g^o,i)^ < KE{tf + CK\\^i{0i)fj'' o0i,^^^^\\o^g^o -^y. (23.1) 
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By definition, 

K 

the sum being restricted to the indices i such that ^i((0, 1)^) and 0i{{O,l)^) have a non 
empty intersection. We then have 

K 

with = {67' o 6i)X' {67^ o {Or' ° ^0,^ (^r' ° Oi)X\ and 

«ll\/^(^i)^llo,a(o,i)^ < C'«;||i^||| sup [v^j? o 6ii] 113,(0,1)2 

+ Cfi;sup||pi ^'i] II 2,(0,1)2 

<C«(||17||| + 

< CKE{tf. (23.3) 

Now, we notice that for Xx G (0, 1) such that 6'j(a;i, 0) G ^i((0, 1)-^), we necessarily have, since 
for ah k e {l,...,K), efc([0,l] x {0}) = 90 n ^fe([0, 1]^), that o ^((a;i, 0) = (/i;(a;i),0), 
showing that on (0, 1) x {0}, wc have \/cn{Oi)a^^i^llf^ = except when ii = 12 = is = i4 = 1- 
Therefore, (23.2) can be expressed as 

K 

V^i{Oi)r o Oiain = V^i{0i)[J2 V^imp^ *h E%nn {6r' o e04!!!ii + A]. (23.4) 

i=l 

Now, from the properties of our convolution by layers, we have (since the derivatives are 
horizontal) that 

K||[pi *h -B**^],!!!! ||o,(0,l)x{0} < C'||-E'*'^,iii ||o,(o,i)x{o}- (23.5) 
Thus, with (23.2), (23.3) and (23.5), we infer 

K.\\^i{6i)f,'' o 61,1111 ||o,a(o,i)= < CKE{t) + CE{t), 
which coupled with (23.1) provides us with 

n\\{y/mvt o 6i),ni o 61,1 ||o,a(o,i)2 < CE{tf + C. 
This provides us the trace estimate: 

KWiv^ivt o ^0,1 {Oi' o irr'yr o oi,i (^r' ° {rr')h,af,^ie,ao,im 

< CE{tf + C. (23.6) 

Consequently with the divergence and curl estimates (22.14) and (22.15) and the trace 
estimate (23.6), we infer by elliptic regularity: 

«||(V^5to0i),i(^r'°(r)"')ll|,^»(e,((o,i)^)) <C^W'+C. (23.7) 

Remark 17. It is the presence o/||r^||9 in the inequalities leading to (23.3) which explains 

9 

the assumption of in H2. It is, however, not essential as will be shown in Section 26. 
One way to see this, is to smooth the initial domain by a convolution with the parameter k 
to form Vf^. Then, by the properties of the convolution, k||0'*|| 9 < C||0|| i . 
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23.2. Asymptotic control of n'^vt in H^{^l). In a similar way as in the previous subsec- 
tion, we would obtain the trace estimate: 

K'\\{y/mvt o eiu {ei' o irr'yr ° oi,i (^r' ° ('?'')-')ii3,a^'=(.,((o,i)^)) 

< CE{tf + C, 
which coupled with (22.16) and (22.17) provides 

K'||(v^^5toe;)a(er'°('7'')"')lli,^^(^K(o,in) <C^W' + C- (23-8) 

24. Basic energy law for the control of v and fjiK independently of k. 
We will use a different type of energy than in [6], namely: 

Definition 24.1. 

K 

^"W=oE/ ^ioOi\{voei),nif, 

2 ^(0,1)^ 

where = ^ ai, S, being defined in Section 20. 

Remark 18. The main differences with respect to the energy of [6] are in the absence in 
our energy of any restriction to the tangent components, allowing a more convenient set of 
estimates, and in a setting in Lagrangian variables. 

We have: 

K 

Ht{t) ^1° ^ivt o ei,ni V o ei,ni 

K 

= - V / 6 o 6i{{a'^)'iq,k ) o 61;, 111 Vj o 61,111 

1=1 Jio,iY 

K 

= -E / 6°^/[(^r)?9o^/,p],iii«,o^Mii, 

where 6f = [W{'n'^ o Oi)]-^ . Next, we see that = -[Hi + H2 + H^], with 

K 

Hi{t) = y^ / mmfpinqoOupVioeuiii, 

K 

Mt) = E / ° 0imrjQ°0i.piii ]vj ° ^i,iii , 



K 

trAo.i)^ 

We immediately have for the third term: 

\Hs{t)\ < CE{tf. (24.1) 
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Next, for H2, since {^iq) o Oi = on d{0, 1)^, 



K 

-E / [^i°Hbirj],pqoeiaiiVjoeian. 



We then notice that from the divergence condition, we have (6^^)^ o 0;,p= in (0, 1)^, 
implying 

K 



Ao,i)^ 

K 
K 

+ ^l°Oiqo 01,111 (bifj,i ij o Oupii 

7^1 J (OA)' 

K 



_ '(0,1) 

Now, in a way similar as (8.15), we have for any / e H^{{Q, 1)^) 

Uio6if,i llHi((o,i)^)'^^ll-^llHi((o,i)=)' 

since the derivative is in the horizontal direction. By applying this result to / = (6f )^,ii for 
the first integral appearing in the equality above, and by using the continuous embedding 
of i?^ into (6 e (1,00)), and of into (3 € (1,4)) for the other integrals, we then 
obtain: 

\H2{t)\<C\\a^\\.\\qUJ\vh<CEitr- (24.2) 

We now come to Hi, which will require more care, and will provide us with the regularity 
of ^(^(fi) in iJ 2 independently of k. We have: 

Hi{t) = / 6 o Giibifjaii [q°Oi],p vj o 61,111 



with 



= Hn{t) + Hi2{t) - I ^10 eiAP [q o ei],p vj o Oum , (24.3) 

JCofV(f?-ogO)^,ni ,, ^. . ^ 
detV(^^o^O 

Hi2it) = - / o .KCofv(,)« o o e^u o , 



and 



^ (CofV(^«o^ (CofV(^-ogO)^,ni [detV(ry-o^O],ni 

J' L detV(^«o^O J'^^' detV(^«o^O ^'^^i [detV(r?« o ^,)]2 ^ 
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SO that 

I / ii^iq) o ei],p Vj o Oiaii I < C\\vh. (24.4) 

^(0,1)2 

We now turn our attention to the other terms of (24.3), and to shorten notations, we will 
set: Qi = qo9i. We first study the perturbation H12, which would not appear if the volume 
preserving condition was respected by our smoothing by convolution. It turns out that we 
do need the double convolution by layers appearing in the definition of v** in order to identify 
time derivatives of space energies. We first notice that since 9i does not depend on t, we 
have: 

iroei)t = u^o{fj-oei), 

from which we infer in (0, 1)^, since 61 is volume preserving, 

[det(V?7'' o ei)]t = divu^ifj" o 61) detVifj" o 61). (24.5) 

24.1. Study of H12. We have after an integration by parts in time, and the use of (24.5): 

rt 3 



/ -ffl2 = / > -^^12 + ^ 
•^0 i=l 



with 






= 


fl 

Jo J(0,l)2 


. „.„^..^ pdiv{t'^(?7'* 6*0 [detV(77'" 06*0], 111/=, ~ ^ 
6 e,(CofV(ry 0,)r. ^^^^^^.^ ^ ^^^^^ ^, 9^,ni 




Jo J(0,l)2 


.(cofv(r . 


H12 = 


-/ ^ 


.o.Kcofv(,-o.o)-[-;:^:::^j- 


and 







\Ri2it)\<CtEitf + C. (24.6) 
24.1.1. Study of -ffj^j* ^'^^ the sake of conciseness, we denote 



We then see, by expanding the third space derivative of the determinant in the integrand 
of HI2, that = X^i=i HI2 + -R12 with the being estimated as HI2 that we precise 
below and RI2 being a remainder estimated as (24.6). By definition of f]'^ we have, if we 
denote 

Hll = f f [Aji{fj^ o 60,2 [[V^m pi i^h Er] {er^ o ei)] ,1111 

Jo 7(0,1)2 

[77, 0,(0-1 o 001,111 ] +RI2, 

with I-R12I < Ct and where, because of the term y/ai{6i), the only indexes i and I appearing 
in this sum are the ones for which 0;((0, 1)^) fl 0j((0, 1)^) ^ 0. Only such indexes will be 
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9 

considered later on when such terms arise. From our assumed regularity on in i?', we 
then have 

Hl^ = / [Mv5 ° [[pi *h Er] {Oi' o Oi)] ,nn 
Jo J(0,l)2 

with |i?{2l ^ CtE{t)^. We next have, since the charts 6i are volume preserving, 

= f I [[AAfi2 o eiUKer' « 0,) [pi ErUj.j.j, 

Jo Je-^{ei{o,i)^) 

fJij2j3j4 iii2i3\ /rrrf.. a a.] . . . 1 -L R 

'-if, nil '-a.iii Lv"»'o ° '^«J«iJ2J3 J T^-"-! 

with \R\ < CtE{tf and 

4tif * = ° ^O^Si {Oi' o eir,r {er' o ety\i {er' o e0^^l ](^r' ° ^0, (24.7a) 
47ii = [(^r' ° ^O'Si (^r' ° ^0'%i (^r' ° ^^)'^l ](^r' ° ^i)- (24.7b) 

Next, we notice that the term Aji {9'[^ o 9i) introduces a factor ai o 9^ which is non zero only 
if a; e 9~'^{9i{Q, if), leading us to 



Hll = f f [[Aji {f,2 o ] {9i' o 9,) p. ^^hnU 
Jo J<o,iy 



'(0,1)2 

Jij2j3j4 iii2i3\ /T^f.. ^ a.] . . . lip 
''i;,llll ''i;,lll LV"«'/J ° "jJmiJ2J3 JT^-"-' 

where 9^^o9i is extended outside of 9^^{9i (0, 1)^) in any fashion. This argument of replacing 
an integral on a subset of (0, 1)^ by an integral on (0, 1)^ will be implicitly repeated at other 
places later on. Now, since p is even, 

HU = f I E^,hh33u Pi ^ [[AjM ° ^0.2 ](^r' ° ^i) 

Jo ./(0,1)2 

4!^!t4:mW^n, o 9.U.,,,]+R. (24.8) 

Now, let us call / = [A,i{fi^ o 0(),2](^r' ° ^i)4'nn'4,m and g = [^m o 9iUi,i,. We 
notice that ||/|l|,(o,i)2 is the natural norm associated to E{t). Here we cannot use directly 

Lemma 21.2 for the case where all the ji = 3, since £^3333 is not necessarily in Hi{{Q, 1)^)' 
a priori. Instead, we write 

rVi 

f{yi,X2) = f{xuX2) + (yi - Xi).f,i {x-i,X2) + / [/,1 {X,X2) - f,l {xi,X2)]dx, 

J XI 

which shows that on (0, 1)^: 

Pi *h [f9]{xi,X2) =f{xi,X2)pL *h9{xi,X2) 

+ f,iixi,X2) / pj^ivi - xi){yi - xi)g{yi,X2)dyi 
Jr " 

+ / Pi(yi-2;i) / [f,i{x,X2) - f,i{xi,X2)]dx g{yi,X2)dyi. 
Jr " Jxi 
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This implies 

Jo ^(0,1)2 

<nn'4VmVi *h [V^^{e^)[m o e.Ui,,, ]] +R-R,- R^, (24.9) 

with 

Jo J(o,i)^ 

/ pi (2/1 - xi){yi - xi)g{yi,X2)dyi] dxidx2, 

Jr " 

R2= I I [-EJ^yijsjaj^ (a;i,a;2) 
Jo J(o,i)^ 

r ryi 

/ Px{yi-xi) / [f,iix,X2) - f,i{xi,X2)]dx g{yi,X2)dyi\ dxidx2. 
Jr ~ Jxi 

Now, for R2, we notice that since 

(a;,a;2) - /,i (a;i,a;2)| < C|/|5_((, i)2|a; - a;i|5 < C\fi''\7^\x - xi\^ , 

we have 

R2<C [ Wfj^W? I [\El'^,hj2j3j4{xi,X2)\ [ pi{yi- xi)Ki\g{yi,X2)\dyi]dxidx2 
Jo J(o,i)2 Jr " 

Jo J(o,i)^ 

<C \\r\\i^H\\Er,nn33jAmo,io,i)^+ \\{Eint,hhhu\\o,{OAr]\\9\\o,io,ir 
<CtKiE{t) + CtE(tf, (24.10) 

where we have used the fact that \/oriVio6i[,i^i^i^i^ ||o, (0,1)2 is contained in the definition 
of E{t). We now turn our attention to Ri. We first remark that 

/,i ={[AAv2 ° 0il2]{9i' o ei)d^^a ci^^^i^ 

4 
n=l 

np^„(^rio^,),f (0-io^,)], 

which imphes that 

4 

n=l 
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with 



Jo J(0,l)2 

/ pi{yi - xi){yi - xi)g{yi,X2)dyi] dxidx2, 



- a;i)g(2/i,a;2)(i2/i] dxidx2. 



Let us study i?}. If wo denote ft-(a:;i,a;2) = J^ pi{yi — xi){y\ — x\)g{y\,X2)dy\, since 
Jo J(o,\Y 

{[A,i{f,'ioei),2m'oe.y^%l)ah] 

Jo J(0,l)2 

Wi!^lt{[Aji{f,^2 o ei),2]{ei' o e^c^^a h]{6r^ o Oi)] . 

Since the derivative of {EY'-ij^jsii ° is in the horizontal direction, we infer similarly 
as in (8.15) that 

Since we have by interpolation i^(o,i)2 < i _(o,i)2, we then infer: 

\R\\ < CtE{tf. (24.11) 
In a similar fashion, for we can identify an horizontal derivative: 

ETmhhu o 6i)i- {ej' o 0,)],i n^=2(^r' ° {07' o Oi) 

which leads for the same reasons as for Rl to \RY\ < CtE{t)'^. Since the other i?^" are 
similar in structure, we have 

\R\"\<CtE{tf. (24.12) 
Consequently from (24.9), (24.10), (24.11) and (24.12), we infer 

Hl'2 = f f [[Aji (772 ° Oi)a ] {97' o 9i) i?ruu..3.4 
Jo J(0,l)2 

4!^iit<^p^ *h [v^mM o 9i],i,i,i, ]]+rii 
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with 

\rll{t)\ < CtE{tf + CtK^E{t). (24.13) 
Since EY", ,,,,,,,, c^^^ifA'' = 4an (^ru2j3j4 ° (^^r' ° ^*)' infer as for i?} that 
\Hl^{t)\ < Ct sup sup WAjiUm^^M^WI + \rll\ < CtE{tf + CtK^ E{t). 

3,1 [0,t] ' 

The other are estimated in the same fashion, leading us to 

\Hl2{t)\ < CtE{tf + CtniEit). (24.14) 

24.1.2. Study of H^2- Next, for iJ^2, we first notice from the asymptotic regularity result 
(22.18) on divit'*(ry"), that can be treated in the same fashion as HI2, leading to 

|i?i2(i)l < CtE{tf + CtniEit). (24.15) 

24.1.3. Study of H^2- We simply write 

^ / ^^flvrmr- ^^^ Jo[detV(^"°^i)]alldiv(^'-°g;) ^ . „ 

^7(0,1)^^'^^'^^^°^^' [dctV(r^'-ogO]^ "^'^ ^.-^^^(i), 

with i?f2 being bounded by a term similar as the right-hand side of (24.15). We also see that 
the first term of this equality can be estimated by a bound similar as the right-hand side 
of (24.15). The third term is treated in a way similar as H\2, in order to put a convolution 
in front of {fjj o 0(),iii. There is no difference linked to the fact that the integral from to 
t does not apply on all terms as for HI2 since pi and the 61 do not depend on time. The 
fourth term follows the same treatment as H12, leading us to 

\Hl2{t)\ < CtE{tf + CtniEit), 

which with (24.14) and (24.15) implies 

\Hi2{t)\ < CtE{tf + CtKiE{t). (24.16) 

24.2. Study of Hu. As for Hu, we have if we still denote = pi -kh. {{\/chv) ° 0%) and 
e"*" the sign of the permutation between {m,n) and (1,2): 



'(0,1)= 
'(0,1)2 

[vn o 6,{6^' o 60]aii_ 



■ / \f:,(ff.\ [g ° r piK] .... 



/(o,i)2- detV(^«o6';)(6'r^°^i 
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with 

\Rii{t)\<CE{tf, 

and 

Therefore, 

J(0,l)2 

with 

^aOi234 ^ r f;-"; ^1. (24.17) 

Similarly as in the study of (from equations (24.8) to (24.13)) we have: 



^(0,1)= 

+ 5ii + Rn, 

with, 

|Sii| <CtE{tf + CtK.^E{t). 
By integrating by parts in space (and using ^iq{6i) = on 9(0, 1)^), 

^^11 = -e™"e''" / (Ci9)(^i)/ifr''''Vi -^^h [y/cHVm ° pi -kh Walvn o 9i],i^i^i^u 

+ Hl^ + Hl^ + 5ii + (24.18) 

with 

>/(0,l)^ 

(24.19a) 

Hn = -e™"e" / g(6'i)[6(6'i)4r^'''']>n Pi [\/^^m o ei],j,j^j,j, pi i^h [s/^iVn o OiU^i^i^ 

(24.19b) 

For H^-i, by taking into account the symmetric role of {12,13,14} and {j2,j3,j4}, we obtain: 

7(0,1)2 

Next, since hl^J^'''' = h!i^^'^'\ this implies 



ilii — — e e 



/ (Ci9)(^i)/ilr ''^''*Pi [\/ai»?m O Qi],iihi2i3ii Pi [v^^n ° 

7(0,1)2 



(0,1)2 

)iiiii2i3i4 Pi *h [y/^^Vn ° Oijjj^j^j^ 



'(0,1)2 

Therefore, by relabeling sr as rs and as jz, we obtain by comparison to (24.19a): 
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and thus = 0. For Hfi, wc have by integrating by parts Hf^ — Hfl + Hfl, with 



J(0,1)^ 

J(0,l)2 

First, for HH we have if we denote Emn = ^"^"^ P^-^h['/^if]m°0i\,hi3u P^-^h[\/ori'iinO0i],i^i^i^. 

rt 1 rt r 

J'ii J'Ji it^mn 



Jo ^ Jo J{0,1)^ 

^ J(0,1)^ 



Therefore 



Jrti J'Ji ]*|| i_(o,i)2ll'?ll| 



l/ Hii\<c j \n\\[[m)[m)h[i:^^^^-] 

Jo Jo 

+cJ)\[[m)[m)hi^:^''^%Ah]t\\i,^^^^^^^ 

With the definition (24.17) and (22.18) for the control of the time derivative of det(V(r/'*)) 
m we then infer: 

1/ Hfl\<CtE{t)^ +N{uo). (24.20) 
Jo 

Next, for H^i, we have by relabeling m and n 

Hll = / qi0i)MOi)hi^:^''''U PL *h [Vmfin o 0i],nj,u [V^v^ o 0iUi,i,u 

= / q{0i)M0i)hif'^'%i^ pi -kh [V^iVn o 0i],j^jsu pi -^h [y/^iVm o 0i],nni3u ■ 

7(0,1)2 

By taking into account the symmetric role of {12,13,14} and {j2,j3,j4}, we then obtain: 



(24.21) 

Consequently, by (24.19b) and (24.21), 

2Hl^ = -e-"e" / q{0i) [^i {0i)h'^:^'''%, pi *h [s/^fim o 0i] pi *h [V^iVn o 0i] 

7(0,1)2 
7(0,1)2 



+ Hll 



+ Hll 



/ q{0mmh^^:^''nii [p 

7(0,1)2 



78 



D. COUTAND AND S. SHKOLLER 



Therefore 



Jo ^ JO ^(0,1)^ 

^ J{o,\Y ^ Jo 

Now, from (24.18) and Hyy = 0, we infer by integrating by parts in time: 

Jo ^ Jo J(o,\Y 

^ 7(0,1)2 

Now, we claim that the only couples (ii, ji) contributing to the sum above are the ones with 
i\ ^ j\. To see that, we notice that if i\ = ji, then by simply relabeling m and n, and using 
the symmetric role of {^25 ^3, ^4} and {72,^3574} 

hUi)i23ifrnn~m .... ... . _ u{ji) 1234 ^nm-n .... ... . 

"rs '/jre)«lj2j3j4 '/iK)JlJ2«3«4 "rs '/iK»lj2j3j4 '/iK»lJ2»3»4 

— fc(j»)i234,nm~n . . . . ^r? 

"rs '/lK»lJ2»3»4 '/lK)Jlj2j3j4 



leading to /iri'''^^*e'""^j",Jii2i3i4 ^j>t«ii2i3j4 = 0) since e™" = -e"™. Consequently if we 
denote 

(ii)l234 



^(ji)234 _ 



Jo ^ Jo -'(0,1)2 



we have 

rt 



^ J(0,l)2 



'(0,1)2 

ViKljlj2j3j4 ViKnii2i3i4 ]]q 



^*[5n+i?ii + ^^?fi]. 



Now, for any fixed (ii, ji), we have e''"e*iJi [(0^^^ o (0,"^ o 6*0,;^ ] = -det(V(6'^"^ o ^,)) 
— 1, leading us to 

_1 mn iiil /" / ffA5V6'-W(-'*)234l^m .... i^" ... . 

^ / / LVStyyV"*/" it'hK'JlJ2J3J4 'llKnil2t3M 



JO 



2 



J(0,l)2 

'iK'3i32J334 ^iK'iii2i3i4 ]n 

(0,1)2 



+ [ / m{0i)d^''^'''vr.m323334 n. 

^ "'(0,1)2 
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Next, by integrating by parts in space: 
rt 1 pt 



Jo ^ Jo -'(0,1)2 

'il ViK^3l32333i %KM2i3»4 Jo 

^ 7(0,1)2 

+ ^ [5ii+iiii + ^//i\], (24.22) 
where we have used the fact that similarly as for , we have 

n _ iiji j(ji)234^m . . . 

U C C U. '/JK'«1J1J2J3J4 '/lKn2J3l4 • 

We now come to the study of the crucial term bringing the regularity of the surface. 

24.2.1. Control of the trace of fjiK on T. Let us study the second term of the right-hand 
side of (24.22): 



'il ViK,:jlj23334 ViKTi-2i3U ' 

(0,1)2 



for which we have 



H = A^^n^i.n f [[[(ag)((?.)d(^')--] o er^ o (rf )-i o r o 9^^, 

^ 7(0,1)2 

= - Je-"e"^'i / [[[(^i9)(^i)'^^''^^'1 ° ° (^")"']nl ° r ° 
^ 7(0,1)2 

Now, for the same reason as before, the couples (i'l, ji) such i'-^ = j[ will not contribute to 
the sum above, leading us to 

^ 7(0,1)2 

[CU2J3J4 ° of,'' O Oi fif^,i2i3U ] 

^ 7(0,1)2 

[Cj2J3J4 ° O??" ° Ci2i3i4 ] 

= _ 1 mn/,,; / [[[(69)(^i)d(^'*)-1 o er^ o irr'U', 

^ 7rj-(e.((0,l)2)) 

[C,.2i3.4 09r' O (^«)-^]„-. Cn2i3i4 O^r' ° (^'^)"'] 

= /+ J, 
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with 



^ Jri'^(Oi((0,'i-)^)) 



^ ^^"(^^((O,!)^)) 



Next, we notice that 



with the perturbation term 

^ Ji)«(e.((o,i)^)) 
= Ji' + Ji, 

where 

MVi^mhu °0-' o (,7«)-i] C»2i3.4 °0r' o (^'')"'] • 

Now, for J|, let us set 
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SO that in order to identify an horizontal derivative on the highest order term we have: 



2^ 



2X 



fu ir ° oi),i {0^' o div[77i.„,,3 o irrx +ri 



= \( fa (div[%,„,,3 o^ri o (f?«)-i] o r o ^0,1 (^r' ° m-') + r\ 

with \r\\ < C||ry'*'||§. Now, we notice that the presence of the factor o {fj^")^^ in fu imphes 
that the integrand in the integral above is zero outside of ij'^{9i{{0,l)'^)). Similarly, the 
presence of pi, -k^ [^^fj o 9i],j.^j.^ o9^^ o {fi'^)~^ implies that x € fi'^{9i{{0, 1)^)) in order for 
this integrand to be non-zero. Therefore, 

Jl = l [ Mr ° ^0 (divfe„,,3 o^ri o o ^« o eOa deWir o 61) + rj. 

^ ^(0,1)2 

Now, since the derivative of div[f)iK;i2ia ° ir)^^] o fj'^ o 0i is taken in the horizontal 
direction, this implies: 

i-^i'i < 1,(0,1)2 II Mvi^,hh °^r' ° {rr']°r o ^iiii,(o,i)^ + i^li 

< C^ll/i^ II 1,(0,1)2 II div[,7..,...a o^r' o (r)-']||i,^»(e,((o,i)2)) + C||Id+ / vwi 

J 

Now, since we have in the same fashion as (22.8): 

II div[%«„,,3 o^ri o (r?'^)-^]|| i,^j.(,,((o,i)2)) < CtE{tf + CK^E{t) + C, 
we then have 

\Jl\ < CtE{tf + C + Cn^E{tf. (24.23) 

Next, for , we notice that fif {ry^],n is a sum of product, each one containing 

a factor (6'^^ o 6';),^" {9^^^ o {'q'^)~^)- This imphes that can be treated in the same way 
as Jl, with the identification of an horizontal derivative on the highest order term of the 
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integrand, leading to the same majorization. We can also treat / in a similar fashion, due 
to the curl estimate (similar as (22.11)): 

II curl[%„,,,3 o^ri o (r?")-^]|| i,^«(e.((o,i)2)) < CtE{tf + N{uo), 

which finally provides us with: 



with 



\h\t)\ < CtE{tf + N{uo) + (Ck5 + s)E{tf + Cs. 
Therefore, by integrating by parts. 



9i7'*(6li((0,l)2)) 



with 



/*)«(e.((o,i)^)) 
SO that, 

|/l'(<)|<[|kll||^||3 + ||^1|l|k||3] M\l + \h\t)\ 

< {Cni + 5)E{tf +Cs + CtE{tf + N{uo). 
Now, since q = Q and = a; on F, we infer 

H = -\l [a,((r)-')bo + fptU [Nm + r«)t] 

with d!lr = ^ o y Therefore, with the initial 

VdetV(77'= o Oi) W 

pressure condition 

Po,m iVm < -C < on r. 
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we infer 

-H<-C I \aii0i)d^''fll,j,j,j, C'Cn.i3i4 ]o07'o (fi-)-^ 

+ tE{tf + \h''\ 

+ {Cni + 5)E{tf + CtE{tY + Cs + N{uo) 
+ {Ck^ + 6)E{tf + CtE{tf + C5 + N{uo) 

JSf)«(0,((O,l)2)) 

+ (Ck5 + S)E{t)^ + CtE{tY + Cs + N{uo). (24.24) 

Now, it is clear that the space integral in front of the first time integral in (24.22) can be 
treated in a similar way, except that we do not have a control on the sign of the boundary 
term as in (24.24), which does not matter since a time integral is applied to it. This therefore 

leads us to 

- /* H,, <-C [ [ai(0OCiii Can ] o Or^ o (^ )-^ 

+ {Cni + S)E{tf + CtE{tf + + N{uo), 

which , with (24.1) and (24.2), finally gives the trace control for each as well as the 

control of v around F: 

H^{t)+c f i[v^(w«,iii]o0r'°r)"T 

Jar)''(ei((o,i)2)) 

< 5E{tf + CstE{tY + CsN{uq). (24.25) 

24.3. Asymptotic regularity of each fjiK- Consequently, we infer that for each I e 
{1, ...,K}, we have the trace control 

||[x/^(^z)^i.,i]oer'°(r)-iWoe,((04n < SEitf + CstE{t)^ + CsNiuo). (24.26) 

Consequently, with the estimates (22.8) and (22.11) on the divergence and curl, we obtain 
by elliptic regularity: 

\\[v^i{ei)fii.a]o o (r)-i|,^«oe,((o,i)=) ^ p{\mi)[mtf + cstm^ + ^^^^(^^0)] 

< 5E{tf + CstEitf + C5N{uo). 

Therefore, 

\\^i{ei)ftu.i |||, (0,1)2 < 5E{tf + CstE{f:f + CsNiua), (24.27) 

which implies that \\^/cn{0i)fiiK.,i2 III, (0,1)2 and \\y/cn{9i)riiK,2 |li,a(o,i)2 are controlled by the 
same right-hand side as in (24.27). Consequently, with (22.8) and (22.11), we infer in the 
same way as we obtained (24.27) from (24.26) that 

\\V^i{0i)fiu,2 |||,(o,i)2 < SE{tf + CstE{tf + CsN{uo), 



84 D. COUTAND AND S. SHKOLLER 

and finally that 

IIVc^('^;)??i«!l|,(o,i)2 < SE{tf + CstEitf + CsN{ua). (24.28) 

24.4. Asymptotic regularity of f]'^. This also obviously implies that for the advected 
domain 

< SE{t)^ + CstE{tf + QiVK), (24.29) 

where Vl^ = f\i^K+\ (suppaj)'^. 

Prom the divergence and curl estimates (22.4) and (22.10), we then infer that 

WPV o ^ill|,(o,i)^ ^ '^^W + Csi^{tf + CsNiuo), (24.30) 
which with (24.28) provides 

lini' <SE{tf + CstE{tf + CsN{uo,K\\n\\9). (24.31) 

2 2 

24.5. Asymptotic regularity of v. The relation (24.25) provides us with the asymptotic 
regularity of v near dQ. For the interior regularity, we notice that if we time-differentiate 
the analog of (22.5) for the cut-off /3 e we obtain 

\\divii(3io9i),.oei'o{r)-')\\^^^^^^^<C. (24.32) 

Similarly, we also have 

\\cmU(f3vo9i),,oer'o{r)-')\\^^-^^^^<C. (24.33) 
From (24.32) and (24.33), elliptic regularity yields 

11/3^0 0,112,(0,1)2 <C, 
which, together with (24.25), provides 

\\v\\l < 5E{tf + C5tE{tf + CsN{uo). (24.34) 

24.6. Asymptotic regularity of q. Prom the elliptic system: 

g = on r, 

we then infer on (0,Tk) (ensuring that (19.2) is satisfied): 

||g|| 7 < CWfj^fv < SEitf + C5tE{tf + CsN{uo). (24.35) 

2 2 

24.7. Asymptotic regularity of n^/oiv o Oi. From (23.7), we have: 

[kW^iv o 0,112,(0,1)2]' < [«||uo|||]' + SE{tf + CstE{tf + CsN{uo). (24.36) 

24.8. Asymptotic regularity of niy/aiv o 9i. From (23.8), we have: 

[K''\\^lVo6l^^^^^^yf < [i^WuoW^f + CstE{tY + CsN{uo), 
which by interpolation leads to 

[wi 11^575 o0,||4,(o,i)2]2 < [Ki\\uoW + 5E{tf + C5tE{tf + CiN{uo). (24.37) 
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24.9. Asymptotic regularity of qt. From the elliptic system: 

g = on r, 

we then infer on (0,T„): 

\\qt\\l < SE{tf + CstE{tf + CsN{uo). (24.38) 

24.10. Asymptotic regularity of Vf Since vl = —{d'^y-q,j, we then infer on (0,T„): 

ll^tlll < C{\\r\\i + \\q\\if < SE{tf + C5tE{tf + CsN{uo). (24.39) 

2 2 2 

25. Time of existence independent of k and solution to the limit problem 

By (24.31), (24.35), (24.34), (24.30), (24.28), (24.36), (24.37) we then infer the control 
on (0,T«): 

E{tf < 5E{tf + CstE{t)^ + CsN{uo), 
which for a choice of 5q small enough provides us with 

Eitf < Cs,N{uo) + Cs,tEit)\ 

Similarly as in Section 9 of [8] , this provides us with a time of existence = Ti independent 
of K and an estimate on (0,Ti) independent of k of the type: 

E{tf < No{uo), 

as long as the conditions (19.2) hold. Now, since ||??(t)||3 < ||Id||3 + / II^^Hs) we see that 

1 

condition (19.2b) will be satisfied for t < — -. The other conditions in (19.2) are satisfied 

No{uo) 

with similar arguments ((24.38) and (24.39) are used for (19.2c) and (19. 2d)). This leads 
us to a time of existence T2 > independent of k for which we have the estimate on (0, T) 

E{tf < Noiuo), 

which provides by weak convergence the existence of a solution {v,q) of (1.1), with a = 0, 
on (0,T). 

26. Optimal regularity 

In this section, we assume that 51 is of class Hi in K'^, that uq G H^{^1), and that the 
pressure condition is satisfied. We denote by N{uo) a generic constant depending on ||wo||3- 
With these requirements, we will only get the regularity of the moving domain r]{Q) 
and not of the mapping 77. 

Due to the fact that if 2 is not continuously embedded in L°° in the case that SI is 
three-dimensional, we cannot directly study the integral terms as in Section 24 as we did 
for the two-dimensional case. Instead, we are forced to also regularize the initial domain, 
by a standard convolution, with a parameter e > fixed independently of k, on the charts 
defining it locally, so that the initial regularized domain 0,^ = 0, obtained in this fashion is of 
class C°° . The regularized initial velocity, by a standard convolution, will be denoted ?io(e). 
We then start at Section 18 in the same way except that the regularity of the functional 
framework is increased by one degree for each quantity. This leaves us with the existence 
of a solution to (17.1) on (O.T^.t), with initial domain fig and initial velocity ?io(e). We 
then perform the same asymptotic analysis as k ^ as we did in Sections 19 to 24, in this 
new framework. We then see that the problematic term is now updated to one which can 
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be treated directly by the Sobolev embedding of H"^ into L°° in 3d. This leads us to the 
existence of a solution to a system similar to (1.1) (with a = 0) with initial domain Qe on 
(OjTe), with S L°°{0,Te; (ilc)), with initial domain fie and initial velocity Mo(e)- 

We then study hereafter the asymptotic behavior of this solution and of as e ^ 0. 
This will be less problematic than in Section 24 since the convolutions by layers with the 
parameter n do not appear in the problem (1.1) with smoothed initial data and domain. 
We will denote the dependence on e this time by a tilde, v standing here for for instance, 
and prove that as e ^ 0, the time of existence and norms of v are e— independent, which 
leads to the existence of a solution with optimal regularity on the initial data, as stated in 
Theorem 1.4. 

Our functional framework will be different than in Sections 19 to 24. Our continuous in 
time energy will be: 

Definition 26.1. 

H{t) =sup[|n|2 f + ll^^lls + 11^5*11 5 + ||g||3 + Wvtth] + 1, (26.1) 

[o,t] 

where h denotes the unit exterior normal to fi{fl). 
Our condition on will be that on (0, T^), 
1 3 ~ 

- < detV?? <-mfl, (26.2a) 

||'7l|3<|f^| + l, i|g||3< IkolU + l, < ||«o||| + l, (26.2b) 

ll^tlb < Ikilb + l, (26.2c) 

yie{l,...,K}, \fio0i,^xfio9i,2\> ^\eiaxei,2 \ on (0,1)^ x {O}, (26.2d) 

where Wi = —Vqo G i?i(f2). We will use a more straightforward approach than in Section 
24, which is enabled by the fact that we have a instead of the convolution by layers a** in 
our equation, by defining the following energy: 

Definition 26.2. 

K 

^'(*) = E/ ^oei\D\vuoei)f, 

1=1 J(OA)' 

where D^f stands for any second space derivative in a horizontal direction, i.e., f,aio:2! 
where ai G {1,2}. Summation over all horizontal derivatives is taken in the expression for 

Remark 19. We also note that this energy is associated with the second time- differentiated 
problem; we thus avoid the use of the curl relation (22.11) for i), which necessitates the 
supplementary condition curluo € H^{Vt) (which we do not have here). 
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With &i = [V(?7 o we have: El = with 



E,{t) = 


1=1 


(0,1)3 


^;)^ 


E2{t) = 


K 

1=1 


/ m)D{[{bi)j]u)D{qo9i),kD^itt 




Esit) = 


K 
1=1 


(0,1)^ 




Ei{t) = 


K 

(=1 


/ i{9i)D[(bi)';]tD{qt o 9i),k D\iu o 

J (0,1)^ 




Esit) = 


K 

1=1 


i(0,l)3 




Eeit) = 


K 

-2E 
1=1 


■/(0,1)3 




Er{t) = 


K 

1=1 •' 


(0,1)^ 




Es{t) = 


K 

-2E 

1=1 


■/(0,1)3 


Y, 


Eg{t) = 


K 

1=1 •' 


(0,1)3 





26.1. Estimate for ijt, q^t and q^tf From the eUiptic system 

a'i{0'iqt,k),j = -[ai{ai]tq,k),j -[a^5',fc a^ii-',; ]t in Q, 
= on 9f2, 

we infer 

\\qth < C[ Ms + \\qh + Wvh + Wvth] < CH{t). (26.3) 
For similar reasons, we also have 

\\qn\\i<CH{t), (26.4a) 
Mmh < CH{t). (26.4b) 
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26.2. Estimate for E2, E4, E5, Ee, E%. We first immediately have thanlcs to the embed- 
ding of into and into L^: 



\E2{t)\ 


< C\\&uh\\q\\^\\mh < C[ ||z)t||2||77||3 + Ml] \\qh\\i>tth < CH{t)\ 


(26.5a) 


\EA{t)\ 


<C\\~ath\\qt\\'.\\vuh<CH{tf, 


(26.5b) 


\E,{t)\ 


<C\\vU\qtU\iuh<CH{tf, 


(26.5c) 


\E^{t)\ 


<C\\ih\\qth\\iu\\2<CH{tf, 


(26. 5d) 


\Es{t)\ 


<q|a||2||g„||5||i}«||2<Cif(<f, 


(26.5e) 



where we have used (26.3) for (26.5c), (26.5d), and (26.4a) for (26.5e). 

26.3. Estimate for E^. By integrating by parts, and using [(6)^],^= 0, we obtain E^ = 
El + El, with 

El = T.I mm%tD\qoei)D\vuo0i),i, 

K 

El = Y. ^(^')''^ m%uD\qoei)D\vuoeiy. 

1=1 

We first have 



El{t)\<C\\au\\i\\qh\\vuh<CH{tr. 



Next, = E?=i ^1% with 



^" = -E / m)Dm'j]uD^qoei)D{d,,o0,),i, 

^3' = -E / Dm)m';]uD'{qoOi)D{iuo0t),i. 

We obviously have |£^3^(t)| < CH{t)'^. Next, we have by integrating by parts in time: 



Jo i(0,l)3 

+ E / mm'j]uD'{qoeOD{i,oeO,i]l, 

1=1 -'(oAr 
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showing, with the continuous embedding of into and of ff ' into L^: 

I I < Ctsup[||a«||i||gt||3||i)t||| + ||a«t||i||gl3||t5t|||] 
Jo [0,t] 

< Ctsnp[ra„\U\Ms\\vt\\i + ||am||i|kl3||i5*|||] 

[0,t] 

+ ||{)«|5 ||g(0)|l3||5„(0)||i + \\it\.tsnp[\\qth\\btt\\i + \\qUbm\\i] + NM 

[o,t] 

<CSH{tf +tH{tf + CsN{uo), (26.6) 
for any S > 0. For the remaining term E^^, 

\El'\<C\\ht\\0h\\vuh<CH{tr. 

Consequently, we have 

I / Es\< C5H{tf + tH{tf + CsN{uq). (26.7) 
Jo 

26.4. Estimate for E-j. By integrating by parts, E7 = E^ + E^, with 

K 



ir[J{o,i)^ 
^ r 



'(0,1) 

We first have 

\EUt)\<C\\aUquh\\itth<CH{t)\ 
Next, we notice by integrating by parts in time and space that 

Jo l^Jo -^(0,1)3 

K 

+ [E / m)D'{b)'^mt o eo,k D^vt o eiy]l 
1^Ao,i)3 

In the same fashion as we obtained (26.6), we then infer 
I Ti^^l <Cisup[||6,||2||gtt||2||^t||| + ||6||2||gt«||2||*t|||] 

Jo [0,t] 

+ ll«tll|llfe(0)||2||6(0)||2 + \\vt\\itsnp[\\q„t\\2\\bh + Il9tt||2||6t||2] + N{uo) 

[0,t] 

< C5H{tf + tH{tf + CiiV(Mo), (26.8) 
where we have used (26.4b) for qut- Consequently, we have 

I / E'j\< C5H{tf + tH{tf + CsN{uq). (26.9) 

^0 
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26.5. Estimate for £^9. We notice by integrating by parts in space that 

K 



m) = E / m)[(h)'}]D\qu o ei)D\vu o 6i),i 

K 



1=1 -^(0,1) 

Next by the divergence condition, 

K 



E9{t) = - E / m)D'[{bi)']]D'{qtt o ei){vu o ei),i 

K 

- 2 E / m)Dm^]D^{qtt o 9i)D{vu o ^04 

J (.0,1)3 

K 

7-:lJ (0.1)3 



1 = 1 -^(0,1)3 

showing that 

\E,{t)\ < CfbhWquhJitth < CH{tY. (26.10) 

26.6. Estimate for Ex. If e^^'' denotes the sign of the permutation between {i,j,k} and 
{1,2,3}, if i,j,k are distinct, and is set to zero otherwise, we obtain 

El = E\ + E1, 

with 



= E / ^(^'),fe 9 ° eiD%h)%D\vtt o eiY, 
ir[J(o,i)3 

1 ^ f 

2 i(0,l)3 

26.6.1. Estimate of Now, for El, since 



AAti + u,j -V[u,^ u,j ] = 0, 

we obtain in O that 

o^Mv,k ),j {t) = A«(0) + / K{a'li,i ),m]l=i - [ a^i,l af{d[v,i ),„ , 

Jo Jo 

and thus, 

ai{a^vt,k ),j = -[ai{a'l]tv,k ),j +[<(af 5,1 a^jV,i ),m ]l=i - a']v,i df{alv,i ),„ . (26.11) 

By elUptic regularity in the interior of f2, we infer that for any iv whose closure is contained 
in O, 

Wvthu, < C^[\\di{a^vt,k),j + ll^^tlb] < C^m. 
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With this estimate and the condition ^ o 6'j,fe = in a neighborhood of (0, 1)^ x {0}, we then 
obtain 

|i;J|<q|gl2[c#(t)||77ll3 + ||i5||^]||5«||2 

< C{H{tf + 1). (26.12) 

26.6.2. Estimate for and the trace regularity. We now study E^, which will be the 
term bringing the asymptotic regularity of the moving domain ^(f2). We have that 

3 

1=1 

with 

K 



K 

Ef = 2 V e'""^ef'''= / £,m)[D{vt o 6i),^ D{fj o 61),^ ]D\vu o Oi^ 

K 

= E l^e'"™^^^'' / m)[D\f,„ o ei),^ D\f,u o 61)^ ]tq oei{fjo 61),^^ . 

r-r2 J(o,i)3 



1=1 

We first notice that 

;i21| I \-u'22\ ^ ri\\x\\ ll„-.l|2||r-. 



| + |i?f|<qkl3||^||^||^«||2 + q|«||3||^t||j||^||3||i)«||2 

< CH{tf. (26.13) 
Now, for the remaining term E^^ , an integration by parts in time provides 



[ j-,23 l p231 I -'- rp232it 
Jq 1 ^ 2^ T 



with 

K ,.t . 

Jo J(0,l)3 

^232 ^ ^mni^pgj f ^(ei)D\vt O Ol),^ D\vt O Oi),) q O 6l {f, O Ol), 



■"9 

'(0,1)3 

First, for the perturbation term E^^^, by integrating by parts in space (and using g = on 

r): 

Jo ^ J(0,l)3 



+ 



Jo -^(0,1)3 



For the first integral, we notice that for any /, g smooth, 

mni pqj f m fi^n mni jqp f m fi^n 
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and, since e^"^^ = —e^''^, this quantity equals zero, leading to 

Jo -^(0,1)3 

Consequently, as in (8.15) since the derivatives in are horizontal, 

\Ef'\ <C f \\vt\\l\\[ ||V(5*Vr?)||i + \\V{qVv)h] 
Jo ^ 

< CtH{tY. (26.14) 

Now for we will introduce the notation 

V{l)=vo6i and E{1) =7)061. 

We then have after a change of variables made in order to get vector fields whose divergence 
and curl are controlled: 

K 



^232 ^ ^mni^^j f „ (^2^^(^) „ E{1)-'),^^ iE{l))E{l) 

1^1 -^(0.1)^ 



Pi 



(7^^y,(0 o E{l)-%^ {E{l))E{l),f m) o E{l)-'),l E{l),l^ \ . 

Now, we notice that any triplet (ii, ji, qi) such that Card{zi, ji, gi} < 3 will not contribute 
to this sum. For instance, if ji = qi, we notice that by relabeling j and q, 



eP''^E{l),l^ E{l),f E{l),il = ^""mVp E{l),i' E{l),f , 
where ji = qi is fixed in the sum above. Now, since e^'-' = —e^^'', this shows that 

ePijE{l)P^ E{l),f E{1),{- = 0. 

By a similar argument, 

ePiiE{l),l- E{l)P/ E{l),f = 0. 
Consequently, only the triplets where Card{ii, ji, (71} = 3 contribute to Ef^^, showing that 

K 

^232 ^ gmm^pi9iii / [^^ „ 9i eP«ePi9Ui£;(Z) Pi E{l),f E{l),l^ 
1=1 ■^(0,1)=' 

{D^mi)oE{i)-');;^ {E{i)){D^Vt{i)oE{i)-%^ [Emu 

Since for each given (pi, ji,gi) we have cPi^p^"^^^ E{1)p^ E{l),f E{l),f = detyE{l) = 1, we 
then infer 

K 

i^232^^e™"V« / _ [^qofi-\D'Vt{l)oE{l)-^),;^{D'Vt{l)oE{ir%6^]. 
By integrating by parts in space, we get by using g = on F: 

K 

£232 ^ _Y^^ram^p,j / ^ o r\D^Vt{l) o E{1)-'):;) {D^Vt{l) o E{l)-^y6^] 

Jfi{h) 

K 

-E^"""^"" / [(?9~or')u {D^Vt{l)oE{l)-^),^ {D^Vt{l)oE{l)-^y5^]. 
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Next, since for any / smooth, 

mni pqj f m _ n 

we infer 

K 

Now, we notice that for ^^e'^^^e^'^-' ^ 0, if p = i, then necessarily j = m. Similarly, if pj^ i, 
then since p ^ n, necessarily, p = m, and thus i = j. Therefore, 

K 

Ef^ = - V e^-e™™ / {D^Vt{l) o £;(/)-^),r {D^VS) o E{l)-y {^q ofj-'),m 

Now, in the same way as we obtained the divergence and curl estimates (24.32) and (24.33), 
we have the same type of estimates for Vt{l) leading us to 

w^m) [(i)Vt(Oo£;(o-i),r-(^'^KO°i?(o-'),™]ii^i(^(,-,)), <c^, 

||V?(0Odiv(i?V,(Ooi?(/)-i)||^i^.^-^^, < c. 

The fact that contains horizontal derivatives once again played a crucial role in these 
estimates. This implies 



1=1 i^m-'fi(^) 

+ E / {D'mi) o E{1)-')^ {D^VS) o E{l)-^f i^q o fj-'h +Ri 

1=1 Jfim 



K 



= E / iD'Vtil)oE{l)-'),liD'Vt{l)oE{l)-y {^qor'),m+Ri, 

with \Riit)\ < CtH{t). Consequently, 

Ef' = -\Y. I \D'V,{l)oE{l)-'\'A{^qor') 
1 ^ f 

+ oE / \D''Vt{l)oE{l)-^Wiqor^),^hm + Ri. 

^ ^ Jdflill) 

If we write q = q{0) + qt, n = N + ht, and use the fact that ^ = 1 on f , we then get 
Jo Jo 

1 ^ r 

^i'' = oE/ \D^Vtil)oE{ir'\^ m,mN^ + R2, (26.15) 
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with \R2{t)\ < 6H{tf + CstH{t)'^ + CsN{uo). Together with (26.5), (26.7), (26.9), (26.10), 
(26.12), (26.13) and (26.14), this provides us on [0,T«] with 

< 5H{tf + CstHitf + CsN{uo). (26.16) 

Similarly as in Section 24, from the pressure condition (1.2), this provides us with an estimate 
of the type: 

ll^^ttll^ + \\vt\\l < 5H{tf + CstHitf + CsN{uo). (26.17) 

2 

The control 

ml < 5H{tf + CstH{tf + CsN{uo), (26.18) 

is then easy to achieve by elliptic regularity on the pressure system. 

Next, we see that (26.17) implies that for any / e {1, ...,K}, vt°Oi, and thus {qo9i),-ifjo 
Oi,i XT] o 01,2, are controlled in H'^{{0,lf x {0}) by the right-hand side of (26.17). This 
implies the same control on {0,T,) for g ° (^0.3 ° xryog;,2 ^2((o, 1)2 x {0}), i.e. 

|go {dl),3.fjo6ui X?7o6';,2 I 

that 

\h\l^^ < 5H{tf + CstH{tf + CsN{uo), (26.19) 

which brings the Hi regularity of the domain fji^)- 

Now, for -5, we notice that from the identity on (0, 1)^ x {0}: 

Vttil) +qoei,s y(/),i xE{l),2 +qoei.:, E{1),, xV{l),2 +qt o Oi,3 E{l)a xE{l),2 = 0, 

we infer by taking the scalar product of the above vector by E{l),i that 

l^(Oa -^ll.co.D^xio} ^ ^^{tf + CstH{tf + CsN{uo), 

which by divergence and curl relations for £,{Oi)V{l),i {E{1)~^) similar in spirit to the ones 
in Section 22, leads to 

\\m)V{l),i 11^,(0,1)3 < 5H{tf + CstH{t)^ + CsN{uo). 

In a similar fashion, 

\\m)V{l),2 ||2,(o,i)^' < ^H(t)^ + CstH{tf + CsN{uo). 

Now, with divergence and curl relations for £^{0i)V {1){E{1)~^) similar in spirit to the ones 
in Section 22, this leads to 

II^SII^ < SH{tf + CstH{tY + CsN{uo), 

and consequently, with the control of the divergence and curl of v inside Q as in Section 22, 
we get 

||w||3 < SH{t)^ + CstH{tf + CsN{uo). (26.20) 

Now, with the estimates (26.17), (26.18), (26.19) and (26.20), we then get similarly as in 
Section 25 the existence of a time T > independent of e such that on (0, T) the estimates 
(26.2) hold, and such that we have H{t) < N{uo) on (0,T) for any e > small enough. 
Therefore, we have a solution to the problem with optimal regularity on the initial data and 
domain as the weak limit as e — > 0. 
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27. Uniqueness 

Let {v, q) and {v, q) be solutions of 1.2 on [0, T]. We denote 5v — v — v and Sq = q — q. 
We then introduce the energy: 

K 

fit) = E / ^(^OI^'K o Oi - v„ o 9i)\^ 
1=1 JiOAr 

where D^v stands for any second order horizontal space derivative v,i^i2. By proceeding in 
the same way as in the previous section, and using the fact that the divergence and curl of 
6v have a transport type structure as well, we obtain an energy inequality similar to (26.17), 
without the presence of N{uo) (since 6v{0) = 0). This establishes uniqueness of solutions. 
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